Ї 


ж 
х * 
* * 
* * 
in 
% * WYDZIAŁ CHEMII 
* « 
* 


UNIWERSYTET MARII CURIE-SKŁODOWSKIEJ 
W LUBLINIE 


Wydział Chemii 
Instytut Nauk Chemicznych 


mgr Magdalena Szaniawska 


BADANIE WŁAŚCIWOŚCI ADSORPCYJNYCH 
I OBJĘTOŚCIOWYCH NIEJONOWYCH 
SURFAKTANTÓW TYPU KOLLIPHOR I ICH 
MIESZANIN Z BERBERYNĄ 


Rozprawa doktorska 
wykonana w Katedrze Zjawisk Międzyfazowych 
pod kierunkiem dr hab. Katarzyny Szymczyk, prof. UMCS 


Lublin 2023 


„Trzeba mieć wytrwałość i wiarę w siebie. Trzeba wierzyć, 
że człowiek jest do czegoś zdolny 


i osiągnąć to za wszelką cenę.” 
Maria Skłodowska-Curie 


Składam podziękowania 

Pani Promotor dr hab. Katarzynie Szymczyk, prof. UMCS 

za cenne uwagi, pomoc w przygotowaniu rozprawy doktorskiej, 
poświęcony czas oraz wyrozumiałość. 

Dziękuję mojej Siostrze Sylwii oraz Przyjaciołom 

za wiarę w moje możliwości nawet wtedy, gdy sama w nie wątpiłam. 


Bez Nich nie osiągnęłabym tego. 


Рапиес! Moich Rodzicow 


Janiny i Arkadiusza 


SPIS TRESCI 
t STRESZCZ N E о DIETS pH das baa ase Gs M RA AER E 5 
2. STRESZCZENIE W JEZYKU ANGIELSKIM ........... eR emere 6 


3. WYKAZ PUBLIKACJI NAUKOWYCH STANOWIĄCYCH PODSTAWĘ 


POSTĘPOWANIA DOKTORSKIBGO sed ka oditaut 7 
4. WYKAZ STOSOWANYCH: SYMBOLI. ian Ай 9 
5. WPROWADZENIE I CEL ROZPRAWY ООКТОҚКЅКІЕЈ............................. 11 
6. OMÓWIENIE WYNIKÓW Gawin e t t es e e 19 

6.1. Właściwości fizykochemiczne ELP, RH40 i Вег............................... 19 


6.2. Właściwości adsorpcyjne i objętościowe dwuskładnikowych mieszanin 


ЕГЧ ага УРИ ЕНИР AARONA AEC 23 


6.3. Właściwości adsorpcyjne i objętościowe ELP + Ber, RH40 + Ber 


12 БЭК 20 Ве cod tee cases A OLA 27 
ZNANIOŚKI zo NI E EE E EE E E A 35 
8- LITERATURA AC ORO Z E E A EG 37 
9. PODSUMOWANIE DOROBKU МАЮКОМЕСО....................................... 48 

ЭГ Püblikacge naukowe: хулан O A A o i steh fed ead 48 

9.2. Konferencje naukowe. iR O LE EE 51 

9.3. Projekty badawcze, -KONKUNSY eat A EE KE r 58 

9.4. Aktywność organizacyjna i członkostwo naukowe .......................... 59 

9.5. Nagrody i wyróżnienia. Was AiR dO 60 
ТОДА АС2ГШЕГ ZO О о АО erat ye 61 


1. STRESZCZENIE 


Tematyka rozprawy doktorskiej obejmuje badanie właściwości adsorpcyjnych 
i objętościowych dwóch surfaktantów niejonowych z grupy Kolliphor: Kolliphoru 
ELP (ELP) i Kolliphoru RH40 (RH40), a także ich dwuskładnikowych mieszanin 
o różnym ułamku molowym w środowisku wodnym bez i w obecności 
naturalnego alkaloidu berberyny (Ber) w szerokim zakresie stężenia surfaktantów 
i różnych temperaturach. Zastosowany do badań alkaloid, berberyna, to 
substancja, która od wielu lat cieszy się zainteresowaniem naukowym ze względu 
na swoją bardzo szeroką aktywność biologiczną. Jednym z głównych problemów 
związanych z jej zastosowaniem jest słaba rozpuszczalność w wodzie, a więc także 
biodostępność. W związku z tym poszukuje się sposobów zwiększenia 
biodostępności berberyny np. przy użyciu surfaktantów. 

W przedstawionej rozprawie doktorskiej za pomocą pomiarów kąta zwilżania 
ELP i RH40 na powierzchni modelowych ciał stałych oraz kątów zwilżania wody, 
formamidu i dijodometanu na sprasowanej berberynie wyznaczono składowe 
i parametry napięcia powierzchniowego ELP, RH40 i Ber oraz przeanalizowano je 
w aspekcie adsorpcji na granicy faz woda-powietrze. Następnie określono 
właściwości adsorpcyjne i objętościowe dwuskładnikowych mieszanin ELP 
i RH40 o ułamku molowym w fazie objętościowej a równym 0,2; 0,4; 0,6; i 0,8 celem 
wyboru do dalszych badań z Ber mieszaniny surfaktantów, dla której istnieje efekt 
synergetyczny w redukcji napięcia powierzchniowego i procesie tworzenia micel. 
Końcowym etapem badań była ścisła charakterystyka fizykochemiczna mieszanin 
Ber z ELP, RH40, a także ELP+RH40 o ułamku molowym ELP w fazie 
objętościowej równym 0,8 w szerokim zakresie stężeń  surfaktantów 
(1.105 - 1-102 mol/ dm3) i temperatur (293-313 К). Na podstawie otrzymanych 
wyników określono zachowanie się berberyny w środowisku wodnym bez 
i w obecności surfaktantów. Stan równowagi pomiędzy roztworem 
monomerycznym a micelarnym surfaktantów, ich dwuskładnikowych mieszanin, 
a także mieszanin z Ber poddano ścisłej analizie termodynamicznej w oparciu 


o otrzymane rezultaty. 


2. STRESZCZENIE W JĘZYKU ANGIELSKIM 


The subject of the doctoral dissertation includes the study of adsorption and 
volumetric properties of two nonionic surfactants from the Kolliphor group: 
Kolliphor ELP (ELP) and Kolliphor RH40 (КН40), as well as their binary mixtures 
with different mole fractions in aqueous environment without and with the 
presence of natural alkaloid - berberine (Ber) in a wide range of surfactant 
concentrations and different temperatures. The alkaloid used for research, 
berberine, is a substance that has been of scientific interest for many years due to its 
very broad biological activity. One of the main problems associated with its use is 
poor solubility in water, and thus also bioavailability. Therefore, ways to increase 


the bioavailability of berberine are sought, e.g. using surfactants. 


In the presented doctoral dissertation, using the measurements of the contact 
angle ELP and RH40 on the surface of model solids and the contact angles of water, 
formamide and diiodomethane on compressed berberine, the components and 
parameters of the surface tension of ELP, RH40 and Ber were determined and 
analyzed in terms of adsorption at the interface between water and air. Then, the 
adsorption and volumetric properties of binary mixtures of ELP and RH40 with 
a mole fraction in the volumetric phase equal to 0.2; 0.4; 0.6 and 0.8 were determined 
in order to select a mixture of surfactants for further research with Ber, for which 
there is a synergistic effect in the reduction of surface tension and the process of 
micelle formation. The final stage of the research was the direct physicochemical 
characterization of mixtures of Ber with ELP, RH40, and ELP + RH40, where the 
ELP mole fraction in the volumetric phase was equal to 0.8, in a wide range of 
surfactant concentrations (1-109 - 1.10? mol/ dm?) and temperatures (293-313 K). 
Based on the obtained results, the behavior of berberine in the aqueous environment 
without and in the presence of surfactants was determined. The equilibrium state 
between the monomeric and micellar solutions of surfactants, their binary mixtures, 
as well as mixtures with Ber subjected to thermodynamic analysis based on the 


obtained results. 


3. WYKAZ PUBLIKACJI NAUKOWYCH STANOWIĄCYCH 
PODSTAWĘ POSTĘPOWANIA DOKTORSKIEGO 


Cykl powiązanych tematycznie 4 artykułów naukowych wchodzących w skład 
rozprawy doktorskiej, opublikowanych w latach 2019-2023 w recenzowanych 
czasopismach naukowych, znajdujących się w aktualnym wykazie czasopism 
sporządzonym zgodnie z przepisami wydanymi na podstawie art. 267 ust. 2 pkt 2 
Ustawy z dnia 20 lipca 2018 roku - Prawo o szkolnictwie wyższym i nauce (Dz. U. 
z 2018 poz. 1668 z późn. zm.), oznaczono numerami D1, D2, D3, D4. Symbolami 
SM2, SM3 i SM4 oznaczono materiał uzupełniający (ang. Supplementary Material) 
do oryginalnych artykułów naukowych. Dodatkowo umieszczono informacje 
naukometryczne: wskaźnik wpływu (ang. Impact Factor, IF) publikacji zgodnie 
z rokiem opublikowania i па przestrzeni ostatnich pięciu lat 
(IF5-letni) oraz punktację wg Ministerstwa Edukacji i Nauki (MEiN). Prowadzenie 
korespondencji z Redakcją Czasopisma oznaczono symbolem „*” 

Dysertacja ma formę zbioru spójnie powiązanych tematycznie następujących 


oryginalnych artykułów naukowych: 


[D1] K. Szymczyk, M. Szaniawska, K. Terpiłowski*, Determination of acoustical 
parameters of aqueous solution of Kolliphors binary mixtures using density, speed 
of sound, viscosity and surface tension measurements, Journal of Surfactants and 
Detergents, 2019, 22(5), 1163-1174 

IF2019/ LF5-letni: 1,654 / 1,822 

МЕТМ:оэз: 40 pkt/40pkt 

Mój wkład w powstanie tej pracy polegał na tworzeniu koncepcji artykułu, zaplanowaniu badań, 
wyborze technik badawczych, wykonaniu wszystkich badań eksperymentalnych i ich wizualizacji, 
analizie i interpretacji otrzymanych danych i ich przygotowaniu do publikacji, wykonaniu części 


obliczeń, przygotowaniu wstępnej wersji publikacji. 


[D2] K. Szymczyk*, M. Szaniawska, J. Krawczyk, Temperature effect on the 
adsorption and volumetric properties of aqueous solutions of Kolliphor ELP®, 


Molecules, 2020, 25(3), 743 


IF2020/ IF5-tetni? 4,411/5,110 

MEiN»oa:: 140 pkt/140 pkt 

Mój wkład w powstanie tej pracy polegał na tworzeniu koncepcji artykułu, zaplanowaniu badań, 
wyborze technik badawczych, wykonaniu wszystkich badań eksperymentalnych i ich wizualizacji, 
analizie i interpretacji otrzymanych danych i ich przygotowaniu do publikacji, wykonaniu części 


obliczeń, przygotowaniu wstępnej wersji publikacji. 


[D3] M. Szaniawska, K. Szymczyk, A. Zdziennicka, B. Jańczuk*, Adsorption 
properties and composition of binary Kolliphor mixtures at the air-water interface 
at different temperatures, Molecules, 2022, 27(3), 877 

IF2022/ IFS5 letni: 4,927 /5,110 

MEiN»o: 140 pkt/140 pkt 

Mój wkład w powstanie tej pracy polegał na tworzeniu koncepcji artykułu, zaplanowaniu badań, 
wyborze technik badawczych, wykonaniu wszystkich badań eksperymentalnych i ich wizualizacji, 
analizie i interpretacji otrzymanych danych i ich przygotowaniu do publikacji, wykonaniu części 


obliczeń, przygotowaniu wstępnej wersji publikacji. 


[D4] M. Szaniawska, К. Szymczyk*, A. Zdziennicka, В. Jańczuk, Thermodynamic 
parameters of berberine with Kolliphor mixtures adsorption and micellization, 
Molecules, 2023, 28(7), 3115 
IF 2023/ IF5-letni: 4,927 /5,110 
MEiN»oa:: 140 pkt/140 pkt 
Mój wkład w powstanie tej pracy polegał na tworzeniu koncepcji artykułu, zaplanowaniu badań, 
wyborze technik badawczych, wykonaniu wszystkich badań eksperymentalnych i ich wizualizacji, 
analizie i interpretacji otrzymanych danych i ich przygotowaniu do publikacji, wykonaniu części 


obliczeń, przygotowaniu wstępnej wersji publikacji. 


Sumaryczna liczba punktów przypisanych czasopismom przez 
МЕНМ:соэз: 460/ 460pkt 
Sumaryczna wartość wskaźnika IF/IF5-ietni prac w cyklu: 15,919/17,152 


4. WYKAZ STOSOWANYCH SYMBOLI I SKRÓTÓW 


ELP 


RH40 


Ber 


PTFE 


PMMA 


CMC 


Kolliphor ELP 

Kolliphor RH40 

berberyna 
politetrafluoroetylen 
poli(metakrylan metylu) 
krytyczne stezenie micelizacji 
napiecie powierzchniowe 


składowa napięcia powierzchniowego wynikająca z oddziaływań 


Liftshitza-van der Waalsa 


składowa napięcia powierzchniowego wynikająca z oddziaływań 


kwasowo- zasadowych w ujęciu Lewisa 


parametr elektronowo - akceptorowy składowej kwasowo - zasadowej 


napięcia powierzchniowego 


parametr elektronowo - donorowy składowej kwasowo - zasadowej 


napięcia powierzchniowego 

maksymalne nadmiarowe stężenie powierzchniowe Gibbsa 
standardowa swobodna energia Gibbsa procesu adsorpcji 
parametr oddziaływań międzycząsteczkowych w mieszanej miceli 


różnica między napięciem powierzchniowym rozpuszczalnika 


i substancji rozpuszczonej w CMC 

standardowa swobodna energia Gibbsa procesu micelizacji 
energia swobodna na granicy faz woda - powietrze 
nadmiarowa energia Gibbsa mieszanej miceli 

energia wymiany cząsteczek w mieszanej miceli 


ciśnienie wewnętrzne 


Ay, 45, ty, t5, y? 


stałe w funkcji wykładniczej stosowanej do opisu izoterm napięcia 


powierzchniowego 

odległość międzycząsteczkowa 

standardowa entropia procesu adsorpcji 
standardowa entalpia procesu adsorpcji 

napięcie międzyfazowe woda - ogon surfaktantu 
napięcie powierzchniowe woda - Ber 

napięcie miedzyfazowe Ber - ogon surfaktantu 
praca adhezji w układzie Ber - woda - ogon surfaktantu 
ułamek molowy 

napięcie powierzchniowe rozpuszczalnika 
impedancja akustyczna 


ciśnienie powierzchniowe surfaktantów 
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5. WPROWADZENIE I CEL ROZPRAWY DOKTORSKIEJ 


Berberyna (Ber, CooH18NO4*) ze względu na budowę chemiczną zaliczana jest do 
grupy alkaloidów, pochodnych izochinoliny. Jest to żółty bezzapachowy proszek 


o gorzkim smaku [1-4]. 


Berberys zwyczajny 


Cynowód chiński 


Gorzknik kanadyjski 


Rys. 1 Główne źródła berberyny (opracowanie własne) 


Nazwa berberyny wywodzi się od rośliny, z której najczęściej jest otrzymywana. 
Substancja ta występuje w kłączach, korzeniach i łodygach roślin stosowanych 
w tradycyjnej medycynie chińskiej, takich jak berberys zwyczajny (Berberis vulgaris), 
berberys indyjski (Berberis aristata) czy cynowód chiński (Coptis chinensis) (Rys. 1) 
[5-8]. Zawartość berberyny w tych roślinach wynosi od 0,05 mg/g do 96,10 mg/g 
[9,10]. Berberyna wykazuje bardzo szeroką aktywność biologiczną (Rys. 2), 
a pierwsze wzmianki na temat medycznego wykorzystania berberyny zawartej 
w Rhizoma coptidis sięgają roku 200 [11,12]. W szerokim spektrum jej działania 
można wyróżnić działanie  przeciwcukrzycowe, |przeciwnowotworowe, 
przeciwbiegunkowe, antyoksydacyjne, przeciwzapalne, przeciwwirusowe (m.in. 
w leczeniu SARS-CoV-2), a także wspomagające układ sercowo - naczyniowy oraz 


obniżające poziom lipidów [13-24]. Zaletą stosowania berberyny jest także fakt, że 


11 


powoduje znikome efekty uboczne i ma niską toksyczność, dzieki czemu może być 


stosowana jako potencjalny lek [25]. 


Działanie 
przeciwcukrzycowe 
. А 4 
Działanie | Działanie 


przeciwbiegunkowe 
ГА 


x | przeciwnowotworowe 


= Działanie 
przeciwzapalne 


Działanie \ Działanie 
(2) antyoksydacyjne ` przeciwdrobnoustrojowe 


Działanie na układ 
а + — 
sercowo-naczyniowy 


Dzialanie d 
rzeciwmiazdzycowe ?ф e. Өө 
P ycowe $8 e o 


Rys. 2 Aktywność biologiczna berberyny [13-24] (opracowanie własne) 


Pomimo szeregu pozytywnych właściwości, stosowanie berberyny jest 
ograniczone, głównie ze względu na jej niską rozpuszczalność w wodzie (5,27+0,29 
i 8,50+0,40 mmol/dm3, odpowiednio w temperaturze 298 i 310 K), jak również 
w buforze fosforanowym o pH=7 (4,05+0,09 i 9,69+0,37 mmol/ dm? w tych samych 
temperaturach) [26-28]. Wynika to głównie z obecności w jej strukturze 
czwartorzędowego kationu amoniowego, a także ze sztywnej i płaskiej struktury 
berberyny [29]. Konsekwencją słabej rozpuszczalności berberyny w wodzie jest jej 
niska biodostępność, np. w pracy Chena i wsp. wynosi ona 0,68 [30], według Liu 
i wsp. 0,36 [31], podobnie w badaniach Fenga i wsp. - 0,37 [32]. Według danych 
literaturowych biodostępność berberyny jest mniejsza niż 1% [33-34]. Tak niska 
wartość tego parametru związana jest głównie z bardzo intensywnym 
metabolizmem, któremu podlega berberyna w jelicie cienkim i wątrobie [35], jak 
również wartość pKa=~15, co oznacza, że w warunkach fizjologicznych występuje 
głównie w formie zjonizowanej [36]. Do wad tego związku zaliczyć można także 
czułość na światło i wysoką temperaturę, które prowadzą do jej degradacji [37]. 


W systemie klasyfikacji leków berberyna zaliczana jest do III klasy związków (Rys. 
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3), charakteryzujacych się bardzo słabą przenikalnością [38]. System ten klasyfikuje 
substancje lecznicze pod względem ich właściwości biofarmaceutycznych, tj. 
przenikalności przez bariery biologiczne i rozpuszczalności. Substancje klasyfikuje 
się sprawdzając, czy dawka przewidziana rozpuszcza się w 250 cm? wody 
(w zakresie pH od 1 do 7,5), a także określając przenikalność tej substancji przez 


bariery biologiczne (np. warstwę komórek w hodowli in vitro) [39]. 


Wysoka przenikalność 


Grupa II 


Czynnik limitujący: 
niska rozpuszczalność 


Grupa I 
Czynnik limitujący: 
brak 


System klasyfikacji 


Wysoka rozpuszczalność J B 
biofarmaceutycznej 


Niska rozpuszczalność 


Grupa IV 
Czynnik limitujący: 
niska rozpuszczalność, 
niska przenikalność 


Grupa III 


Czynnik limitujący: 
niska przenikalność 


Niska przenikalność 


Rys. 3 System klasyfikacji biofarmaceutycznej (opracowanie własne) 


Głównym wyzwaniem w zastosowaniu berberyny jako substancji leczniczej jest 
więc poprawa jej rozpuszczalności, a co za tym idzie biodostępności za pomocą 
różnych metod. W literaturze można znaleźć informacje o wykorzystywaniu do 
tego celu wielu technik, m.in. tworzenia nowych pochodnych berberyny [40], 
mikro- i nanoemulsji [41,42], stałych lipidowych nanocząstek [43] czy liposomów 
[44]. Należy mieć na uwadze, że metody te poza zaletami mają także wady 
i ograniczenia. W metodzie współrozpuszczania stosuje się organiczne, toksyczne 
rozpuszczalniki, a substancja rozpuszczana może mieć mniejszą stabilność, co 
w późniejszym etapie po procesie może niekorzystnie wpływać na wchłanianie 
[45]. Z kolei tworzenie mikroemulsji jest procesem skomplikowanym i trudnym do 
walidacji. Stosowanie roztworów olejowych może spowodować, że dany preparat 
nie może być podany dożylnie. W przypadku mikroemulsji efekt terapeutyczny 


często jest uzależniony od oddziaływań substancji leczniczej ze związkiem 
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chemicznym tworzacym micele czy liposomy [46,47]. Tworzenie soli takze nie jest 
metodą odpowiednią dla każdego rodzaju związków [48]. Jedną 
z wad tworzenia soli jest to, że zastosowanie przeciwjonu o dużej masie 
cząsteczkowej z lekiem o niskiej aktywności może spowodować, że wymagana jest 
wyższa dawka. Kolejnym problemem jest tendencja do tworzenia hydratów 
i odmian polimorficznych. Istnieją również przypadki, w których sole mocnych 
kwasów lub zasad są dobrze rozpuszczalne w wodzie, będąc jednocześnie bardzo 
higroskopijne, co prowadzi do niestabilności formulacji [49]. Poszukiwanie 


prostych i skutecznych metod jest zatem pożądane i w pełni uzasadnione. 


Stosunkowo nowym rozwiązaniem jest wykorzystanie związków 
powierzchniowo czynnych (surfaktantów) do ekstrakcji aktywnych składników 
pochodzenia roślinnego o wysokim stężeniu [50]. Zastosowanie surfaktantów do 
poprawy rozpuszczalności w wodzie niektórych substancji jest obecnie 
przedmiotem badań wielu naukowców np. Krstonośić i wsp. zastosowali 
Poloxamer 407 i Brij 20 do ekstrakcji kwasu chlorogenowego [51], Giacomo i wsp. 
użyli Tritonu X-114 oraz Genapolu X-080 do ekstrakcji micelarnej papainy [52], 
Shestopalova z zespołem zastosowali Triton X-100 w celu poprawy 
rozpuszczalności azorubiny [53], a Wu i wsp. wykorzystali Span 20 do zwiększenia 
rozpuszczalności rytonawiru [54]. Rozwiązanie to wydaje się jak najbardziej 
uzasadnione, gdyż od dawna przemysł spożywczy, kosmetyczny 
i farmaceutyczny wykorzystuje związki powierzchniowo czynne (surfaktanty) do 
tworzenia różnych układów np. emulsji lub pian [55]. Surfaktanty są substancjami 
bardzo często stosowanymi w przemyśle kosmetycznym i farmaceutycznym, ze 
względu na swoje właściwości [56-58], zdolność gromadzenia się na granicy faz 
przy bardzo niskich stężeniach oraz wpływ na właściwości powierzchniowe cieczy, 
w której są rozpuszczone [59]. Są to substancje o budowie asymetrycznej, składają 
się z części hydrofilowej (jonowej lub niejonowej) oraz hydrofobowej. Po 
przekroczeniu granicznego stężenia nazywanego krytycznym stężeniem micelizacji 
(CMC), cząsteczki surfaktantów łączą się w agregaty nazywane micelami. Micele 
są wykorzystywane w zastosowaniach biologicznych ze względu na ich unikalne 


właściwości, takie jak biokompatybilność i zdolność do ochrony leków oraz 
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ukierunkowanego ich dostarczania [60]. Ogromną zaletą miceli jest umiejętność do 
solubilizowania w ich rdzeniu słabo rozpuszczalnych w wodzie lub 


hydrofobowych leków, co zwiększa ich biodostępność [61]. 


Surfaktanty mogą być stosowane pojedynczo lub w mieszaninach o różnym 
ułamku molowym. Od charakteru chemicznego składników zależą oddziaływania 
poszczególnych składników, a więc synergetyczne i antagonistyczne działanie. 
Kiedy oddziaływania międzycząsteczkowe są mocniejsze w mieszaninie niż 
w roztworach pojedynczych składników, zjawisko takie nazywa się 
synergetyzmem [59,62]. Powszechnie w przemyśle farmaceutycznym 
i kosmetycznym stosuje się mieszaniny surfaktantów o działaniu synergetycznym, 
ponieważ są one skuteczne w znacznie niższych stężeniach niż pojedyncze 
składniki [63]. Efekt taki często jest rozważany w oparciu o teorie roztworów 
nieidealnych, np. przybliżoną teorię roztworów rzeczywistych (ang. regular solution 
approximation), szczególnie za pomocą parametru oddziaływań 


miedzyczasteczkowych, В [64]. 


W formulacjach farmaceutycznych wykorzystuje sie surfaktanty o nazwie 
handlowej Kolliphor. Ich rola w takich układach związana jest głównie 
z hamowaniem aktywności P-glikoproteiny, dzięki czemu poprawia się 
rozpuszczalność substancji, a co za tym idzie wzrasta biodostępność i skuteczność 
[65]. Są to związki zatwierdzone przez amerykańską Agencję Żywności i Leków 
(FDA) do doustnego, domięśniowego i dożylnego stosowania w formulacjach 


farmaceutycznych leków [66]. 


Najpopularniejszym z tej grupy związków jest Kolliphor ELP (oczyszczona 
forma Kolliphoru EL, Cremophor ELP, etoksylowany kwas (R)-12-hydroksy-cis-9- 
oktadekenowy). Surfaktant ten otrzymywany jest w reakcji kwasu (R)-12-hydroksy- 
cis-9-oktadekenowego z tlenkiem etylenu w stosunku molowym 1:35 [94]. 
Głównym składnikiem tego związku jest rycynooleinian glikolu polietylenowego 
i glicerolu, który razem z estrami tłuszczowymi glikolu polietylenowego tworzy 


grupę hydrofobową tego surfaktantu (Rys. 4a). 
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Kolliphor ELP (ELP) jest stosowany jako środek wspomagający przenikalność 
i biodostępność wielu leków i substancji biologicznie aktywnych [95,96]. 
W literaturze znaleźć można wiele prac, w których wykorzystywano ten surfaktant 
do poprawy biodostępności i rozpuszczalności, m.in. wardenafilu [67], pochodnych 
artemizyny [68], cyklosporyny [69], diazepamu, propofolu [70], paklitakselu [71], 
symwastatyny [72], czy indirubiny [73]. Jest także składnikiem szczepionek mRNA 
przeciw COVID-19 firm Pfizer BioNTech i Moderna [74]. 


х+у+2=40 


Rys. 4 Struktury Kolliphoru ЕІР (а) і КН40 (b) 


Kolejny surfaktant z tej grupy, Kolliphor RH40, (RH40) (Rys. 4b) (uwodorniony 
kwas  (R)-12-hydroksy-cis-9-oktadekenowy), jest związkiem chemicznym 
otrzymywanym w wyniku reakcji 1 mola uwodornionego kwasu (R)-12-hydroksy- 
cis-9-oktadekenowego z 40 molami tlenku etylenu [94]. Kolliphor RH40 
wykorzystywany również jest jako solubilizator i środek poprawiający 
biodostępność substancji [75,76]. Kolliphor RH40 z powodzeniem jest stosowany 
jako składnik samoemulgujących systemów dostarczania leków (ang. Self- 


Emulsifying Drug Delivery Systems, SEDDS), m. in. w formulacjach zawierających 
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skwalen [77] czy sildenafil [78]. Kolliphor RH40 jest stosowany także jako środek 
poprawiający biodostępność kamptotecyny [79], itrakonazolu [80], karwedilolu 
[81], czy celekoksybu [82]. 


Mimo, że w literaturze występuje wiele danych dotyczących surfaktantów typu 
Kolliphor [83-92], to trudno jednak doszukać się informacji dotyczących ich 
współdziałania z berberyną. Istnieją badania dotyczące mieszanin berberyny 
np. z surfaktantami, Brij-S20 [93], SLS czy Tween 80 [26]. Brak jest jednak 
kompleksowej analizy tego typu układów w szerokim zakresie stężeń 
i temperatur. Dodatkowo właściwości fizykochemiczne samej berberyny, jak 
również ELP i RH40 nie zostały dotychczas jednoznacznie zbadane. Zatem 
z praktycznego i teoretycznego punktu widzenia poznanie właściwości 
adsorpcyjnych i objętościowych ELP, RH40 i Ber, a także ich mieszanin 
w środowisku wodnym jest jak najbardziej uzasadnione i stało się celem 
przedstawionej rozprawy doktorskiej [D1-D4]. W pierwszej kolejności za pomocą 
pomiarów napięcia powierzchniowego przeprowadzono badania wodnych 
roztworów Ber (1-10? mol/dn? - 1.104 mol/ ат) w temperaturach 293-313 K oraz 
ELP i RH40 w szerokim zakresie ich stężenia (1.10% - 1.10? mol/ dm?) i w różnych 
temperaturach (T = 293 - 318 K). Zmierzono także kąty zwilżania na powierzchni 
modelowych ciał stałych jakimi są apolarny politetrafluoroetylen (PTFE) oraz 
monopolarny polifmetakrylan metylu) (PMMA) (293 K). Uzyskane wyniki 
umożliwiły przeprowadzenie termodynamicznej analizy procesu adsorpcji 
składników mieszanin na granicy faz woda-powietrze, ale również ciało stałe-woda 
i ciało stałe-powietrze, ustalenie relacji pomiędzy ich adsorpcją na tych granicach 
faz oraz określenie składu mieszanych monowarstw adsorpcyjnych. Następnie, 
stosując dodatkowo pomiary gęstości, lepkości, prędkości dźwięku określono 
właściwości adsorpcyjne i objętościowe dwuskładnikowych mieszanin ELP i RH40 
o ułamku molowym w fazie objętościowej a = 0,2; 0,4; 0,6; 0,8 w zakresie stężeń 
1:106- 1.10? mol/dm$ i temperaturach 293-318 K celem wyboru do dalszych badań 
z Ber mieszaniny surfaktantów, dla której istnieje efekt synergetyczny w redukcji 
napięcia powierzchniowego i procesie tworzenia micel. Zmierzono także kąty 


zwilżania wody, formamidu i dijodometanu na sprasowanej Ber w celu 
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wyznaczenia składowych i parametrów napięcia powierzchniowego Ber. 
Końcowym etapem badań była ścisła charakterystyka fizykochemiczna mieszanin 
Ber z surfaktantami ELP, RH40, a także ELP + RH40 na podstawie pomiarów 
napięcia powierzchniowego, przewodności właściwej, gęstości, widm emisyjnych 
Ber w wodnych roztworach ELP i RH40 oraz ELP + RH40 w szerokim zakresie 
stężeń surfaktantów (1:10% - 1.102 mol/dm?) i temperatur (293-313 K). Ма 
podstawie otrzymanych wyników określono zachowanie sie berberyny 
w środowisku wodnym bez i w obecności surfaktantów. Stan równowagi 
pomiędzy roztworem monomerycznym a micelarnym surfaktantów, ich 
dwuskładnikowych mieszanin, a także mieszanin z Ber poddano analizie 


termodynamicznej w oparciu o otrzymane rezultaty. 
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6. OMOWIENIE WYNIKOW 
6.1. Właściwości fizykochemiczne ELP, RH40 i Ber 


Zachowanie się surfaktantów, ich mieszanin, a także mieszanin surfaktantów 
z innymi związkami organicznymi w środowisku wodnym zależy nie tylko od 
rodzaju grup chemicznych obecnych w cząsteczkach tych związków, ale również 
od ich liczby i polarności. Na właściwości adsorpcyjne i agregacyjne tego typu 
roztworów szczególny wpływ ma wielkość i powierzchnia kontaktu części 
hydrofilowej i hydrofobowej surfaktantów i dodatków z innymi cząsteczkami, 


a także składowe i parametry ich napięcia powierzchniowego. 


Według van Ossa i wsp. [97-100] napięcie powierzchniowe cieczy i ciał stałych 
można podzielić na dwie składowe. Pierwszą, która stanowi napięcie 
powierzchniowe każdej cieczy lub ciała stałego wynika z oddziaływań 
miedzyczasteczkowych Lifshitza - van der Waalsa (y’”) i drugą, która wynika 
z oddziaływań kwasowo - zasadowych (y^P). Składowa kwasowo-zasadowa 
zależy od parametrów eletronowo-akceptorowego ( y* ) i elektronowo - 
donorowego (y`). Ponieważ van Oss i wsp. sugerują, że udział oddziaływań dipol- 
dipol i dipol-dipol indukowany w składowej Lifshitza-van der Waalsa napięcia 
powierzchniowego jest mniejszy niż 2%, wynika z tego, iż składowa y^ napięcia 
powierzchniowego ciała stałego i cieczy zależy tylko od dyspersyjnych 


oddziaływań międzycząsteczkowych. 


Fowkes [101] stwierdził, że dyspersyjne oddziaływania międzycząsteczkowe na 
granicy faz można wywnioskować na podstawie oddziaływań jednego 
elementu / fragmentu z dwunastoma jego sąsiadami oraz odległości pomiędzy 
nimi. W przypadku wodnych roztworów ELP i RH40 na granicy faz roztwór- 
powietrze można wyróżnić cząsteczki wody oraz grupy —CH3, —CH; —, = СН —, 
= CO i—OH. Ponieważ oddziaływania Liftshitza - van der Waalsa w pierwszym 
przybliżeniu są równe oddziaływaniom dyspersyjnym, to zależą one od rodzaju 
grup funkcyjnych znajdujących się w cząsteczce surfaktantu i napięcie 


powierzchniowe surfaktantu jest związane z orientacją jego cząsteczek na danej 
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granicy faz. Biorąc pod uwagę ten fakt van Oss i Constanzo [98] stwierdzili, ze 
napięcie powierzchniowe surfaktantu jest różne w zależności od orientacji jego 
cząsteczek na granicy faz surfaktant-powietrze. Stąd też napięcie powierzchniowe 
związku powierzchniowo czynnego, którego cząsteczki są zorientowane grupą 
hydrofobową (ogonem) w kierunku fazy gazowej można nazwać napięciem 
powierzchniowym ogona surfaktantu (ут, T - tail). Dla surfaktantu, którego 
cząsteczki są zorientowane w stronę fazy gazowej grupami hydrofilowymi (głową), 
jego napięcie powierzchniowe to napięcie powierzchniowe głowy surfaktantu 


(Yu, H - head). 


Ze wzgledu na bardzo silne oddziatywania hydrofobowe pomiedzy ogonami 
czasteczek ELP i RHAO ich zmierzone napiecie powierzchniowe mozna traktowaé 
jako napiecie powierzchniowe glowy [D4]. W zwiazku z tym prawdopodobienstwo 
orientacji czasteczek ELP i RH40 na granicy faz surfaktant-powietrze glowami 
skierowanymi w strone fazy gazowej jest najbardziej realne. Biorac to po uwage 
wyznaczono skladowe i parametry napiecia powierzchniowego glowy ELP i RH40 
w temperaturze 293K poprzez pomiary ich kata zwilzania na powierzchni 
modelowych ciał stałych, politetrafluoroetylenu (PTFE) i poli(metakrylanu) metylu 
(PMMA) [D4]. Okazało się, że wyznaczone w ten sposób składowe napięcia 
powierzchniowego ELP i RH40 (Tab. 1, [D4]) są zbliżone do składowych napięcia 
powierzchniowego głowy Tritonu X-165 (TX165) [102], co potwierdza założenie, że 
cząsteczki ELP i RH40 na granicy faz woda-powietrze są skierowane głowami 


w stronę fazy gazowej. 


W cząsteczce berberyny (Ber), w przeciwieństwie do ELP i RH40, trudno jest 
wyróżnić część hydrofobową i hydrofilową, dlatego też należy oczekiwać, że 
cząsteczki Ber na granicy faz woda-powietrze będą ułożone równolegle. Napięcie 
powierzchniowe berberyny (Ber) określono na podstawie kąta zwilżania 
zmierzonego dla wody (56,2), formamidu (38,4) i dijodometanu (46,1?) па 
sprasowanej berberynie metodą van Ossa i wsp. [97-99]. Znając składową Lifshitza 


LW 


- van der Waalsa (y*' ) napięcia powierzchniowego wody, formamidu 


i dijodometanu oraz parametry elektronowo - akceptorowe (y*) i elektronowo - 


20 


donorowe (y ) tego napięcia [103], otrzymano napięcie powierzchniowe Ber 
(46,52 mN/m), składową y*” (36,42 mN/m) oraz parametry ү? (1,56 mN/m) 
iy (16,38 mN/m) tego napięcia [D4]. Z przeprowadzonych obliczeń wynika, że 
napięcie powierzchniowe Ber niewiele różni się od tego dla np. TX165 [104], jeśli 
cząsteczki TX165 skierowane są grupą hydrofilową w stronę fazy gazowej. Należy 
jednak podkreślić, że udział składowej y"" w napięciu powierzchniowym Ber jest 
mniejszy, a y^? większy od tego dla TX165. Dodatkowo, udział składnika y^ 
w napięciu powierzchniowym Ber wynosi 78% ito powoduje słabą rozpuszczalność 
berberyny w wodzie [D4]. Należy podkreślić, że berberyna w niewielkim stopniu 
obniża napięcie powierzchniowe wody (Rys. 1, [D4]), a maksymalne obniżenie 
odpowiada stężeniu Ber w wodzie bliskiemu 1-104 mol/dm3. Przy większym 
stężeniu Ber w roztworze wodnym napięcie powierzchniowe roztworu jest stałe, 
jednak punktu przegięcia izotermy napięcia powierzchniowego wodnego roztworu 


berberyny nie można traktować jako CMC [D4]. 


Rozpuszczalność związków organicznych w wodzie zależy nie tylko od ich 


LW oraz y^P, ale także od wielkości 


napięcia powierzchniowego i wielkości y 
i powierzchni kontaktu grup hydrofilowych i hydrofobowych ich cząsteczek. 
Wielkość cząsteczek ELP, RH40 i Ber obliczono na podstawie długości wiązań, kąta 
pomiędzy wiązaniami i odległości między cząsteczkami [D4]. Wykazano, 
że objętość jednej cząsteczki ELP, RH40 i Ber można określić na podstawie objętości 
sześcianów, w które można wpisać poszczególne fragmenty danej cząsteczki. 
Objętość otrzymanych w ten sposób cząsteczek ELP, RH40 i Ber wyniosła 
odpowiednio 4378,64 i 4825,21 i 428,25 A3. Biorąc pod uwagę otrzymane w ten 
sposób wielkości badanych cząsteczek wyznaczono ich powierzchnie kontaktu 
z cząsteczkami wody. Powierzchnia ta wyniosła odpowiednio 3779,9 A2 dla ELP, 
4210,2 A2 dla RH40 i 277,45 A2 dla Ber. Pamiętając, że powierzchnia kontaktu 
cząsteczek wody jest równa 10 A2 [105], teoretycznie wykazano, że 378 cząsteczek 
wody może kontaktować się z ELP, 421 z RH40 a 28 z Ber [D4]. Ponieważ 
powierzchnia kontaktu ogona cząsteczek ELP i RH40 wynosi 967,7 À? [D3], 
jednocześnie może kontaktować się z nim 97 cząsteczek wody. Z drugiej jednak 


strony z danych literaturowych wynika, że grupa oksyetylenowa w cząsteczce 
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związku powierzchniowo czynnego moze być połączona silnymi wiązaniami 
wodorowymi z dwiema cząsteczkami wody oraz słabymi wiązaniami z trzema 
cząsteczkami [D4, 106, 107]. Biorąc to pod uwagę wykazano, że w zależności od 
konfiguracji łańcuchy oksyetylenowe w cząsteczce ELP mogą być otoczone 
najwyżej około 245 cząsteczkami wody, a łańcuchy RH40 około 280 cząsteczkami, 
przy czym silnymi wiązaniami wodorowymi z atomami tlenu, 
z wyjątkiem tego w grupie —0H, z ELP połączone jest 82 cząsteczki wody, 
a 92 z RH40 [D3]. Obliczenia te są zgodne z wartościami liczb hydratacji ELP i RH40 
wyznaczonymi na podstawie pomiarów prędkości dźwięku i gęstości [D1]. 
Dodatkowo znaczącą hydratację grup hydrofilowych ELP przy jego stężeniu 
w roztworze mniejszym niż CMC i w temperaturze 293K wykazano na podstawie 
pomiarów lepkości i zarejestrowanych widm emisyjnych pirenu [D2]. 
Udowodniono także, że hydratacja ta w tym zakresie stężeń istotnie wpływa na 


konformację i zwartą strukturę łańcuchów ELP. 


Liczba cząsteczek wody otaczających grupy  hydrofilowe związków 
powierzchniowo czynnych sprawia, że są one w dużej mierze rozpuszczalne 
w wodzie. Z drugiej jednak strony liczba cząsteczek wody otaczających grupy 
hydrofobowe cząsteczek ELP i RH40 jest siłą napędową ich adsorpcji na granicy faz 
woda-powietrze i procesu micelizacji. Siła ta zależy od składowych 
i parametrów napięcia powierzchniowego wody i surfaktantów, a także od 
wielkości powierzchni kontaktu. Wyznaczona powierzchnia kontaktu cząsteczek 
ELP i RH40 przez fazę wodną wynosi odpowiednio 951,62 i 1054,58 A2 [D4]. 
Powierzchnia kontaktu pomiedzy dwiema czasteczkami Ber wynosi 108,96 A? ijest 
mniejsza od maksymalnej powierzchni kontaktu pomiędzy ogonami dwóch 


cząsteczek ELP i RH40, która równa jest 324,71 А2 [D4]. 


Duży rozmiar głowy cząsteczek ELP i RH40 gwarantuje dobrą rozpuszczalność 
tych związków w wodzie, natomiast ogona - powoduje dobrą aktywność 
powierzchniową i skłonność do tworzenia miceli w środowisku wodnym [D2, D4]. 
Przykładowo maksymalne nadmiarowe stężenie powierzchniowe Gibbsa ELP 


i RH40 w temperaturze 293 K jest większe od tego dla np. Tritonów, zaś proces 
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tworzenia micel ELP i RH40 rozpoczyna sie przy stezeniach mniejszych niz dla 
Tritonów [D2, D3, 108]. Warto podkreślić, ze nadmiarowe stężenie 
powierzchniowe Gibbsa berberyny na granicy faz woda-powietrze wynosi 
odpowiednio 1,5-10%, 1,48.10% i 1,46-10-% mol/ dm? w temperaturze 293, 303 i 313 K, 
a obliczona na podstawie tych wartości minimalna powierzchnia zajmowana przez 
cząsteczkę Ber na tej granicy faz jest zbliżona do powierzchni kontaktu berberyny 


w przypadku jej równoległej orientacji [D4]. 


6.2. Właściwości adsorpcyjne i objętościowe dwuskładnikowych 


mieszanin ELP+RH40 


Przeprowadzone pomiary napięcia powierzchniowego y,ywodnych roztworów 
dwuskładnikowych mieszanin ELP i RH40 o ułamku molowym ELP w fazie 
objętościowej a = 0,2; 0,4; 0,6; 0,8 wykazały, że przy danym stężeniu mieszaniny 
napięcie to w większości przypadków maleje niemal liniowo w funkcji temperatury 
[Rys. 1, D3]. Świadczy to o tym, że w badanym zakresie temperatur (293-318 K) 
dehydratacja RH40 i ELP w mieszaninie jest nieznaczna. Z drugiej jednak strony 
brak jest liniowej zależności pomiędzy napięciem powierzchniowym wodnego 
roztworu mieszaniny RH40 i ELP przy danym stężeniu w funkcji jej składu (Rys. 4, 
[D3]) - obserwuje się dodatnie i ujemne odchylenia wartości y;y od zależności 
liniowej. Świadczy to o możliwości wystąpienia zmian liczby hydratacji grupy 
hydrofilowej surfaktantów i/lub zmian w upakowaniu cząsteczek surfaktantów 
w wyniku zmiany ich konfiguracji. Z danych literaturowych jasno wynika, że 
możliwość tworzenia się różnych konfiguracji surfaktantów о dużych 
i rozgałęzionych cząsteczkach, typu ELP i RH40, jest większa niż w przypadku 
klasycznych surfaktantów liniowych [109, 110]. Zjawisko to jest szczególnie 
widoczne w przypadku surfaktantów zawierających w swoich cząsteczkach grupy 


oksyetylenowe [111, 112]. 
Z przeprowadzonych pomiarów wynika, że izotermy napięcia 
powierzchniowego wodnych roztworów mieszanin ELP i RH40 o danym składzie 


w różnych temperaturach są prawie równoległe (Rys. 1b-e, [D3]). Wykazano, że 


wszystkie te izotermy można opisać funkcją wykładniczą drugiego rzędu (Rys. S1- 
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56, [SM3]), a stałe w równaniu tej funkcji powiązać ze składowymi i parametrami 
napięcia powierzchniowego głowy i ogona ELP oraz RH40. Wniosek ten dotyczy 
zwłaszcza stałej Уур, której wartości dla wszystkich mieszanin i badanych 
temperatur są zbliżone do minimalnych wartości y,y danej mieszaniny. Ponieważ 
minimalne wartości y;, zależą od napięcia powierzchniowego ogona surfaktantu, 
można stwierdzić, że yo jest również związana z napięciem powierzchniowym 


grupy hydrofobowej surfaktantu [D3]. 


Wyznaczone eksperymentalnie izotermy napięcia powierzchniowego mieszanin 
ELP i RH40 opisano także za pomocą równania Szyszkowskiego i przeprowadzono 
próbę ich obliczenia teoretycznego na podstawie wartości napięcia 
powierzchniowego ELP i RH40, przy założeniu, że badane mieszaniny zachowują 
się jak roztwory idealne. Wykazano, że przy takim samym stężeniu ELP i RH40 dla 
wszystkich badanych mieszanin obliczone przy takim założeniu wartości y;y są 
większe od zmierzonych (Rys. 52-55, [5М3]). Różnica jednak pomiędzy 
zmierzonymi a obliczonymi wartościami Ууц, zależy od składu mieszaniny 
i np. w przypadku mieszaniny ELP i RH40 o ułamku molowym ELP w fazie 
objętościowej równym 0,6, obliczone wartości y,y są prawie identyczne 
z wartościami zmierzonymi (Rys. 54, [SM3]). Dodatkowo dla mieszaniny o danym 
składzie różnica pomiędzy obliczoną a zmierzoną wartością y;, praktycznie nie 
zmienia się w funkcji temperatury. Stosując równanie Fainermana i Millera [113- 
114] również teoretycznie wyznaczono izotermy napięcia powierzchniowego 
mieszanin ELP i RH40 [D3]. Równanie to oparte na izotermach napięcia 
powierzchniowego roztworów pojedynczych składników mieszanin pozwala na 
przewidywanie izoterm napięcia powierzchniowego prostych mieszanin 
dwuskładnikowych o stałym składzie a różnym stężeniu, dla których powierzchnie 
zajmowane przez jeden mol indywidualnych składników są do siebie zbliżone. 
Wykazano, że równanie Fainermana i Millera z powodzeniem można zastosować 
dla badanych mieszanin w zakresie ich stężenia w fazie objętościowej stężeń od 


0 do CMC. 


24 


Różnica pomiędzy wyznaczonymi eksperymentalnie a obliczonymi teoretycznie 
wartościami y,y Świadczy о tym, że roztwory mieszaniny ELP i RH40 nie 
zachowują się jak roztwory idealne. Potwierdzeniem tego wniosku są np. wartości 
standardowej swobodnej energii Gibbsa procesu adsorpcji, AG24, (Tab. 52, [D3]) 
i maksymalnego nadmiarowego stężenia powierzchniowego Gibbsa ( Г" ) 
zastosowane do obliczeń үг, za pomocą równania Szyszkowskiego [D3]. Wartości 
r™* nie zmieniają się liniowo w funkcji składu mieszaniny surfaktantów (Rys. 513, 
[SM3]) i mieszczą się w przedziale od 2,44-10 do 3,12:10:9 mol/m?. Odpowiadają 
także minimalnej powierzchni zajmowanej przez jedną cząsteczkę na granicy faz, 
która mieści się w przedziale od 53,2 do 68 А2, Świadczy to zatem o istnieniu 
zależnego od składu mieszaniny efektu synergetycznego w redukcji napięcia 
powierzchniowego wody. Dla badanych mieszanin efekt ten został oszacowany za 
pomocą wartości parametru oddziaływań miedzyczasteczkowych В“ obliczonego 
na podstawie modelu Hua i Rosena [64]. W badanym zakresie stężeń i temperatur 
wyznaczony parametr ВМ przyjmuje zarówno wartości ujemne (efekt 


synergetyczny), jak i dodatnie (efekt antagonistyczny) (Tab. 51, [SM3)). 


W przypadku mieszaniny ELP i RH40 o ułamku molowym ELP w fazie 
objętościowej równym 0,6 wszystkie wartości B" są dodatnie i istnieje niemal 
liniowa zależność pomiędzy 8“ a napięciem powierzchniowym roztworów 
i temperaturą (Rys. S16d, [SM3]). Wskazuje to na brak efektu synergetycznego tej 
mieszaniny w redukcji napięcia powierzchniowego wody. Dla pozostałych 
badanych mieszanin RH40 i ELP efekt ten występował, ale nie w każdej 
zastosowanej w badaniach temperaturze. Interesującym jest, że w przypadku 
mieszaniny o ułamku molowym ELP w fazie objętościowej równym 0,8 wraz ze 
wzrostem stężenia i spadkiem napięcia powierzchniowego wartości 8“ maleją 
(Tab. S1, [SM3]). Ujemne w całym zakresie temperatur wartości parametru 
В zaobserwowano dla tej mieszaniny przy y;y = 45 mN/m, a więc w pobliżu CMC 
[D3]. Wykazano, że w przypadku badanych mieszanin ELP i RH40 efekt 
synergetyczny w redukcji napięcia powierzchniowego wody wynika 
prawdopodobnie ze zmian konfiguracji cząsteczek ELP i RH40 w mieszanej 


monowarstwie adsorpcyjnej w porównaniu z monowarstwą utworzoną przez 
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pojedyncze surfaktanty i/lub zmiany liczby wiązań wodorowych na skutek 


dehydratacji grup hydrofilowych i hydrofobowych surfaktantów. 


Istnienie tego efektu jest ściśle związane ze zmianami składu i stężenia mieszanej 
monowarstwy adsorpcyjnej na granicy faz woda-powietrze w funkcji stężenia 
mieszaniny i temperatury przy stałym składzie w fazie objętościowej. Skład 
mieszanej monowarstwy adsorpcyjnej utworzonej przez mieszaniny ELP i RH40 
obliczono na podstawie modelu Hua i Rosena (Tab. 51, [SM3]) oraz wartości 
ciśnienia powierzchniowego surfaktantów (7). Wykazano, że dla danego stężenia 
i temperatury wartości te są zbliżone, a w zakresie stężeń mieszaniny w fazie 
objętościowej odpowiadającym nienasyconej monowarstwie adsorpcyjnej, tj. od 
zera do 1.109 mol/dm? występuje niezależna adsorpcja ELP i RH40 i ułamek 
molowy ELP w monowarstwie jest mniejszy od tego w fazie objętościowej (Tab. S1, 
[SM3]). Przy stężeniach odpowiadających nasyconej monowarstwie adsorpcyjnej 


ułamek molowy ELP jest zbliżony lub większy od a [D3]. 


Wyznaczone na podstawie izoterm napięcia powierzchniowego w temperaturze 
293K wartości krytycznego stężenia micelizacji СМС, Псмс (różnicy między 
napięciem powierzchniowym rozpuszczalnika yo i napięciem powierzchniowym 
substancji rozpuszczonej y w punkcie CMC), AG, (standardowa swobodna 
energia Gibbsa procesu micelizacji), AGs n (energia swobodna na granicy faz woda 
- powietrze), ВМ (parametr oddziaływań międzycząsteczkowych w mieszanej 
miceli), g” (nadmiarowa energia Gibbsa mieszanej miceli) i £&" (energia wymiany 
cząsteczek w mieszanej miceli) potwierdziły także istnienie efektu synergetycznego 
w procesie tworzenia mieszanych micel mieszaniny ELP i RH40 o ułamku 
molowym ELP w fazie objętościowej równym 0,8 (Tab. 1, [D1]). Również wzrost 
obliczonych na podstawie pomiarów gęstości, lepkości i prędkości dźwięku 
roztworów tej mieszaniny wartości impedancji akustycznej (2) i ciśnienia 
wewnętrznego (л;) powyżej stężenia 10? mol/dm? (Rys. 1-Rys. 4, [D1]) wskazały 
na silne oddziaływania pomiędzy składnikami badanych binarnych roztworów 


micelarnych. Warto podkreślić, że te oddziaływania międzycząsteczkowe, jak 


również gęstsze upakowanie surfaktantów w mieszanych micelach ELP i RH40 
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przy ułamku molowym ELP w fazie objętościowej równym 0,8 powyżej 
C = 10? mol/dm3, stwierdzono również na podstawie wartości odległości 
międzycząsteczkowej (molecular free lenght, Le) obliczonych z modelu Chauhana 
i Sharmy [107]. Wartości Lę, czyli odległości pomiędzy powierzchniami dwóch 
cząsteczek, zależą zarówno od oddziaływań międzycząsteczkowych, jak 
i wewnątrzcząsteczkowych pomiędzy składnikami roztworu i maleją znacznie 
powyżej CMC mieszaniny, potwierdzając istnienie efektu synergetycznego 
w procesie tworzenia mieszanych micel. Jest to bardzo istotny wniosek, który 
połączony z rozważaniami na temat zmiany liczby hydratacji czy stopnia wiązania 
wody z micelami ELP, RH40 i ich mieszaniny [D1-D3] ma duże znaczenie 
w wyjaśnieniu oddziaływań surfaktant - berberyna i może decydować 


o potencjalnym praktycznym zastosowaniu tego typu układów. 


6.3. Właściwości adsorpcyjne i objętościowe ELP + Ber, RH40 + Ber 
i ELP + RH40 + Ber 


Biorąc pod uwagę udział oddziaływań międzycząsteczkowych Lifshitza - van 
der Waalsa w napięciu powierzchniowym wody, ELP, RH40 i Ber, można 
stwierdzić, że obniżenie napięcia powierzchniowego wody na skutek adsorpcji 
badanych związków na granicy faz woda-powietrze zależy jedynie od redukcji 
składowej kwasowo-zasadowej napięcia powierzchniowego wody [D4]. Składowa 
y" napięcia powierzchniowego ELP (27,51 mN/m), RH40 (27,38 mN/m) i Ber 
(36,42 mN/m) [D4] jest większa niż ta składowa napięcia powierzchniowego wody 
(26,85 mN/m). W przypadku napięcia powierzchniowego ogona ELP i RH40 
składowa у” tego napięcia jest tylko nieznacznie większa niż ta dla wody 
(Tab. 1, [D4]). Jak wspomniano powyżej, w cząsteczkach Ber trudno jest rozróżnić 
ogon i głowę, i w związku z tym w monowarstwie adsorpcyjnej na granicy faz 
woda-powietrze przyjmują one równoległą orientację. Wobec tego, minimalne 
napięcie powierzchniowe wodnego roztworu Ber nie może być mniejsze od jej 


napięcia powierzchniowego (46,52 mN/m) (Tab 1. [D4]). 
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Zakładając, ze tylko napięcie powierzchniowe ogona ELP i RH40 decyduje 
o redukcji napięcia powierzchniowego wody, można stwierdzić, że minimalne 
napięcie powierzchniowe roztworu wodnego powinno być znacznie mniejsze niż 
w przypadku roztworu wodnego Ber. 7 analizy różnic pomiędzy napięciem 
powierzchniowym Ber oraz ogona ELP i RH40 wynika, że jeśli na granicy faz woda- 
powietrze powstaje mieszana monowarstwa adsorpcyjna z roztworu zawierającego 
ELP, RH40 i ELP + RH40 z Ber, napięcie powierzchniowe wodnych roztworów tych 
mieszanin przy tym samym stężeniu ELP i RH40 oraz ELP + RH40 powinno być 
większe niż dla roztworu bez Ber. Izotermy napięcia powierzchniowego wodnych 
roztworów ELP, RH40 i ELP + RH40 w obecności i przy braku Ber w roztworze 
(Rys. 2 - 4, [D4]) potwierdzają, że usunięcie cząsteczek ELP lub RH40 z mieszanej 
monowarstwy adsorpcyjnej na granicy faz woda - powietrze przez cząsteczki Ber 
zwiększa napięcie powierzchniowe roztworu. Izotermy te mają prawie taki sam 
kształt i w danej temperaturze istnieje niemal prostoliniowa zależność pomiędzy 
napięciem powierzchniowym roztworu a logarytmem ze stężenia surfaktantów lub 
ich mieszaniny. Wykazano, że podobnie jak w przypadku izoterm ELP, RH40 
i ELP+RH40 [D3], izotermy ELP, RH40 i ich mieszaniny z berberyną można opisać 
funkcją wykładniczą drugiego rzędu, a stałe w równaniu tej funkcji powiązać ze 
składowymi i parametrami napięcia powierzchniowego wszystkich składników 
roztworu. Dodatkowo zmiany stałych А;, А, ё i tą tej funkcji ze zmianą 
temperatury (Rys. S2, [SM4]) są bardziej skomplikowane niż zmiany stałej y? (Rys. 
S2a, [SM4]) i są związane ze zmianami stopnia hydratacji cząsteczek ELP, RH40 
i Ber i odległości pomiędzy tymi cząsteczkami i wodą. Zmiany stopnia hydratacji 
oraz składowych i parametrów napięcia powierzchniowego mieszaniny wpływają 


z kolei na stężenia i skład mieszanej monowarstwy adsorpcyjnej. 


Skład mieszanej monowarstwy adsorpcyjnej utworzonej przez dwu- 
i trójskładnikowe mieszaniny surfaktantów w wielu przypadkach można 
wyznaczyć za pomocą zmodyfikowanych równań Rosena i Rubingha [59, 64, D4]. 
Równania te jednak można zastosować do wyznaczenia ułamka molowego 
mieszanej monowarstwy tylko w ograniczonym składzie mieszaniny w fazie 


objętościowej. Z tego względu w przypadku badanych mieszanin surfaktantów 
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z berberyną określono skład mieszanej monowarstwy tylko przy niewielkim 
stężeniu surfaktantów [D4]. Wykazano, że skład obliczonej w ten sposób mieszanej 
monowarstwy jest zbliżony do tego wyznaczonego na podstawie izotermy napięcia 
powierzchniowego wodnych roztworów poszczególnych składników mieszaniny, 
znacznie różni się od tego w fazie objętościowej i zależy od temperatury (Rys. S2 — 
S5, [SM4]). Interesującym jest, że przy dużym stężeniu ELP, RH40 lub mieszaniny 
ELP + RH40 ułamek molowy Ber w mieszanej monowarstwie jest większy od tego 
w fazie objętościowej [D4]. Sugeruje to, że cząsteczki berberyny adsorbują się na 
granicy faz woda-powietrze łącznie z cząsteczkami surfaktantu w wyniku silnych 
oddziaływań hydrofobowych między cząsteczkami Ber i surfaktantu i cząsteczki 
Ber w tej monowarstwie nie są zorientowane równolegle do granicy faz, lecz 
prostopadle lub pod pewnym kątem. Efektem tego jest wzrost stężenia Ber 


w mieszanej monowarstwie (Rys. 52 — S5, [SM4]). 


Stężenie danego składnika mieszaniny  surfaktantów czy mieszaniny 
surfaktantów z dodatkami można wyznaczyć pośrednio m.in. za pomocą izotermy 
adsorpcji Gibbsa i Frumkina [104]. Niestety w badanym zakresie stężeń mieszanin 
odpowiadającym nasyconej mieszanej monowarstwie na granicy faz woda- 
powietrze za pomocą równania izotermy adsorpcji Gibbsa nie udało się wyznaczyć 
maksymalnego stężenia mieszaniny ELP, RH40 i ELP + RH40 w obecności Ber. 
Należy podkreślić, że słaba rozpuszczalność berberyny w wodzie uniemożliwia 
pomiar napięcia powierzchniowego wodnych roztworów mieszaniny surfaktantu 
i Ber w szerokim zakresie stężenia surfaktantu. Z tego powodu stężenie danego 
surfaktantu w badanych pojedynczych i mieszanych monowarstwach na granicy 
faz woda — powietrze wyznaczono za pomocą równania Frumkina, stosując do 
obliczeń wkład poszczególnych składników roztworu w redukcję napięcia 
powierzchniowego [105, D4]. Wykazano, że dla wszystkich badanych mieszanin 
stężenie poszczególnych składników w mieszanej monowarstwie na granicy faz 
woda-powietrze jest mniejsze niż poszczególnych składników osobno (Rys. S6 — 
S8, [SM4]). Ponadto suma stężeń ELP + Ber, RH40 + Ber czy ELP + RH40 + Ber 
w mieszanej monowarstwie adsorpcyjnej jest mniejsza niż w przypadku 


surfaktantów lub ich mieszaniny bez berberyny (Rys. 56 — S8, [SM4]). 
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Prawdopodobną przyczyną tego faktu jest adsorpcja Ber łącznie z cząsteczkami 
surfaktantu. Cząsteczki berberyny mogą łączyć się z ogonem cząsteczki surfaktantu 
zmniejszając jego stopień hydratacji. W konsekwencji powstaje kompleks Ber - 
surfaktant o mniejszej tendencji do adsorpcji na granicy faz woda-powietrze niż 
pojedyncza cząsteczka ELP czy RH40. Z drugiej jednak strony zaadsorbowany 
kompleks cząsteczka berberyny + ogon surfaktantu w mieszanej monowarstwie 
adsorpcyjnej może być skierowany w stronę fazy wodnej zmieniając jej napięcie 
powierzchniowe. Ze względu na budowę cząsteczki Ber, która posiada ładunek 
dodatni, mogą pojawiać się także odpychające oddziaływania elektrostatyczne 
pomiędzy kompleksami Ber + ogon surfaktantu. Tego typu oddziaływania mogą 
zmniejszać upakowanie mieszanej monowarstwy adsorpcyjnej, zmniejszając 
jednocześnie stężenie surfaktantu w porównaniu do jego stężenia w roztworze bez 


berberyny. 


Standardowa swobodna energia Gibbsa (AG2,,), standardowa entalpia (AH24,) 


i entropia procesu adsorpcji ( 450,,) są przydatne do określenia tendencji 
surfaktantów do adsorpcji na różnych granicach faz i do poznania przyczyny tej 
tendencji [D4]. W przypadku badanych mieszanin Ber + surfaktant wartości АС, 
dla ELP, RH40 i Ber, w mieszaninach ELP + Ber, RH40 + Ber, ELP + RH40 + Ber 
obliczono za pomocą zaproponowanego równania przy założeniu, że 
współczynniki aktywności są bliskie jedności [D4] i porównano z tymi obliczonymi 
na podstawie zmodyfikowanego równania Langmuira (Rys. S9, [SM4]). Otrzymane 
wyniki wskazują, że w zakresie stężeń surfaktantów, w których występują one 
w roztworze w postaci monomerycznej, wartości АС: obliczone 
z zaproponowanego równania są praktycznie stałe i nie zmieniają się ze wzrostem 
T (Tab. 4, [D4]). Przy stężeniu surfaktantów większym niż CMC wartości АС, 
rosną w funkcji stężenia. Należy jednak podkreślić, że cząsteczki surfaktantów 
adsorbują się na granicy faz woda-powietrze tylko w postaci monomerycznej, 
a przy stężeniu większym niż CMC ich stężenie w formie monomerycznej jest stałe. 
W związku z tym obliczone w ten sposób wartości AG2,, dla stężeń powyżej CMC 


nie są rzeczywiste. Zatem na podstawie wartości AG?,,, ale także ДНО; można 
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stwierdzić, ze obecność Ber zmniejsza tendencję ELP i RH40 do adsorpcji na granicy 


faz woda - powietrze ze względu na adsorpcję kompleksów Ber + surfaktant. 


Obok adsorpcyjnych, również właściwości agregacyjne badanych układów 
Ber + surfaktant mogą mieć bardzo istotne znaczenie np. w procesie solubilizacji. 
W związku z tym zrodziło się pytanie, czy berberyna jest obecna w micelach ELP, 
RH40 i ich mieszaniny. W celu znalezienia odpowiedzi na to pytanie w pierwszej 
kolejności za pomocą pomiarów napięcia powierzchniowego, przewodnictwa, 
gęstości, widm emisyjnych berberyny [D1, D2, D3, D4] wyznaczono krytyczne 
stężenie micelizacji badanych mieszanin. Wykazano, że wartości CMC ELP, RH40 
i ich mieszaniny w obecności Ber są większe od tych bez Ber (Tab. 2, [D4]). Należy 
pamiętać, że tworzenie micel w wodnym roztworze surfaktantów występuje na 
skutek oddziaływań hydrofobowych pomiędzy ogonami surfaktantów poprzez 
fazę wodną. Wartości tych oddziaływań są dodatnie, w przeciwieństwie do 
oddziaływań pomiędzy głowami surfaktantów, które są ujemne niezależnie od 
tego, czy micele tworzą niejonowe, czy jonowe związki powierzchniowo czynne. 
Zgodnie z równaniami zaproponowanymi przez Jańczuka i wsp. [106, 107] siła 
oddziaływań hydrofobowych zależy od napięcia międzyfazowego woda - ogon 
i powierzchni kontaktu ogona surfaktantu. Biorąc pod uwagę podejście van Ossa 
i wsp. oraz składowe i parametry napięcia powierzchniowego wody, berberyny 
i surfaktantów obliczono napięcie międzyfazowe woda - ogon surfaktantu (wr), 
woda-Ber (Ywg ) i Ber-ogon surfaktantu ( урт ) [D4]. Wartości te wynoszą 
odpowiednio 46,0; 6,0 i 10,7 mN/m. Na ich podstawie obliczono pracę adhezji 
ogona surfaktantu do ogona surfaktantu, ogona surfaktantu do Ber i Ber do ogona 
surfaktantu przez fazę wodną. Wykazano, że praca adhezji ogona surfaktantu do 
ogona surfaktantu przez fazę wodną wynosi 2ywr= 92 mJ/m?, a Ber do Ber jest 
równa Ywg= 12 mJ/m*. Praca adhezji w układzie Ber - woda- ogon surfaktantu 


(WET) wynosi 41,3 mJ/ m? [D4]. 


Ze względu na fakt, że napięcie międzyfazowe surfaktant - woda dla ELP i RH40 
jest ujemne i bliskie -18 mJ/m?, ich całkowita praca adhezji poprzez fazę wodną jest 


bliska 56 mJ/m?. Nie istnieje więc istotna różnica w tendencji do kontaktu 
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surfaktantu z drugą cząsteczką surfaktantu i surfaktantu z cząsteczką berberyny 
poprzez fazę wodną. Jednakże prawdopodobieństwo wiązania poprzez fazę 
wodną kompleksów Ber - związek powierzchniowo czynny jest mniejsze niż 
cząsteczek samego surfaktantu i prawdopodobnie z tego powodu CMC mieszanin 
ELP, RH40 i ELP + RH40 z Ber jest większe niż bez Ber (Tab. 2, [D4]). Biorąc pod 
uwagę powyższe obliczenia można założyć, że proces micelizacji mieszaniny 
badanych surfaktantów zachodzi w wyniku nie tylko łączenia się cząsteczek 
surfaktantów, ale także surfaktantów z cząsteczkami Ber. Wydaje się zatem, że 
obecność Ber w micelach ELP, RH40 i ich mieszaninie nie jest spowodowana 
adsorpcją czy też penetracją cząsteczek Ber w micelach, ale raczej jest wynikiem 


wspólnej agregacji (Rys. 5). 


W przypadku mieszanych micel ELP + Ber i RH40 + Ber trudno jest określić 
ułamek molowy związku powierzchniowo czynnego i berberyny w micelach na 
podstawie ich wartości CMC. W związku z tym ułamki te wyznaczono za pomocą 
zmodyfikowanej koncepcji Hua i Rosena [59, 64] i wykazano, że ułamek molowy 
Ber nie różni się istotnie od 0,5. Świadczy to o tym, że oddziaływania kompleksów 
berberyna + surfaktant poprzez fazę wodną odgrywają wiodącą rolę w procesie 
agregacji. Stwierdzono także, że ułamek molowy poszczególnych składników 
mieszaniny ELP + RH40 + Ber zmienia się w niewielkim stopniu w funkcji 
temperatury (Tab. 3, [D4]). Należy jednak podkreślić, że zmiana wartości ułamków 
molowych składników mieszaniny, jak i samej wartości CMC w funkcji 
temperatury wynika z jednej strony ze zmiany energii kinetycznej, a z drugiej 
strony ze zmiany stopnia hydratacji, zwłaszcza części hydrofilowej surfaktantu. 


Dodatkowo wpływać może także zmiana konfiguracji cząsteczek surfaktantu. 


Wpływ wymienionych czynników na skład i wielkość badanych micel znajduje 
także odzwierciedlenie w wartościach parametru oddziaływań 
miedzyczasteczkowych (£^). Ponieważ Ber nie tworzy samoistnie micel, nie było 
możliwe określenie parametru oddziaływań pomiędzy cząsteczkami Ber 
i surfaktantów w micelach mieszanin ELP + Ber i RH40 + Ber bezpośrednio 


z równania zaproponowanego przez Rosena i Hua [59]. W związku z tym założono, 
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ze kompleks ELP + Ber i RH40 + Ber to pojedynczy zwiazek i wyznaczono za 
pomocą tej koncepcji BM (Ber + ELP) - RH40 i (Ber + RH40) - ELP w mieszanych 
micelach ELP + RH40 + Ber (Tab. 3, [D4]). Z przeprowadzonych obliczeń wynika, 
że parametr 8“ obliczony na podstawie wartości CMC wyznaczonego z pomiarów 
napięcia powierzchniowego jest ujemny dla kompleksu RH40 + Ber 
w każdej temperaturze( —0,4159 w 293 K, —1,5256 w 303 K i —1,1598 w 318 K), ale 
dla kompleksu ELP + Ber tylko w T = 293 K (—5,0834). Podobny przebieg zmian 
zaobserwowano w przypadku wartości nadmiarowej energii Gibbsa tworzenia 
miceli przypadającej na mol mieszaniny surfaktantów (G™) (Tab. 3, [D4]). Ujemna 
wartość parametru ВМ sugeruje, że istnieje efekt synergetyczny w procesie 
tworzenia micel. Niestety jednak, w przypadku badanych mieszanin niemożliwe 
było potwierdzenie tej sugestii za pomocą drugiego warunku istnienia 


synergetyzmu określonego przez Rosena i Hua [D4]. 
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7. WNIOSKI 


Ма podstawie wyników uzyskanych z przeprowadzonych badań oraz ich 


termodynamicznej analizy stwierdzono m.in., że: 


1. Berberyna w niewielkim stopniu redukuje napięcie powierzchniowe wody a jej 
maksymalne nadmiarowe stężenie Gibbsa w warstwie powierzchniowej 


w temperaturze 293K jest zdecydowanie mniejsze od tego dla ELP i RH40. 


2. Napięcie powierzchniowe berberyny wynika zarówno ze składowej Lifshitza- 
van der Waalsa, jak i kwasowo - zasadowej, jednak udział składowej Lifshitza-van 
der Waalsa jest dużo większy wpływając na słabą rozpuszczalność berberyny 


w wodzie. 


3. Adsorpcja cząsteczek ELP i RH40 na granicy faz woda-powietrze powoduje 


redukcję składowej Lifshitza-van der Waalsa napięcia powierzchniowego wody. 


4. Składowe Lifshitza-van der Waalsa napięcia powierzchniowego Ber oraz ogona 
ELP i RH40, a także stopień hydratacji cząsteczek Ber oraz ogona i głowy 
cząsteczek ELP i RH40 w największym stopniu decydują o stężeniu i składzie 


nasyconej monowarstwy adsorpcyjnej na granicy faz woda-powietrze. 


5. Stężenie poszczególnych składników w mieszanej monowarstwie adsorpcyjnej 
na granicy faz woda-powietrze dla wszystkich badanych mieszanin z Ber jest 


mniejsze od tego dla roztworów pojedynczych związków. 


6. Przy dużych stężeniach ELP, RH40 lub mieszaniny ELP + RH40 ułamek molowy 


Ber w mieszanej monowarstwie jest większy od tego w fazie objętościowej. 


7. Stężenie berberyny w mieszanych micelach jest większe niż w fazie 
objętościowej, co wskazuje, że tendencja Ber do solubilizacji w micelach ELP, 
RH40 i ELP + RH40 jest większa niż jej tendencja do adsorpcji na granicy faz woda- 


powietrze. 


8. Skład mieszanej monowarstwy adsorpcyjnej ELP +RH40 na granicy faz woda- 


powietrze można określić na podstawie wkładu poszczególnych składników 
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mieszaniny w redukcje napiecia powierzchniowego wody i jest on bliski temu 


obliczonemu za pomocą koncepcji Rosena i Hua. 


9. Istnieje efekt synergetyczny w procesie tworzenia mieszanych micel ELP+RH40. 
Efekt ten wykazano za pomocą wartości parametru oddziaływań 
międzycząsteczkowych obliczonych z pomiarów napięcia powierzchniowego 
roztworów, a także parametrów akustycznych wyznaczonych z pomiarów 


gęstości i prędkości dźwięku. 


10. Biorąc pod uwagę wartości nadmiarowego stężenia powierzchniowego Gibbsa 
Ber w mieszaninie ELP i RH40, współczynniki aktywności tych związków 
w warstwie adsorpcyjnej na granicy faz woda-powietrze i w micelach obliczone 
metodą Rosena i Hua i wkład Ber, ELP i RH40 w redukcję napięcia 
powierzchniowego wody, można za pomocą zaproponowanych równań określić 


energię Gibbsa oddziaływań Ber z ELP i RH40 poprzez fazę wodną. 
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Abstract Measurements of density, speed of sound, and 
surface tension as well as viscosity of aqueous Kolliphor® 
ELP (ELP) and Kolliphor® RH 40 (КН40) solutions ав 
well as binary mixtures at different surfactant mole frac- 
tions were made at 293 K to investigate their aggregation 
behavior. The free volume, internal pressure, and molar 
cohesive energy were calculated and compared to the spe- 
cific acoustic impedance and intermolecular free length to 
obtain qualitative information about the character of inter- 
actions between the surfactant molecules in the mixture 
through the water phase. 
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Introduction 


Micellar systems are widely applied in chemical, biochemi- 
cal, and pharmaceutical areas as solubilizing agents and also 
constitute a basis for drug delivery or solvent systems in sep- 
aration processes. Moreover, some aqueous solutions with 
nonionic surfactants can be applied in extraction processes 
being of significant importance in biotechnological product 
purification (Ingram et al., 2012; Lam et al., 2005). In phar- 
macy, nonionic Kolliphors (known before as Cremophors) 
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find application in lipid-based delivery aimed at improving 
poorly water-soluble drugs (Berthelsen et al, 2015; 
Gelderblom et al, 2001). The triricinoleate ester of 
ethoxylated glycerol is the main constituent of Kolliphors. 
Polyethylene glycol ricinoleate and the corresponding free 
glycols are the other constituents (Christiansen et al., 2010). 
Kolliphor EL (EL), the most common of all Kolliphors, can 
be used as a solubilization vehicle for some hydrophobic 
drugs including cyclosporin A, diazepam, propofol, and pac- 
litaxel (Bakonyia et al., 2018; Bali et al., 2011; Chen et al., 
2008; Constantinides and Wasan, 2007; Cuiné et al., 2008; 
Gao et al., 2003; Sadurni et al., 2005; Seelig and Gerebtzoff, 
2006; Wang et al., 2004). Moreover, mixed micelles of non- 
ionic surfactants are common in industry (Rosen, 2004). 
Interactions between individual surfactants can result in 
lower critical micelle concentration (CMC) of the binary sur- 
factant mixture vs. either individual surfactant or a thermo- 
dynamically more stable one. Thus, better knowledge about 
the origin of nonideal and synergistic behavior contributes to 
more efficient surfactant mixture design so as to reduce the 
amount of surfactant in a formulation. Determination of vol- 
umetric and viscometric properties of aqueous solutions of 
binary Kolliphor? ELP (ELP), purified grade of EL, and 
Kolliphor? RH 40 mixtures from the measurements of den- 
sity, viscosity, speed of sound, and surface tension at 203 K 
was the aim of this article. These data were used to obtain 
information about the nature of the solute-solute, solute—sol- 
vent, and solvent-solvent interactions. 


Experimental 


Kolliphor® ELP (ELP) (Cremophor® ELP, Polyoxyl 
35 Hydrogenated Castor oil, and Polyoxyl-35 castor oil) 
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Fig. 1 A plot of the values of p of the aqueous solutions of RH40, ELP, and their binary mixtures at the mole fraction of ELP equal to 0.2, 0.4, 0.6, 
and 0.8 at T = 293 K vs. log C as well as the values of p of the studied solutions at C= 10 2,2 x 107°, 5 x 107°, 8 x 107°, and 1072 M vs. a 


and Kolliphor® RH 40 (RH40) (Cremophor® RH 
40, Macrogolglycerol hydroxystearate, PEG-40 castor oil, 
and Polyoxyl 40 hydrogenated castor oil) were provided by 
Sigma (Poznań, Poland) and used without further purifica- 
tion. RH40 is mainly composed of trihydroxystearate ester 
of ethoxylated glycerol with a smaller amount of polyethyl- 
ene glycol trihydroxystearate as well as the corresponding 
free glycols. ELP is a nonuniform product without a single 
well-defined structure (Shaukat, 2011). Tricinoleate esters 
of ethoxylated glycerols, polyethylene glycol ricinooleates 
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as well as their derivatives with  polyethoxylated 
12-hydroxy ricinoleic acid residues are the main compo- 
nents of ELP. Thus, ELP can be treated as a mixture. Fur- 
thermore, as the double bond at C-9 in ELP molecules is 
cis and C-12 atoms of the fatty acid chains are achiral or 
racemic, the proportions of both (R) and (S) enantiomers 
are equal (Shaukat, 2011). Doubly distilled and deionized 
water used for preparation of the aqueous solutions of indi- 
vidual surfactants and their mixtures with the bulk phase 


mole fraction of ELP (a) equal to 0.2, 0.4, 0.6, and 0.8 at 


log C 


Fig. 2 A plot of the values of u of the aqueous solutions of RH40, ELP, and their binary mixtures at the mole fraction of ELP equal to 0.2, 0.4, 0.6, 
and 0.8 at T = 293 K vs. log C as well as the values of u of the studied solutions at C = 107°, 2 х 107°, 5 x 107, 8 x 107°, and 10 ^ M vs. а 
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Fig. 3 A plot of the values of 7 of the aqueous solutions of RH40, ELP, and their binary mixtures at the mole fraction of ELP equal to 0.2, 0.4, 
0.6, and 0.8 at T = 293 K vs. log C as well as the values of у of the studied solutions at C = 107°, 2 x 107, 5 x 107, 8 x 107°, and 10^ M vs. а 


the concentration from 1079 to 10? M was obtained from 
a Destamat Bi18E distiller (Heraeus, Germany). 

Speed of sound and densities of various solutions were 
measured using a vibrating tube densitometer and the speed 
of sound analyzer (Anton Paar DSA 5000M) equipped with 
the automatic viscosity correction and two integrated Pt 
100 thermometers. Temperature was maintained to 0.01 K 
using a temperature controller. The analyzer was calibrated 


210°" (Pa) 


with triple-distilled water and dry air. The standard density 
measurement uncertainties were +2 x 1077 kg m ? and 
those of speed of sound were +0.1 m 877. An Anton Paar 
viscometer (AMVn) at T = 293 K + 0.01 K with a preci- 
sion of 0.0001 mPa s and uncertainty 0.396 was used for 
viscosity measurements of aqueous solutions. 

Surface tension measurements were made with a Krüss 
K100C tensiometer under atmospheric pressure using the 
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Fig. 4 A plot of the values of z; of the aqueous solutions of RH40, ELP and their binary mixtures at the mole fraction of ELP equal to 0.2, 0.4, 
0.6, and 0.8 at T = 293 K vs. log C as well as the values of z, of the studied solutions at C = 8 x 1077 and 10 ? M vs. a 
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internal pressure is made from the following equation 
(Chauhan et al., 2013b; Kumar et al., 2012): 


ктк)? 
A М0/6 


(1) 


where р is the density, и is the speed of sound, 7 is the vis- 
cosity, R is the gas constant, T is the temperature, M is the 
the molar mass, b is the the packing factor that is assumed 
to be equal to 2 in the liquid system and K’ = 4.28 x 10° is 
the constant regardless of the liquid. An increase of z; is 
found with increasing concentration of both individual sur- 
factants and their mixtures (Fig. 4) indicating that inter- 
molecular interactions increase as a result of formation of 
solvent molecule aggregates around the solute affecting the 
structural arrangement of the solvent system. On the other 
hand at high concentration of surfactants in the solution, it 
can be associated with the changes of the size and/or shape 
of the Kolliphor micelles. 

The surface tension (y; y) of the surfactant aqueous solu- 
tions (Fig. 5) was used to determine the CMC, Псмс (the 
difference between the surface tension of the solvent, yo, 
and the solute at the CMC, усмс), АС"... (the standard 
Gibbs energy of micellization) (Szymczyk, 2011), and 
AG; 


min 


Ay (mN m) 
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Fig. 5 A plot of the values of yry of the aqueous solutions of RH40, 
ELP, and their binary mixtures at the mole fraction of ELP equal to 
0.2, 0.4, 0.6, and 0.8 at T = 293 K vs. log C 


ring method. Measurements of the surface tension of pure 
water at 293 + 0.1 K were performed to calibrate the tensi- 
ometer and to check the cleanliness of the glassware. In all 
cases, more than 10 measurements were made and the stan- 
dard deviation did not exceed +0.2 mN m. 


(the free energy at the water-air interface) (Sugihara 


et al., 2003), yo/ycmc as well as ХИ (the mole fraction of 
surfactant in the mixed micelle), p~“ (molecular interaction 
parameter), ЇГ and ЇГ (the activity coefficients of surfac- 
tants 1 and 2 in the mixed micelle), p" (the excess Gibbs 
energy of the mixed micelle) and eM (the molecule 
exchange energy in the mixed micelle) based on the 
Rubingh and Rosen theory (Hua and Rosen, 1982; Rosen, 


Results and Discussion 


The internal pressure, л;, was calculated from the values of 
density (Fig. 1), speed of sound (Fig. 2), and viscosity 


(Fig. 3) of ELP, RH40 aqueous solutions, and their binary 
mixtures. The internal pressure is a measure of total disper- 
sion, ionic, and dipolar forces for a liquid and includes all 
attractive and repulsive forces. The calculation of the 


2004; Rubingh, 1979; Wang et al., 2005). The calculated 
parameters are presented in Table 1. Values of CMC, 
Temc.AG?,., АСЗ... Й, gM, and еМ for the mixture with 


the mole fraction of ELP (a) in the bulk phase equal to 0.8 


Table 1 Values of CMC, Псмс, AGO, AGE, Yolyemc, ХМ, В“, fM, fM. 9M. and eM for the aqueous solutions of RH40, ELP, and their binary 
mixtures at the mole fraction of ELP (a) equal to 0.2, 0.4, 0.6, and 0.8 at T = 293 K 


RH40 а-02 a=0.4 a=0.6 a=0.8 ELP 
CMC (mol dm *) 6.64 x 1075 3.94 x 1075 2.88 x 1077 2.45 x 1077 1.92 х 10^? 2.14 x 1079 
Iemc(mN m^!) 31.64 31.70 32.55 33.16 34.8 32.91 
AG? (kJ mol!) —23.44 —24.71 —25.47 —25.86 —26.46 —26.19 
АС in (KJ mol!) 21.56 19.16 15.30 13.82 11.99 12.59 
Yol cMc 1.77 1.77 1.81 1.84 1.92 1.82 
ху = 0.452 0.620 0.79 0.895 = 

pM - —0.694 —0.980 —0.342 —2.097 = 

ЎЎ = 0.867 0.686 0.807 0.186 = 

ję — 0.812 0.868 0.985 0.977 ES 

g" (kJ шог) = —0.434 —0.700 —0.422 —3.338 = 
eM(kJ mol!) — —0.282 —0.398 —0.139 —0.851 — 
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Fig. 6 A plot of the values of MCE of the aqueous solutions of 
RH40, ELP and their binary mixtures at the mole fraction of ELP 
equal to 0.2, 0.4, 0.6, and 0.8 at T = 293 K vs. log C as well as the 
values of MCE of the studied solutions at C = 107°, 2 x 107°, 
5 x 107°, 8 x 107°, and 10^ M vs. а 


clearly indicate that at the studied temperature, synergism 
in mixed micelle formation causes a significantly stronger 
driving force for micellization. This mixture shows the 
smallest value of ДС, which is regarded as the free 
energy change accompanied by the transition from the bulk 
phase to the surface phase for all solution components and 
the biggest value of the yo/rcwc parameter, which is con- 
nected with the solvophobic effect (Rodriguez et al., 2008; 


Sachin et al., 2019). On the other hand, the measure of 


a 


r (ps) 


attraction between different solution components is the 
molar cohesive energy (MCE), which can be calculated 
from the relation: 


МСЕ= mV, (2) 


where Vm is the molar volume (Chauhan et al., 2013a). 
Figure 6 shows that MCE values increase for all solutions 
when C is higher than 107° M due to formation of micelles. 
Values of MCE for RH40 molecules are larger consistent 
with smaller viscous relaxation time (7) evaluated based on 
(Chauchan et al., 2018): 
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(3) 


Figure 7 shows that the т values increase with the 
increasing surfactant concentration and ELP mole fraction 
at the studied temperature. Smaller values of 7 are found 
for RH40 solutions as the viscous relaxation time is propor- 
tional to the solution viscosity (Fig. 3). However, at high 
concentrations, much larger than CMC determined from the 
surface tension measurements, a linear increase of tvalues 
with the mole fraction of ELP is not observed thus 
suggesting the presence of synergetic intermolecular inter- 
actions in the molecules of ELP and КН40 in the mixture. 

This observation is consistent with the acoustic imped- 
ance (Z), which are possible to be calculated from the rela- 
tion (Hua and Rosen, 1982): 


Z-up (4) 


r (ps) 


log C 


Fig. 7 A plot of the values of t of the aqueous solutions of RH40, ELP and their binary mixtures at the mole fraction of ELP equal to 0.2, 0.4, 
0.6, and 0.8 at T = 293 K vs. log C as well as the values of t of the studied solutions at C = 8 x 107? and 10 2 M vs. а 
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Table 2 Values of ks, Ly, Z, Va, and R, for the aqueous solutions of RH40, ELP, and their binary mixtures at the mole fraction of ELP (a) equal 
to 0.2, 0.4, 0.6, and 0.8 at T = 293 K 


С ks (107? m? №!) 

RH40 1079 4.5582 
2x 1075 4.5582 
5x 1075 4.5582 
8 x 1079 4.5581 
1079 4.5581 
2х 107 4.5581 
5х107 4.5579 
8x 1077 4.5578 
107* 4.5576 
2x 1077 4.5572 
5x 10-4 4.5557 
8 x 10+ 4.5545 
1077 4.5527 
2х 107 4.5439 
5x10? 4.5245 
8 x 107° 4.5044 
107? 4.4911 
a=0.2 1079 4.5572 
2x 1075 4.5572 
5x 1079 4.5572 
8 х 1075 4.5572 
1077 4.5572 
2x10? 4.5571 
5x10? 4.5571 
8 x 10> 4.5570 
107* 4.5568 
2x 107 4.5563 
5x107% 4.5547 
8 x 1074 4.5527 
107° 4.5497 
2x 107 4.5378 
5x10? 4.5168 
8 x 107° 4.4987 
1072 4.4855 
a=0.4 1079 4.5560 
2x 1075 4.5560 
5x 1079 4.5559 
8 х 1079 4.5558 
1077 4.5558 
2x 1077 4.5557 
5х107 4.5556 
8 x 1077 4.5554 
107* 4.5551 
2x 10-4 4.5543 
5x10* 4.5524 
8 x 10-4 4.5502 
10? 4.5478 
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Le (107? m) 


0.435003 
0.435003 
0.435003 
0.435003 
0.435002 
0.435000 
0.434993 
0.434986 
0.434977 
0.434958 
0.434888 
0.434828 
0.434743 
0.434321 
0.433392 
0.432429 
0.431792 
0.434958 
0.434958 
0.434958 
0.434957 
0.434955 
0.434954 
0.434952 
0.434948 
0.434937 
0.434914 
0.434838 
0.434742 
0.434599 
0.434032 
0.433025 
0.432158 
0.431521 
0.434899 
0.434899 
0.434895 
0.434891 
0.434891 
0.434886 
0.434880 
0.434872 
0.434858 
0.434817 
0.434726 
0.434621 
0.434511 


Z (10° kgm * s 1) 


1.4798 
1.4798 
1.4798 
1.4798 
1.4798 
1.4799 
1.4799 
1.4799 
1.4799 
1.4800 
1.4803 
1.4806 
1.4809 
1.4826 
1.4864 
1.4903 
1.4929 
1.4800 
1.4800 
1.4800 
1.4800 
1.4800 
1.4800 
1.4800 
1.4800 
1.4801 
1.4802 
1.4805 
1.4809 
1.4814 
1.4836 
1.4877 
1.4913 
1.4939 
1.4802 
1.4802 
1.4802 
1.4802 
1.4802 
1.4802 
1.4803 
1.4803 
1.4804 
1.4805 
1.4809 
1.4813 
1.4817 


V, (107° m? шог) 


29.4256 
29.4256 
29.4256 
29.4255 
29.4255 
29.4246 
29.4197 
29.4149 
29.4075 
29.3964 
29.3496 
29.3113 
29.2476 
28.9309 
28.2686 
27.5750 
27.1150 
29.3879 
29.3879 
29.3879 
29.3867 
29.3859 
29.3852 
29.3845 
29.3823 
29.3743 
29.3588 
29.3071 
29.2385 
29.1255 
28.6853 
27.9577 
27.3488 
26.8910 
29.3385 
29.3385 
29.3355 
29.3318 
29.3318 
29.3280 
29.3247 
29.3192 
29.3091 
29.2782 
29.2154 
29.1403 
29.0556 


R, (1074 т!9 5713 mol!) 


45.6813 
45.6813 
45.6812 
45.6812 
45.6811 
45.6809 
45.6807 
45.6806 
45.6807 
45.6798 
45.6777 
45.6753 
45.6753 
45.6751 
45.6624 
45.6513 
45.6437 
45.6826 
45.6826 
45.6826 
45.6827 
45.6826 
45.6824 
45.6821 
45.6819 
45.6820 
45.6812 
45.6793 
45.6781 
45.6800 
45.6846 
45.6739 
45.6582 
45.6499 
45.6844 
45.6843 
45.6844 
45.6845 
45.6845 
45.6845 
45.6842 
45.6839 
45.6839 
45.6838 
45.6820 
45.6807 
45.6813 
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Table 2 Continued 


С ks (107? m? №!) Le (1079 m) Z (106 kgm ? s^ ) V, (1075 m? шог) R, (1075 m 95 5715 mot!) 
2x 1077 4.5359 0.433940 1.4839 28.6161 45.6848 
5x103 4.5148 0.432930 1.4881 27.8934 45.6714 
8 x 10? 4.4961 0.432030 1.4919 27.2644 45.6540 

1072 4.4817 0.431340 1.4947 26.7744 45.6430 
a=0.6 1079 4.5558 0.434893 1.4802 29.3325 45.6848 
2x 1075 4.5557 0.434886 1.4802 29.3271 45.6849 
5х 1079 4.5555 0.434877 1.4803 29.3196 45.6852 
8 x 1079 4.5554 0.434873 1.4803 29.3157 45.6854 
1077 4.5554 0.434871 1.4803 29.3140 45.6854 
2x 107 4.5552 0.434862 1.4803 29.3076 45.6854 
5x10? 4.5550 0.434851 1.4804 29.2997 45.6853 
8 x 1075 4.5547 0.434838 1.4804 29.2895 45.6854 
107* 4.5544 0.434824 1.4805 29.2792 45.6854 
2х 1077 4.5535 0.434779 1.4806 29.2460 45.6852 
5x 107 4.5516 0.434689 1.4810 29.1829 45.6835 
8x107* 4.5493 0.434580 1.4815 29.1050 45.6822 
1077 4.5469 0.434464 1.4819 29.0155 45.6829 
2x 107 4.5347 0.433880 1.4841 28.5658 45.6864 
5x 107 4.5132 0.432854 1.4884 27.8322 45.6727 
8 x 10? 4.4947 0.431964 1.4921 27.2155 45.6536 
1072 4.4800 0.431257 1.4950 26.7043 45.6454 
a=0.8 1079 4.5566 0.434930 1.4801 29.3630 45.6838 
2x 1075 4.5566 0.434929 1.4801 29.3630 45.6838 
5x 1079 4.5566 0.434929 1.4801 29.3629 45.6837 
8 x 1079 4.5566 0.434929 1.4801 29.3630 45.6838 
1077 4.5566 0.434929 1.4801 29.3629 45.6837 
2x 107 4.5562 0.434911 1.4802 29.3487 45.6839 
5x10? 4.5554 0.434870 1.4803 29.3148 45.6849 
8x 10> 4.5548 0.434842 1.4804 29.2921 45.6855 
107* 4.5545 0.434829 1.4804 29.2820 45.6857 
2x107* 4.5532 0.434768 1.4807 29.2328 45.6869 
5x 10-4 4.5507 0.434647 1.4811 29.1414 45.6870 
8x 1074 4.5475 0.434493 1.4817 29.0215 45.6883 
1077 4.5446 0.434357 1.4822 28.9132 45.6903 
2x 107 4.5318 0.433745 1.4846 28.4441 45.6931 
5x103 4.5115 0.432773 1.4886 27.7658 45.6744 
8 x 10? 4.4943 0.431944 1.4921 27.1790 45.6610 
1072 4.4797 0.431242 1.4950 26.6683 45.6538 
ELP 1079 4.5580 0.434998 1.4799 29.4249 45.6801 
2x 1075 4.5580 0.434998 1.4799 29.4249 45.6801 
5x 1079 4.5580 0.434998 1.4799 29.4249 45.6801 
8x 1075 4.5580 0.434998 1.4799 29.4249 45.6801 
1077 4.5580 0.434997 1.4799 29.4248 45.6801 
2x 1077 4.5580 0.434997 1.4799 29.4248 45.6800 
5x10? 4.5579 0.434991 1.4799 29.4219 45.6795 
8x 10> 4.5578 0.434984 1.4799 29.4167 45.6795 
1074 4.5577 0.434979 1.4800 29.4139 45.6791 
2x 1077 4.5573 0.434962 1.4800 29.4053 45.6777 
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С ks (107? m? №!) Le (10779 m) 
5 x 107+ 4.5563 0.434916 
8 x 1074 4.5539 0.434800 
107° 4.5523 0.434725 
2x 107 4.5440 0.434326 
5x10? 4.5226 0.433305 
8x10? 4.5026 0.432344 
1072 4.4873 0.431609 


V; (m3 тог) 


Z (10° kgm ? s7!) 


Va (1076 m? mol!) К, (1075 m s 15 mol!) 


1.4803 29.3801 45.6744 
1.4807 29.2971 45.6732 
1.4810 29.2427 45.6725 
1.4826 28.9527 45.6690 
1.4869 28.2297 45.6535 
1.4909 27.5557 45.6362 
1.4939 27.0216 45.6286 


V, (m? mol) 


log C 


Fig. 8 A plot of the values of V; of the aqueous solutions of RH40, ELP, and their binary mixtures at the mole fraction of ELP equal to 0.2, 0.4, 
0.6, and 0.8 at T = 293 K vs. log C as well as the values of V; of the studied solutions at C = 8 x 10? and 107° M vs. а 


The acoustic impedance probes cohesive forces between the 
surfactant molecules in solution enabling explanation of 
molecular interactions. Table 2 shows that Zvalues of the aque- 
ous solutions of surfactant mixtures are greater than those for a 
single surfactant at a given concentration. Acoustic impedance, 
density, and speed of sound increase with increasing surfactant 
concentration in mixtures indicating possible strong molecular 
interactions among the mixture components. EJ.A. Suys 
showed that the diameter of EL aggregates measured by DLS 
and SAXS was 12-15 nm with little micellar growth observed 
on increasing EL concentration (Suys et al., 2019). For charac- 
terization of aggregation, the distance between the surfaces of 
two molecules i.e. molecular free length, L; can be determined. 
І; values depend on both intermolecular and intramolecular 
interactions among the solution components and can be calcu- 
lated from (Chauhan and Sharma, 2015): 


Li=Kyk; (5) 
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where K = [(93.875 + 0.375)T x 107%] and ks is the isen- 
tropic compressibility possible to be determined from the 
speed of sound (u) and density (p) from the Newton- 
Laplace equation (Kiełek and Marczak, 2010): 


На. 
"= 


(6) 


KS 


Values of Lp and ks were calculated based on those of u and 
p of the aqueous solutions of ELP, RH40, and their mixtures 
(Figs 1 and 2) and are presented in Table 2 along with the 
values of available volume (V,) and molar sound velocity (R,) 
(Sannaningannavar et al., 2013), which were calculated from: 


l= 
ү, = == ŻY (7) 
R, = Vnu!’ (8) 
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Fig. 9 A plot of the values of m, of the aqueous solutions of RH40, ELP, and their binary mixtures at the mole fraction of ELP equal to 0.2, 0.4, 
0.6, and 0.8 at C from 2 x 107° to 10 ? M and T = 293 K vs. log C as well as the values of л, of the studied solutions at C = 10 ? M vs.a 


where uœ = 1600 ms ! and V4, is the molar volume. Lr 


values decrease significantly at the value of C higher than 
1077 M for all solutions indicating closer molecular pack- 
ing and is a result of specific intermolecular interactions 
among the molecules of surfactants in micelles. The syner- 
getic effect between different surfactant molecules is 
evidenced by the smaller values of Lẹ at a given concentra- 
tion for all mixtures compared to that of single surfactants. 
Comparison of the distance between two molecules (L;) 
and available volume of surfactant molecule in the solution 
(V4) (Table 2) with the values of free volume (Vj), that is 
the average volume in which a single molecule can move 
inside a hypothetical shell due to surrounding molecules 
repulsion seemed worthwhile (Chauchan et al., 2018). The 
free volume also referred to as a void space between mole- 
cules caused by irregular molecular packing, and can be 
calculated using the formula (Chauchan et al. 2018; 
Chauhan and Sharma, 2015): 


Ze : 


where K’ = 4.28 x 10? is a constant independent of liquid 
nature. Figure 8 shows that V; values decrease with concen- 
tration at high concentration of surfactants, and the addition 
of more ELP mixture increases V;. A free volume increase 
results in greater liquid disorder due to increasing molecu- 
lar mobility (Agnihorti and Adgaonkar, 1989) and V; is 
dependent on the relative strength of contractive and expan- 
sive forces. Thus, У; increases with increasing ELP fraction 
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may be a result of dispersion forces, steric hindrance as 
well as unfavorable geometric fitting. Vy increase can also 
be caused by more disorder of water molecules around sur- 
factant molecules as well as decrease in oxyethylene chain 
hydration. Ultrasonic can be used to obtain the primary 
hydration number, np, from the expression (Burakowski 
and Glinski, 2007; Singh and Kumar, 2011): 


ње (1-5) (10) 


where nw and ns are the numbers of moles of water and 
solute, respectively, and ks, o is the pure water isentropic 
compressibility. m, is the number of water molecules near 
the solute the properties of which are changed compared to 
those of the bulk solvent. Water molecules in this shell can 
be considered to be incompressible due to being tightly 
trapped. Figure 9 shows that the calculated values of ny for 
mixtures at C more than 1077 M are higher than those for 
single surfactants indicating that mixtures have more struc- 
tured water molecules, which results in lower compressibil- 
ities. With surfactant aggregate formation, highly structured 
and low compressibility water molecules released from the 
vicinity of the hydrophobic part of surfactant become bulk 
water. The water molecules are highly structured around 
the hydrophobic part being of a rather low compressibility 
compared to the bulk water (Duboué-Dijon and Laage, n. 
d.; Desai and Dixit, 1996; El Eini et al., 1976). 

The viscosity B coefficient obtained by fitting the experi- 
mental viscosity data with the Jones—Dole equation based 
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Fig. 10 A plot of the values of (7, — 1)C~°° of the aqueous solutions of RH40, ELP and their binary mixtures at the mole fraction of ELP equal 
to 0.2, 0.4, 0.6, and 0.8 at T = 293 K vs. С05 as well as B-coefficients determined from the Jones—Dole equation vs. a 


on the (4, — 116779 plots against C?? where 4, is the rela- 
tive viscosity and C is the molar concentration is regarded 
as a measure of solute molecules size and shape as well as 
structural ramifications of the solute-solvent interactions 
(Ali et al., 2005; Jones and Dole, 1929). Figure 10 shows 
that B was linear with the mole fraction of ELP for the post 
micellar region. Moreover, Bvalues obtained based on this 
equation were positive consistent with the water—water 
structure-breaking character of surfactant molecules. The 
greater kosmotropic effect in aqueous solutions of surfac- 
tant mixtures as well as more intensive solute-solvent inter- 
actions for the mixtures and ELP solutions are indicated by 
the higher positive values of B for ELP and mixtures 
with RH40. 


Conclusions 


Stronger interactions between ELP and RH40 in the mix- 
ture than in the single-surfactant solutions as well as the 
synergetic effect are indicated by the calculated acoustical 
parameters, especially molecular free length and free vol- 
ume, and molecular interaction parameter. Moreover, more 
conversion to bulk water of the structured water molecules 
is found for the mixtures, which results in lower compress- 
ibilities and higher values of hydration numbers. The posi- 
tive values of the Jones—Dole coefficient B for the studied 
mixtures determined from the viscosity measurements are 
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also of significant importance as they indicate the water— 
water structure breaking character of surfactant molecules. 


Conflict of Interest The authors declare that they have no conflict 
of interest. 


References 


Agnihorti, P. K., & Adgaonkar, C. S. (1989) Theoretical evaluation of 
ultrasonic velocity in binary liquid mixtures. Ultrasonics, 27: 
248-251. 

Ali, A., Hyder, S., & Akhtar, Y. (2005) Viscometric studies of a- 
amino acid in aqueous NaCl and MgCl, at 303 K. Indian Journal 
of Physics, 79:157-160. 

Bakonyia, M., Berkóa, S., Kovácsa, A., Budai-Szńcsa, M., Kisa, N., 
Erös, G., ... Csányi, E. (2018) Application of quality by design 
principles in the development and evaluation of semisolid drug car- 
rier systems for the transdermal delivery of lidocaine. Journal of 
Drug Delivery Science and Technology, 44:136-145. 

Bali, V., Ali, M., & Ali, J. (2011) Nanocarrier for the enhanced bio- 
availability of a cardiovascular agent: In vitro, pharmacodynamic, 
pharmacokinetic and stability assessment. International Journal of 
Pharmaceutics, 403:46-56. 

Berthelsen, KR., Holm, R., Jacobsen, J., Kristensen, J., 
Abrahamsson, B., & Müllertz, A. (2015) Kolliphor surfactants 
affect solubilization and bioavailability of fenofibrate. Studies of 
in vitro digestion and absorption in rats. Molecular Pharmaceutics, 
12:1062-1071. 

Burakowski, A., & Gliński, J. (2007) Hydration numbers of non- 
electrolytes — Application of the acoustic method of Pasynski. 
Chemical Physics, 332:336-340. 

Chauchan, S., Sharma, V., Kaur, M., & Chaudhary, P. (2018) Tem- 
perature-dependent aggregation of bio-surfactant in aqueous 


J Surfact Deterg (2019) 


J Surfact Deterg 


solutions of galactose and lactose: Volumetric and viscometric 
approach. Chinese Journal of Chemical Engineering, 26: 
1119-1131. 

Chauhan, S., Chaudhary, P., Sharma, К., & Kumar, K. K. (2013a) 
Temperature-dependent volumetric and viscometric properties of 
amino acids in aqueous solutions of an antibiotic drug. Chemical 
Papers, 67:1442-1452. 

Chauhan, S., Kumar, K., & Patial, B. S. (2013b) Study of acoustical 
parameters f proline in lecithin-ethanol mixture at varying tempera- 
ture. Indian Journal of Pure and Applied Physics, 51:531—541. 

Chauhan, S., & Sharma, K. (2015) Extended studies on molecular 
interactions of SDBS and DTAB in aqueous solutions of amino 
acid at T = 293.15-313.15 K. Journal of Molecular Liquids, 211: 
675—685. 

Chen, Y., Li, G., Wu, X., Chen, Z., Hang, J., Qin, B., ... Wang, R. 
(2008) Self-microemulsifying drug delivery system (SMEDDS) of 
vinpocetine: Formulation development and in vivo assessment. Bio- 
logical & Pharmaceutical Bulletin, 31:118—125. 

Christiansen, A., Backensfeld, T., & Weitschies, W. (2010) Effects of 
nonionic surfactants on in vitro triglyceride digestion and their sus- 
ceptibility to digestion by pancreatic enzymes. European Journal 
of Pharmaceutical Sciences, 41:376—382. 

Constantinides, P. P., & Wasan, K. M. (2007) Lipid formulation strat- 
egies for enhancing intestinal transport and absorption of P- 
glycoprotein (P-gp) substrate drugs: In vitro/in vivo case studies. 
Journal of Pharmaceutical Sciences, 96:235—248. 

Cuiné, J. F., McEvoy, C. L., Charman, W. N., Рошоп, C. W., 
Edwards, G. A., Benameur, H., & Porter, C. J. (2008) Evaluation 
of the impact of surfactant digestion on the bioavailability of dana- 
Zol after oral administration of lipidic self-emulsifying formulations 
to dogs. Journal of Pharmaceutical Sciences, 97:995-1012. 

Desai, T. R., & Dixit, S. G. (1996) Interaction and viscous properties 
of aqueous solution of mixed cationic and nonionic surfactant. 
Journal of Colloid and Interface Science, 177:471—477. 

Duboué-Dijon, E., & Laage, D. (2015) Characterization of the local 
structure in liquid water by various order parameters. Journal of 
Physical Chemistry B, 119:8406—8418. 

El Eini, D. I., Barry, B. W., & Rhodes, J. (1976) Micellar size, shape, 
and hydration of long-chain polyoxyethylene nonionic surfactants. 
Journal of Colloid and Interface Science, 54:348—351. 

Gao, Р., Rush, B. D., Pfund, УУ. P., Huang, T., Bauer, J. M., 
Morozowich, W., ... Hageman, M. J. (2003) Development of a 
supersaturable SEDDS (S-SEDDS) formulation of paclitaxel with 
improved oral bioavailability. Journal of Pharmaceutical Sciences, 
92:2386-2398. 

Gelderblom, H., Verweij, J., Nooter, К., Sparreboom, А., & 
Cremophor, E. L. (2001) The drawbacks and advantages of vehicle 
selection for drug formulation. European Journal of Cancer, 37: 
1590-1598. 

Hua, X. Y., & Rosen, M. J. (1982) Calculation of the coefficient in 
the Gibbs equation for the adsorption of ionic surfactants from 
aqueous binary mixtures with non-ionic surfactants. Journal of Col- 
loid and Interface Science, 87:469-477. 

Ingram, T., Storm, S., Glembin, P., Bendt, S., Huber, D., 
Mehling, T., & Smirnova, I. (2012) Aqueous surfactant two-phase 
systems for the continuous countercurrent cloud point extraction. 
Chemie Ingenieur Technik, 84:840—848. 

Jones, G., & Dole, M. (1929) The viscosity of aqueous solutions of 
strong electrolytes with special reference to barium chloride. Jour- 
nal of the American Chemical Society, 51:2950—2964. 

Kiełek, K., ść Marczak, W. (2010) Hydration of non-electrolytes in 
H,O and D5O investigated by Passynski’s method. International 
Journal of Thermophysics, 31:77-84. 

Kumar, R., Mahesh, R., Shanmugapriyan, B., & Kannappan, V. 
(2012) Volumetric, viscometric, acoustic and refractometric 
studies of molecular interactions in certain binary systems of 


J Surfact Deterg (2019) 


o-chlorophenol at 303.15 K. Indian Journal of Pure and Applied 
Physics, 50:633—640. 

Lam, H., Кауооя, M., Haynes, C. A., Wang, D. І C, & 
Blankschtein, D. (2005) Affinity-enhanced protein partitioning in 
decyl beta-Dglucopyranoside two-phase aqueous micellar systems. 
Biotechnology and Bioengineering, 89:381-392. 

Rodriguez, A., del Mar, G. M., & Moya, M. L. (2008) Effects of addi- 
tion of polar organic solvents on micellization. Langmuir, 24: 
12785-12792. 

Rosen, J. M. (2004) Surfactants and interfacial phenomena (3rd ed.). 
New York, NY: Wiley Interscience. 

Rubingh, D. N. (1979) Mixed Micelle Solutions. In K. Mittal (Ed.), 
Solution chemistry of surfactants (pp. 337-354). New York, NY: 
Plenum Press. 

Sachin, K. M., Karpel, S. A., Singh, M., & Bhattarai, A. (2019) Self- 
assembly of sodium dodecylsulfate and dode- 
cyltrimethylammonium bromide mixed surfactants with dyes in 
aqueous mixtures. Royal Society Open Science, 6:181979. 

Sadurni, N., Solans, C., Azemar, N., & Garcia-Celma, M. J. 
(2005) Studies on the formation of O/W nano-emulsions, by 
low-energy emulsification methods, suitable for pharmaceutical 
applications. European Journal of Pharmaceutical Sciences, 26: 
438-445. 

Sannaningannavar, F. M., Navati, B. S., & Ayachit, N. H. (2013) 
Studies on thermo-acoustic parameters of poly(ethylene glycol)- 
400 at different temperatures. Journal of Thermal Analysis and 
Calorimetry, 112:1573-1578. 

Seelig, A., & Gerebtzoff, G. (2006) Enhancement of drug absorption 
by noncharged detergents through membrane and P-glycoprotein 
binding. Expert Opinion on Drug Metabolism & Toxicology, 2: 
733-752. 

Shaukat, A. (2011) Koliphor (TM) RH 40. In T. Reintjes (Ed.), Solu- 
bility enhancement with BASF pharma polymers—solubilizer com- 
pendium (pp. 89-94). Germany: Lamperthe, BASF SE Pharma 
Ingredients & Services. 

Singh, T., & Kumar, A. (2011) Thermodynamics of dilute aqueous 
solutions of imidazolium based ionic liquids. The Journal of Chem- 
ical Thermodynamics, 43:958-965. 

Sugihara, G., Miyazono, A., Nagadome, S., Oda, T., Hayashi, Y., & 
Ko, J. S. (2003) Adsorption and micelle formation of mixed surfac- 
tant systems in water II: A combination of cationic gemini-type sur- 
factant with MEGA-10. Journal of Oleo Science, 52:449-461. 

Suys, E. J. A., Warren, D. B., Pham, A. C., Nowell, С. J., 
Clulow, A. J., Benamuer, H., ... Chalmers, W. P. (2019) A non- 
ionic polyethylene oxide (PEO) surfactant model: Experimental 
and molecular dynamics studies of Kolliphor EL. Journal of Phar- 
maceutical Sciences, 108:193—204. 

Szymczyk, K. (2011) The properties of binary mixtures of ethoxylated 
octyl phenols with ethoxylated fluorinated alkanols at the water/air 
interface. Journal of Surfactants and Detergents, 14:415-423. 

Wang, S. W., Monagle, J., McNulty, C., Putnam, D., & Chen, H. 
(2004) Determination of P-glycoprotein inhibition by excipients 
and their combinations using an integrated high-throughput pro- 
cess. Journal of Pharmaceutical Sciences, 93:2755—2767. 

Wang, Z. H., Li, G. Z., Zhang, G. Y., Diao, Z. Y., Chen, L. S., & 
Wang, Z. W. (2005) Molecular interaction in binary surfactant mix- 
tures containing alkyl polyglycoside. Journal of Colloid and Inter- 
face Science, 290:598—602. 


Biographies 


Katarzyna Szymczyk studied Chemistry at Maria Curie Skłodowska 
University in Lublin, Poland from 1997 to 2002. Then after gradua- 
tion, she was employed in the Department of Interfacial Phenomena, 


WILEY AOCS% 


Faculty of Chemistry, UMCS in Lublin where she obtained her 
degrees (M.Sc. 2002, Ph.D. 2007, Assistant Professor 2015). She has 
published over 60 original articles in professional journals and books. 
Her research areas are interfacial phenomena at the solid—liquid, 
solid-gas, and liquid-liquid interfaces, adsorption, wettability, adhe- 
sion, and surface as well as volumetric properties of surfactants and 
their mixtures with the active substances of plant origin. 


Magdalena Szaniawska is a doctoral course student in the Depart- 
ment of Interfacial Phenomena, Maria Curie-Sklodowska University 


WILEY AOCS% 


J Surfact Deterg 


in Lublin (Poland). Her main research interests are the volumetric and 
adsorption properties of aqueous solution of nonionic surfactants and 
their mixtures with biologically active compounds at different 
temperatures. 


Konrad Terpiłowski was born in Poland in 1979. He studied chemis- 
try at Maria Curie-Skłodowska University in Lublin and graduated in 
2003. At present, he is an assistant professor in the Department of Inter- 
facial Phenomena UMCS Lublin. His research work concerns apparent 
surface free energy of solids and stability of dispersed systems. 


J Surfact Deterg (2019) 


ZALACZNIK 2 


PUBLIKACJA [D2] 


K. Szymczyk*, M. Szaniawska, J. Krawczyk, Temperature effect on the adsorption and 
volumetric properties of aqueous solutions of Kolliphor ELP. Molecules, 2020, 25(3), 
743. 


{ molecules (МӨР 


Article 


Temperature Effect on the Adsorption and Volumetric 
Properties of Aqueous Solutions of Kolliphor®ELP 


Katarzyna Szymczyk *©, Magdalena Szaniawska and Joanna Krawczyk 
Department of Interfacial Phenomena, Institute of Chemical Sciences, Faculty of Chemistry, 
Maria Curie-Skłodowska University in Lublin, Maria Curie-Skłodowska Sq. 3, 20-031 Lublin, Poland; 
magdalena.szaniawskaQpoczta.umcs.lublin.pl (M.S.); j.krawczykG&poczta.umcs.lublin.pl (J.K.) 
* Correspondence: katarzyna.szymczykQpoczta.umcs.lublin.pl; Tel.: +48-81-537-5538; Fax: +48-81-533-3348 


check for 
Received: 10 January 2020; Accepted: 6 February 2020; Published: 9 February 2020 © updates 


Abstract: Density, viscosity and surface tension of Kolliphor® ELP, the nonionic surfactant aqueous 
solutions were measured at temperature T = 293-318 K and at 5K interval. Steady-state fluorescence 
measurements have been also made using pyrene as a probe. On the basis of the obtained results, 
a number of thermodynamic, thermo-acoustic and anharmonic parameters of the studied surfactant 
have been evaluated and interpreted in terms of structural effects and solute-solvent interactions. 
The results suggest that the molecules of studied surfactant at concentrations higher than the critical 
micelle concentration act as structure makers of the water structure. 


Keywords: Kolliphor® ELP; volume expansivity; thermo-acoustic parameters; Gruneisen parameter; 
fractional free volume; coiling and micropolarity index 


1. Introduction 


Considering the landscape of current drug development, it can be stated, that40% of NCE (new 
chemical entities) are characterized by poor water solubility [1]. Therefore, there is a need for excipients 
to solubilize such candidates in both the early preclinical and clinical evaluation, as well as for the 
development of the marketed drug dosage forms. Solubilization is the process of drug uptake through 
complex formation into, e.g., oligomers of dextrose and fatty acids, through the cosolvent systems 
(such as ethanol, polyethyleneglycol and glycerol), or through the surfactant systems [2]. Contrary to the 
expectation that pharmaceutical excipients are pharmacologically inactive, there is abundant evidence 
that they can influence drug metabolism and efflux transport [3-5]. Some surfactants, particularly 
nonionic Kolliphors which were known earlier as Cremophors, find application in pharmaceutical 
formulations as the forms of solid dosage and delivery systems based on lipids aimed at improvement 
of poorly water-soluble drugs bioavailability [2,6]. The triricinoleate ester of ethoxylated glycerol 
is the main component of Kolliphors. The others are polyethylene glycol ricinolates, as well as 
the corresponding free glycols [7]. Kolliphor® EL (EL), being the most often used Kolliphor, was 
applied as a vehicle in the case of solubilization of some hydrophobic drugs, including cyclosporin 
A, diazepam, propofol and paclitaxel [8]. However, some researchers stated that EL is not an inert 
vehicle but exerts a range of biological effects, some of which have important clinical implications. Its 
use was associated with severe anaphylactoid hypersensitivity reactions, hyperlipidemia, abnormal 
lipoprotein patterns, aggregation of erythrocytes and peripheral neuropathy. The pharmacokinetic 
behavior of EL is dose-independent, although its clearance is highly influenced by the duration of 
the infusion [9,10]. This is particularly important since EL can affect the disposition of various drugs 
by changing the unbound drug concentration through micellar encapsulation. From this point of 
view, a clear understanding of the biological and pharmacological role of surfactants is essential to 
help oncologists avoid side-effects associated with the use of paclitaxel or other agents. On the other 
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hand to describe the behavior of surfactants in relation to two non-miscible phases and the range of 
temperatures at which they are active, some adsorption and volumetric properties of the solution 
of such surfactants are needed. Thus, the purpose of the presented studies was to determine some 
adsorption, volumetric, thermo-acoustic as well as anharmonic properties of aqueous solutions of 
Kolliphor® ELP (ELP), a purified grade of EL, by surface tension, density and viscosity measurements 
at T = 293-318 K with 5 K interval. Moreover, the properties of the solutions were studied by means 
of steady-state fluorescence measurements. Based on the results, the analysis was applied to study 
different molecular interactions in the solutions taking into account the change of these properties 
depending on concentration and temperature. 


2. Results and Discussion 


Table 1 presents the calculated different parameters determining the surface and bulk properties of 
aqueous solutions of ELP based on the values of their surface tension (y; y) (Figure 1). These parameters 
include: CMC (critical micelle concentration), surfactant efficiency for water surface tension reduction, 
which means the concentration of the surfactant for 20 mN/m reduction of the surface tension (C29), 
preference of the surfactant for adsorption in relation to the micelle formation with the pointer 
indicating possible reduction of water surface tension due to the surfactant presence (CMC/ C20), 
the surface pressure at the CMC (IIcyc), the surface excess concentration at the surface saturation (Tm), 
the minimum surface area per molecule (An), as well as the standard free enthalpy (AG? and AGO 4), 
enthalpy (AR and AH? 1) and entropy (AS? i and AS? q) of micellization and adsorption [11-13]. 
For studied ELP solutions, I'm and Am were calculated from the equations: 


Cdyw 1 dyw 1 dyLv 


Ty = – = 1 

Е RTdC RT dlnC 2.303RT dlogC 0) 
1 

An = ——— 2 

m NT, ( ) 


where C represents the concentration of surfactant, R is a gas constant, Т is temperature and N is an 
Avogadro number [11,12]. Next, the values of AGO. and AG? j were determined: 


CMC 


0 _ 
AG ic = RT In E (3) 
II 
AG. = AGinic — ( смс) (4) 
m 


where w is the number of water moles in 1 dm?. 
Knowing the values of AG? at different temperatures, it was possible to calculate AH? and AS”. 
If it is assumed that in a range of temperature, AH? is constant, then: 


p ZY (5) 


On the other hand, if AS? is constant, it is obtained [11]: 


= -AH? (6) 


Molecules 2020, 25, 743 


Table 1. Values of CMC, Сэр, CMC/ Coo, Псмс, AGO 


mic” 


AS? 


mic” AH io! 


the aqueous solutions of ELP at the temperatures from 293 to 318 K. 


3 of 13 


Lm, Am, AG? di AS? „and AHO for 


T = 293K 298K 303K 308K 313K 318K 
M 3 2.4x10? . 209x10? 203х105 1.97х 10-5 1.91 x 107? 1.85 x 1079 
(mol/dm?) 
C20 3 3.83 х 1076 3.82 х 1076 3.79 x 1076 3.74 x 1076 3.65 x 1076 3.53 x 1076 
(mol/dm?) 
CMC/ Co 5.59 5.47 5.35 5.26 5.23 5.24 
Псмс 32.91 32.82 32.72 32.59 32.43 32.29 
(mN/m) 
AG". 
(fc -26.19 -26.70 -27.22 -27.75 -28.28 -28.81 
А50. 
(kJ/mol K) 

(Јат) 4.58 4.59 4.59 4.59 4.59 4.58 
Tin (mol/m?)  317х10:5 3.11х10%6 3.06 х 1076 3.01 х 106 2.95 х 1076 2.90 х 1076 
Am (nm?) 0.524 0.533 0.542 0.552 0.563 0.572 

0 
CA —36.58 —37.23 -37.91 -38.58 -39.27 —39.93 
(kJ/mol) 
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d 
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Figure 1. A plot of the values of the surface tension (угу) of the aqueous solution of Kolliphor © ELP 
(ELP) at T = 293 K, 298 K, 303 K, 308 K, 313 K and 318 K vs. the logarithm of the surfactant concentration, 


log C. 


As follows from the table, the calculated values of AG? i and AG? qare negative, indicating that the 


ELP molecules have a trend to adsorb at the water-air interface and to form micelles in the bulk phase, 
as well as the two processes are spontaneous. The adsorption free enthalpy values are more negative 
in comparison with the micellization values indicating greater propensity of ELP for adsorption at the 
interface than micelles formation in the bulk phase. The AG? and AG? , Values become more negative 
with the increasing temperature probable due to the greater stability of the adsorbed and micellized 
molecules of ELP compared to the freely dispersed in the aqueous phase. However, the positive values 
of ASQ ds show that, after adsorption and micellization, the studied solutions become more random [11]. 
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On the other hand, АН. and AH? 4, positive values the indicate predominance of bond breaking 
when micellization and adsorption proceed. Figure 1 shows the gradual decrease of ELP surface 
tension values with the increasing temperature being 293-318 K, which is similar to the I „I lIcyc and 
CMC/C values (Table 1). As follows from the data presented in Table 1, the ELP CMC values drop 
slightly with the temperature rise from 293 to 318 K. This indicates that temperature increase can break 
down the intra-hydrogen bonds between the surfactant molecules and weaken the hydration action of 
the hydrophilic groups, which favors the micelle formation [14]. 

The surface tension isotherm (угу = f (log C)) is known to not be the only way to determine 
CMC. What is more, the values of CMC for a surface-active agent very often differ depending on 
its determination method [15]. Considering the values of the dynamic viscosity (7) of ELP aqueous 
solutions (Figure 2 as an example), it is evident that the great increase of the л values take place at 
293 K and the concentration exceeding 1077 M. 


7 (mPa.s) 
п (mPa.s) 


T(K) 


Figure 2. A plot of the values of 7 of aqueous solutions of ELP at C from 107* to 107? M vs. the 
temperature, T, as well as the values of 7 of the aqueous solutions of ELP at T = 293 K vs. log C. 


Based on the relations between the relative viscosity (7/0) and concentration C [16] of ELP at 
a given T (Figure 3 as an example), the CMC values are as follows: 9.87 х 1076 M, 9.84 x 1076 M, 
9.73 x 1075 M, 9.66 x 107 M, 9.52 x 1076 M and 9.32 x 1075 M in the temperature range 293-318 K. 
They are similar to the values given in Table 1. However, besides the determination of CMC from the 
viscosity measurements, it was attempted to relate the surface tension with viscosity [17]. Moreover, 
the basic equation proposed by Pelofsky is the linear relation and can be used for organic and inorganic 
phases of pure and mixed components [18]. It is interesting that, taking into account the measured 
values of yry and r at a given temperature, and C corresponding to the saturated monolayer at the 
water-air interface, there is the linear dependence proposed by Pelofsky that is between In y; y and 
1/1 for the ELP solutions (Figure 3). 

On the other hand, the constant B in the relation presented by the author as a function of the 
molecular weight and thermal conductivity possesses a positive value at each temperature contrary to 
the values for n-alkanes, n-alcohols, water as well as some aqueous solutions [17]. In our opinion, the B 
parameter values can be better explained by conducting additional measurements for a larger amount 
of surfactant concentrations. Therefore, it is certain that the measured dynamic viscosity (7) of ELP 
solutions is greatly sensitive to changes in temperature (exemplary Figure 2) following the Arrhenius 
law [19,20]: 


Eq 
л = Bexp эт (7) 
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where B is the is pre-exponential factor and Е, is shear activation energy. The calculated from 
Equation (7) values of E, (Figure 4), that is energy necessary for individual micelles motion in an 
environment of surrounding micelles determined on the basis of this law, increases significantly at ELP 
concentrations in the bulk phase higher than 1077 M. Moreover, the highest value of E, being 16.35 
kJ/mol is found at C = 107? M. Furthermore, the enthalpy of activation (AH") values, as well as the 
change in heat capacity of activation (AC;), are affected to a great extent by temperature (exemplary 
Figures 5 and 6). The values of AH* and AC, were determined based on the relations proposed by 


Mukherjee et al. [21,22]: 


АН" dinn 8) 
RT? dT 
АС" ААН" 0) 
P dT 
C 
0.00001 0.000020 0.00008 000004 0.00005 
SS Se 1.0116 
415 F- T= 203K 
1.0114 
4.10 
dE 1.0112 
5 
Е 
> = 
ё 400 1.0110 
3.95 
1.0108 
3.90 
1.0106 
3.85 А 1 a 1 А 1 А 1 : 1 


0.993890 0.993895 0.993900 0.993905 0.993910 0.993915 
1/n 


Figure 3. A plot of the values of the In yzy vs. 1/7 as well as 1/0 vs. C for the aqueous solutions of 
ELP at T = 293 K. 
T(K) 
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Figure 4. A plot of the values of shear activation energy (E,) of the aqueous solutions of ELP vs. log 
C, as well as the values of the partial molar volume (Vy) of the aqueous solutions of ELP vs. the 


temperature, T. 
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Figure 5. A plot of the values of the enthalpy of activation (AH*) of the aqueous solutions of ELP at C 
from 1074 to 107? M vs. the temperature, T, as well as the values of AH* of the aqueous solutions of 
ELP at T = 293 K vs. log C. 
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Figure 6. A plot of the values of the change in heat capacity of activation (AC) of the aqueous solutions 
of ELP at C from 107* to 1072 M vs. the temperature, T, as well as the values of АС, of the aqueous 
solutions of ELP at T = 293 K vs. log C. 


What is more, the АН” values calculated from Equation (8) for the ELP solutions are positive and 
expectedly declined with increasing temperature (Figure 5) indicating that the processes are connected 
with the heat absorption in the solutions under considerations. From this reason, the АС, values are 
negative and decrease linearly with T (Figure 6). 

As the ELP molecules are strongly hydrated by water, and it is possible that the H3O* ions 
can be associated to the oxyethylene groups, similarly to other nonionic surfactants studied by 
us [15,22-25], it was also interesting to analyze the viscosity data according to the Jones-Dole equation 
and A and B coefficients [26,27]. This equation and the relationship between (n, — 1)C~°° and C9? 
(Figure 7) for C equal to and higher than 5 x 1073 M, where nr is the relative viscosity, were used 
for the determination of B coefficients values which proved to be positive (Figure 8), disclosing the 
water-structure-breaking nature of the ELP molecules. It is interesting that the value of dB/dT, which 
by some authors is designed as a better criterion for determining any solute effect on the structure of 
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solutions depending on temperature [28], for the studied concentrations of ELP solutions, is equal 
to —0.21, and according to Hugue et al. [28], indicates that the solute is a structure maker. To solve 
an idea about the structure-making or breaking role of the ELP in the solution, the apparent molar 
volume values, py [29,30], were determined from the density measurements (p) (exemplary Figures 9 
and 10) and the equation which has the form: 


1000(po — 
(= M | | p) (10) 


where M is the molecular weight of the surface-active agent and ро is the density of the “pure” solvent. 
The calculated values of py were then analyzed based on the Hepler equation [31]. It proved that 
at C equal to and greater than 5 x 1077 M, the values of (9200 / ar), where p? is the apparent 
molar volume at infinite dilution, are positive, indicating the structure-making properties of the 
investigated surfactant. 
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Figure 7. A plot of the values of (n, — He? of the aqueous solutions of ELP at T = 293 K, 298 K, 303 K, 
308 K, 313 K and 318 K vs. СО 
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Figure 8. A plot of the values of A and B coefficients determined from the Jones—Dole equation and 
viscosity of the aqueous solutions of ELP vs. T. 
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Figure 9. A plot of the values of p of the aqueous solutions of ELP at C from 1074 to 1072 M vs. the 
temperature, T, as well as the values of p of the aqueous solutions of ELP at T = 293 K vs. log C. 
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Figure 10. A plot of the values of the apparent molar volume, py, of the aqueous solutions of ELP at C 
from 10735 to 1072 M vs. the temperature, T. 


However, the changes of the py values with T before and after C = 107^ M (shown in Figure 10) are 
different pointing out the great structural changes of ELP molecules due to the increase in temperature 
and/or those of intermolecular interactions between the surfactant and water. They, in turn, can be 
determined, among others, based on the volume expansivity values a (Figure 11), which are also 
regarded as the thermal expansion coefficient being a measure of volume change with the temperature 
and can be calculated based on the measurements of density from the equation [18,28]: 


a= ze] (11) 


where Vn is the partial molar volume. It follows from Figure 11 that for ELP, the calculated values of a 
increase with the rise of T and C, which is consistent with the increasing tendency of Vy for studied 
solutions (Figure 4) at the concentrations higher than CMC of the studied surfactant determined from 


the surface tension measurements (Table 1). 
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Figure 11. A plot of the volume expansivity a of the aqueous solutions of ELP at C from 1074 to 1072 M 
vs. the temperature, T, as well as the values of a of the aqueous solutions of ELP at T = 293 K vs. log C. 


Having the ELP a values at various temperatures and concentrations, there could be calculated the 
thermodynamic parameters such as: reduced volume (V), Moelwyn-Hughes parameter (C1), reduced 
compressibility (B), isochoric temperature coefficient of internal pressure (X), Sharma parameter (50), 
Huggin's parameter (F), isochoric temperature coefficient of volume expansivity (X’), anharmonic 
microscopic isothermal Gruneisen parameter (I), fractional free volume (f), Gruneinsen parameter 
(L5), isobaric thermo-acoustic parameter (К) and isochoric thermo-acoustic parameter (K") [32-35]. 
These parameters for ELP solutions are presented in Table $1 (Supplementary Materials) and were 
calculated from the following equations: 
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Table 51 shows the variable character of бо and T which differs from the suggestions given 
by Sharma et al. [20,36]. They stated that Sp is the constant in any liquid or solid-state system. 
However, fractional free volume (f) values, which are expressed in terms of the repulsive exponent of 
intermolecular potential, are characterized by the non-linear rise with T pointing out that the surfactant 
molecules mobility results in greater liquid (surfactant + water) disorder, which is due to irregular 
molecular packing [35]. Simultaneously the C1, X”, F, Г and Г, the decreasing values with the increasing 
temperature and concentration of the surfactant indicating the molecular ordering increase with the 
increasing values of T and C [37]. This is consistent with the conclusions drawn based on the viscosity 
measurements namely that the ELP molecules act as a structure maker at its high concentrations. 
One should bear in mind that ELP, as a nonionic surfactant possessing three amphiphilic chains, is a 
nonconventional surfactant that can be treated as a polymer. The coiling index was determined based 
on the pyrene fluorescence spectroscopy in order to reveal the changes of the chains' conformation. 
It is the ratio of the excimer intensity Ir (at ~480 nm) and the monomer, ly = (I4 + 13)/2, where П 
and Із are the intensity of the first and the third vibrionic peaks in the pyrene emission spectra [38]. 
The high value of Ig/Im indicates the medium with a greater possibility of excimer formation and 
coiled surfactant chains. As can be seen in Figure 12, the highest coiling index values are found at 
the concentrations close to the ELP CMC values calculated from the surface tension measurements as 
well as the C values with the polarity index (I; /I3) characterized by a minimal value. Considering the 
concentration values higher than CMC, there is observed the Ig / I; values’ decrease with the increasing 
number of micelles providing more pyrene locations or micellar size or shape changes. At C higher 
than 107? M,I E / Iu is still constant with the lowest values indicating stretching of the surfactant chains, 
which is consistent with the f values which increase remarkably at C greater than 1073 M at a given 
temperature (Table S1). 
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Figure 12. Plot of the values of Ig /Ім and I4 /I3 determined from the fluorescence spectra of pyrene in 
the aqueous solutions of ELP vs. log C. 


3. Materials and Methods 


Kolliphor? ELP (ELP) (Cremophor®ELP, Polyoxyl 35 Hydrogenated Castor Oil, Polyoxyl-35 
Castor Oil) purchased from Sigma were used without further purification. Its aqueous solutions 
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were obtained at the concentrations 1079 to 107? M applying the doubly distilled and deionized 
water provided by a Destamat Bi18E distiller. The surface tension was measured by means of a 
Kriiss K100C tensiometer according to the platinum ring tensiometer method (du Nouy’s method). 
The measurements of pure water surface tension at 293 K were aimed at the calibration of the 
tensiometer and glassware cleanliness control. There were performed ten successive measurements 
with the standard deviation not exceeding +0.2 mN/m. The controlled temperature was within +0.1 K. 
The measurements of surfactant aqueous solutions density were performed by means of a U-tube 
densitometer (DMA 5000 Anton Paar) at 293-318 K. The density and temperature measurement 
precision given by the producer was +0.000005 g cm Запа +0.001 К. The calculated uncertainty was 
0.01%. The densitometer calibration with distilled and deionized water was regular. 

The surfactant aqueous solutions viscosity was determined using the Anton Paar viscometer 
(AMVn) at 293-318 K +0.01 K. Its precision was 0.0001 mPa s and the uncertainty was 0.3%. 
The steady-state fluorescence was measured at 293-318 K by means of a Hitachi F-2700 Fluorescence 
Spectrometer where pyrene was as a luminescence probe (C Py = 4 х 1079 M). The fluorescence 
excitation for pyrene was induced at 335 nm and the range of the emission spectra recording was 
350-600 nm at a scan speed of 300 nm/min. The widths of the excitation and emission slit widths were 
5 nm. 


4. Conclusions 


In this paper, some physicochemical properties of the aqueous solutions of the nonionic surfactant, 
Kolliphor®ELP (ELP), are singled out, discussed and compared on the basis of the surface tension, 
density and viscosity measurements as well as fluorescence spectra. From the presented data and 
calculations, it is evident that ELP displays a greater propensity to be absorbed at the interface than 
form micelle in the bulk phase and that bond breaking predominates in the micellization and adsorption 
process. The ELP structure-making tendency at concentrations higher than 107? M was confirmed 
by the Hepler's theory as well as the values of viscosity B coefficients and apparent molar volume. 
In addition, at this concentration range and a given temperature, the stretching of the surfactant 
chains were proved on the basis of pyrene emission spectra and the values of fractional free volume. 
On the other hand, at a concentration of ELP smaller than CMC, the polar head of the surfactant is 
strongly hydrated resulting in very compact conformation, probably because of water-bridging, which 
promotes gauche conformations. 


Supplementary Materials: The following are available online, Table S1: Values of V, C, В, X, So, F, X', V, f, Tp, К 
and K” for the aqueous solutions of ELP at the temperatures from 293 to 318 K. 
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Table S1. Values of V, Co Ë ,Х,5,.Р, X,L, f TK and K for the aqueous solutions of 
ELP at the temperatures from 293 to 318 K. 


C 293K 298K 303K 308K 313K 318K 
10? 1.0681 1.0848 1.1026 1.1215 1.1416 1.1629 
2x10? 1.0681 1.0848 1.1026 1.1215 1.1416 1.1629 
5x10? 1.0681 1.0848 1.1026 1.1215 1.1416 1.1630 
8x10? 1.0681 1.0848 1.1026 1.1215 1.1416 1.1630 
104 1.0682 1.0848 1.1027 1.1215 1.1417 1.1630 
2x104 1.0683 1.0850 1.1028 1.1216 1.1418 1.1631 
V  5x104 1.0685 1.0852 1.1030 1.1218 1.1419 1.1632 
8x104 1.0687 1.0854 1.1032 1.1221 1.1422 1.1634 
103 1.0690 1.0856 1.1034 1.1223 1.1423 1.1635 
2х107 1.0698 1.0865 1.1042 1.1230 1.1430 1.1642 
5x103 1.0723 1.0889 1.1067 1.1255 1.1455 1.1667 
8x103 1.0745 1.0912 1.1089 1.1278 1.1478 1.1690 
102 1.0764 1.0929 1.1105 1.1291 1.1488 1.1698 


10° 20.3750 17.4834 15.4432 13.9331 12.7715 11.8545 

2x10? 20.3730 17.4822 15.4421 13.9329 12.7711 11.8543 
5x10? 20.3710 17.4809 15.4406 13.9327 12.7706 11.8538 
8x10? 20.3657 17.4792 15.4399 13.9325 12.7702 11.8534 

104 20.3583 17.4786 15.4364 13.9321 12.7688 11.8529 

2x104 20.3303 17.4569 15.4230 13.9275 12.7625 11.8501 
C. 5x10^ 20.2835 17.4284 15.4094 13.9107 12.7575 11.8454 
8x10^ 202410 17.4045 15.3876 13.8912 12.7446 11.8369 

103 20.1879 17.3682 15.3695 13.8833 12.7383 11.8318 

2x10? 20.0034 17.2532 15.2946 13.8331 12.7040 11.8083 
5x10? 19.5278 16.9378 15.0743 13.6735 12.5852 11.7178 
8x10? 19.1156 16.6615 14.8801 13.5320 12.4793 11.6371 

102 18.7756 16.4544 14.7511 13.4517 12.4308 11.6098 

10° 3.8265 4.1502 4.5198 4.9423 5.4277 5.9861 

2x10? 3.8267 4.1504 4.5200 4.9423 5.4279 5.9862 
5x10? 3.8268 4.1505 4.5204 4.9424 5.4282 5.9866 
8x10? 3.8273 4.1508 4.5205 4.9425 5.4284 5.9869 

104 3.8279 4.1509 4.5213 4.9426 5.4291 5.9873 

2x10^ 3.8304 4.1540 4.5244 4.9442 5.4323 5.9893 

B 5x10^ 3.8345 4.1582 4.52775 4.9500 5.4349 5.9928 
8x104 — 3.8383 4.1617 4.5325 4.9568 5.4414 5.9990 
10? 3.8430 4.1670 4.5367 4.9595 5.4446 6.0028 
2x10? 3.8597 4.1842 4.5542 4.9773 5.4624 6.0202 
5x10? 3.9052 4.2335 4.6077 5.0353 5.5255 6.0891 
8x10? 3.9477 4.2794 4.6574 5.0893 5.5842 6.1529 
10? 3.9850 4.3157 4.6919 5.1211 5.6120 6.1750 
10° -0.5882 -0.5558 -0.5230 -0.4902 -0.4575 -0.4251 
2x10? -0.5882 -0.5558 -0.5230 -0.4902 -0.4575 -0.4251 


5x105 -0.5881 -0.5557 -0.5230 -0.4902 -0.4575 -0.4251 
8x105 -0.5881 -0.5557 -0.5230 -0.4902 -0.4575 -0.4251 
104 -0.5880 -0.5557 -0.5229 -0.4902 -0.4574 -0.4251 
2x104 -0.5878 -0.5554 -0.5227 -0.4901 -0.4572 -0.4250 
5x104 -0.5873 -0.5550 -0.5224 -0.4897 -0.4571 -0.4248 
8x104 -0.5869 -0.5547 -0.5220 -0.4892 -0.4567 -0.4245 
103 -0.5864 -0.5542 -0.5216 -0.4890 -0.4565 -0.4243 
2x10? -0.5847 -0.5526 -0.5202 -0.4877 -0.4554 -0.4233 
5x10? -0.5799 -0.5480 -0.5158 -0.4836 -0.4515 -0.4197 
8x10? -0.5756 -0.5439 -0.5119 -0.4798 -0.4479 -0.4164 
102 -0.5718 -0.5406 -0.5091 -0.4776 -0.4463 -0.4153 
10? 0.955981 0.918851 0.879754 0.839011 0.796822 0.753574 
2x10? 0.95596 0.918832 0.879727 0.839005 0.796807 0.753567 
5x10? 0.95594 0.918811 0.879693 0.838999 0.796784 0.75354 
8x10? 0.955884 0.918784 0.879679 0.838993 0.796768 0.753519 
10^ 0.955808 0.918773 0.879597 0.838979 0.796708 0.753491 
2x10^ 0.955518 0.918427 0.879288 0.838834 0.796447 0.753344 
So 5x104 0.95503 0.91797 0.878974 0.8383 0.796235 0.753088 
8x10^ 0.954583 0.917585 0.878467 0.83768 0.795695 0.752637 
10? 0.954023 0.916997 0.878046 0.837429 0.795432 0.752362 
2x10? 0.952041 0.91511 0.876283 0.835813 0.793983 0.751097 
5x10? 0.946694 0.909743 0.870947 0.830561 0.788862 0.746152 
8x10? 0.941754 0.904795 0.866039 0.825739 0.784168 0.74163 
10? 0.937447 0.900923 0.862671 0.822931 0.78197 0.740074 
10° 1.5046 1.4536 1.4017 1.3493 1.2964 1.2435 
2x10? 1.5046 1.4536 1.4017 1.3493 1.2964 1.2435 
5x10? 1.5046 1.4535 1.4016 1.3493 1.2964 1.2434 
8x10? 1.5045 1.4535 1.4016 1.3493 1.2964 1.2434 
104 1.5044 1.4535 1.4015 1.3492 1.2963 1.2434 
2x104 1.5040 1.4530 1.4011 1.3490 1.2960 1.2432 
F 5х10* 1.5033 1.4524 1.4007 1.3484 1.2957 1.2429 
8x104 1.5027 1.4519 1.4000 1.3476 1.2950 1.2423 
103 1.5019 1.4511 1.3995 1.3473 1.2947 1.2420 
2х107 1.4991 1.4485 1.3972 1.3452 1.2929 1.2405 
5x103 1.4917 1.4413 1.3902 1.3386 1.2866 1.2345 
8x10? 1.4848 1.4347 1.3839 1.3325 1.2808 1.2290 
10? 1.4789 1.4296 1.3795 1.3290 1.2781 1.2272 
10° -1.1253 -1.1533 -1.1820 -1.2114 -1.2416 -1.2725 
2x10? -1.1253 -1.1533 -1.1820 -1.2114 -1.2416 -1.2725 
5x105 -1.1254 -1.1533 -1.1821 -1.2115 -1.2417 -1.2725 
8x105 -1.1254 -1.1533 -1.1821 -1.2115 -1.2417 -1.2725 
104 -1.1255 -1.1533 -1.1821 -1.2115 -1.2417 -1.2726 
2x104 -1.1257 -1.1536 -1.1823 -1.2116 -1.2419 -1.2727 
X 5x104 -1.1261 -1.1539 -1.1826 -1.2120 -1.2420 -1.2728 
8x104 -1.1264 -1.1542 -1.1829 -1.2124 -1.2424 -1.2732 
103 -1.1268 -1.1547 -1.1833 -1.2126 -1.2426 -1.2734 
2x10? -1.1283 -1.1561 -1.1845 -1.2137 -1.2437 -1.2743 
5x10? -1.1324 -1.1600 -1.1884 -1.2175 -1.2473 -1.2778 
8x10% -1.1361 -1.1637 -1.1920 -1.2210 -1.2507 -1.2810 
102 -1.1394 -1.1665 -1.1944 -1.2230 -1.2522 -1.2822 
105 5.9944 5.6159 5.3479 5.1480 4.9923 4.8671 
2x105 5.9941 5.6157 5.3477 5.1480 4.9923 4.8670 
5x10? 59939 5.6155 5.3475 5.1480 4.9922 4.8670 
8x10? — 5.9932 5.6153 5.3475 5.1480 4.9921 4.8669 
104 5.9922 5.6152 5.3470 5.1479 4.9919 4.8668 
2x104 5.9886 5.6124 5.3452 5.1473 4.9911 4.8665 
I 5x104 59824 5.6087 5.3434 5.1450 4.9904 4.8658 
8x10^ — 5.9769 5.6055 5.3406 5.1424 4.9887 4.8646 
10? 5.9699 5.6008 5.3382 5.1414 4.9878 4.8639 
2x10? 5.9458 5.5857 5.3283 5.1347 4.9832 4.8607 
5x10? 5.8836 5.5443 5.2992 5.1134 4.9671 4.8481 
8x10? — 5.8296 5.5081 5.2736 5.0945 4.9527 4.8369 
102 5.7851 5.4809 5.2565 5.0838 4.9461 4.8331 


Ё 105 0.142971 0.151152 0.157533 0.162654 0.166881 0.170443 


2x10? 0.142977 0.151155 0.157536 0.162654 0.166882 0.170444 
5x10? 0.142982 0.151159 0.157541 0.162655 0.166884 0.170446 
8x10? 0.142996 0.151164 0.157543 0.162656 0.166886 0.170447 
10^ 0.143016 0.151166 0.157555 0.162657 0.166891 0.17045 
2x10^ 0.143091 0.151231 0.157599 0.162673 0.166914 0.170461 
5x10^ 0.143217 0.151317 0.157643 0.162733 0.166934 0.17048 
8x10^ 0.143331 0.151388 0.157715 0.162802 0.166982 0.170514 
10? 0.143474 0.151498 0.157774 0.162829 0.167006 0.170534 
2x10? 0.143972 0.151844 0.15802 0.163007 0.167135 0.170629 
5x10? 0.145274 0.152804 0.158749 0.163574 0.167587 0.170995 
8x10? 0.146421 0.153655 0.159398 0.164082 0.167991 0.171323 
10? 0.147382 0.1543 0.159833 0.164371 0.168178 0.171435 
105 10.0209 8.5751 7.5549 6.7999 6.2191 5.7606 
2х105 10.0198 8.5744 7.5544 6.7998 6.2189 5.7605 
5x10? 10.0188 8.5738 7.5536 6.7997 6.2186 5.7603 
8x10? 10.0162 8.5729 7.5533 6.7996 6.2184 5.7600 
104 10.0125 8.5726 7.5515 6.7994 6.2177 5.7598 
2x104 9.9985 8.5618 7.5448 6.7971 6.2146 5.7584 
5x104 — 9.9751 8.5475 7.5381 6.7887 6.2121 5.7560 
8x104 — 9.9538 8.5356 7.5271 6.7789 6.2056 5.7518 
103 9.9273 8.5174 7.5181 6.7750 6.2025 5.7492 
2x10? 9.8350 8.4599 7.4806 6.7499 6.1853 3.7375 
5x10? 9.5973 8.3022 7.3705 6.6701 6.1259 5.6923 
8x10? — 93911 8.1641 7.2734 6.5993 6.0730 5.6519 
102 9.2211 8.0605 7.2089 6.5592 6.0487 5.6382 
10° 9.6875 8.2417 7.2216 6.4666 5.8858 5.4272 
2x10? 9.6865 8.2411 7.2210 6.4665 5.8856 5.4272 
5x10? 9.6855 8.2405 7.2203 6.4664 5.8853 5.4260 
8x10? 9.6828 8.2396 7.2200 6.4663 5.8851 5.4267 
104 9.6791 8.2393 7.2182 6.4660 5.8844 5.4265 
2x104 9.6652 8.2285 7.2115 6.4638 5.8813 5.4251 
5x104 9.6418 8.2142 7.2047 6.4553 5.8787 5.4227 
8x104 9.6205 8.2023 7.1938 6.4456 5.8723 5.4185 
10? 9.5940 8.1841 7.1848 6.4417 5.8692 5.4159 
2x10? 9.5017 8.1266 7.1473 6.4165 5.8520 5.4041 
5x10? 9.2639 7.9689 7.0372 6.3367 5.7926 5.3589 
8x10? 9.0578 7.8308 6.9400 6.2660 5.7397 5.3186 
10? 8.8878 7.7272 6.8756 6.2259 5.7154 5.3049 
10° -3.6931 -2.6259 -1.8737 -1.3185 -0.8934 -0.5602 
2x105 -3.6924 -2.6254 -1.8733 -1.3185 -0.8933 -0.5601 
5x105 -3.6916 -2.6249 -1.8727 -1.3184 -0.8931 -0.5600 
8x105  -3.6897 -2.6243 -1.8725 -1.3183 -0.8930 -0.5598 
104 -3.6869 -2.6240 -1.8712 -1.3182 -0.8924 -0.5596 
2x104 -3.6766 -2.6161 -1.8663 -1.3165 -0.8902 -0.5586 
5x104 -3.6593 -2.6055 -1.8613 -1.3103 -0.8883 -0.5569 
8x104 -3.6436 -2.5967 -1.8532 -1.3031 -0.8836 -0.5538 
103 -3.6240 -2.5833 -1.8466 -1.3003 -0.8813 -0.5520 
2x10? -3.5559 -2.5409 -1.8190 -1.2818 -0.8688 -0.5435 
5x10? -3.3804 -2.4246 -1.7379 -1.2233 -0.8255 -0.5108 
8x10? -3.2282 -2.3227 -1.6664 -1.1715 -0.7870 -0.4817 
102 -3.1027 -2.2463 -1.6190 -1.1421 -0.7693 -0.4718 
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Abstract: The studies on the adsorption properties and composition of the adsorbed monolayer at the 
water-air interface of the binary Kolliphor® ELP (ELP) and Kolliphor® RH 40 (КН40) mixtures based 
on the measurements of the surface tension (yzy) of their aqueous solution in the temperature range 
from 293 to 318 К were carried out. The угу isotherms were described by the exponential function 
of the second order and the Szyszkowski equation as well as predicted by Fainerman and Miller 
equation. The obtained yzy isotherms were analyzed using the exponential function of the second 
order, the Szyszkowski, Fainerman and Miller as well as independent adsorption equations. The 
'YLy isotherms were also used for determination of the Gibbs surface excess concentration of RH40, 
ELP and their mixture (T) at the water-air interface as well as the mixed monolayer composition. 
Based on T and the constant a in the Szyszkowski equation, the standard thermodynamic functions of 
adsorption were considered. From the consideration dealing with the yzy isotherms obtained by us, 
it results, among others, that these isotherms for the non-ideal solution of macromolecular surfactants 
mixture can be predicted using the Fainerman and Miller equation. From this consideration, it also 
results that a simple method proposed by us, based on the isotherms of RH40 and ELP, allows us to 
predict the composition of their mixed monolayer in the whole concentration range of RH40 and ELP 
in the bulk phase. 


Keywords: Kolliphor? ELP; Kolliphor® RH 40; adsorption; thermodynamic parameters of adsorption; 
monolayer composition 


1. Introduction 


The wide application of surfactants in many industrial areas, agriculture, pharmacy, 
medicine and everyday life results from their ability to adsorb at the solid-liquid, solid-air 
and liquid-air interfaces and to form micelles in the bulk phase of the solutions [1—4]. 
Indeed, in practice the mixtures of the surfactants, rather than the individual ones, are 
used [4-7]. This is due to the specific properties of the mixed adsorption layers at the 
different interfaces or the resulting mixed micelles that cannot be predicted on the basis of 
the single surfactant properties. 

The properties of the adsorption layers at the interfaces depend on the kind of the 
surfactants, orientation of their molecules towards the interface, layers thickness and 
pacing [4]. In some practical applications, the flexibility of the adsorption layers at the water- 
air or water-oil interfaces plays an important role. It can be expected that such conditions 
can be satisfied by the surfactant mixtures, the molecules of which are large and branched. 
Such surfactants include nonionic Kolliphors (known before as Cremophors) which are 
represented by Kolliphor? ELP (ELP) and Kolliphor? RH 40 (RH40). These surfactants 
are applied, among others, in pharmacy [8,9]. Their main constituent is triricinoleate ester 
of ethoxylated glycerol (Scheme S1 in Supplementary Materials). Besides this constituent 
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the others are polyethylene glycol ricinoleate and the corresponding free glycols [10,11]. 
Three hydrocarbon chains with one -OH group each constitute the hydrophobic part of the 
RH40 and ELP molecules. One difference between the tail of RH40 and ELP is that in the 
hydrophobic ELP tail the double bond exists between the carbon atoms. The hydrophilic 
character of the molecule head results mainly from the presence of oxyethylene groups. 
Thus, at the first approximation the properties of the hydrophilic part of RH40 and ELP 
can be compared to those of Triton’s [4,12,13]. However, the RH40 and ELB molecules 
have a more complicated structure than Triton’s, and for this reason more cases of these 
molecule orientations are possible at the water-air interface [10,11]. These orientations can 
influence on the reduction in the water surface tension by the adsorbed RH40 and/or ELP 
molecules in the surface region. This may be the reason that the literature does not report 
the unambiguous opinion dealing with the influence of RH40 and ELP on the water surface 
tension. It is practically impossible to find description or prediction of the isotherm of the 
surface tension of the aqueous solution of these surfactants at different temperatures. 

As it is known, in practice the mixtures of surfactants are used, but not single surfac- 
tants. Therefore, it is important to predict the surface tension isotherms, composition and 
packing of the mixed monolayer at the water-air interface as well as the thermodynamic 
parameters of the adsorption of the surfactants mixture at this interface based on the surface 
tension isotherms of the aqueous solution of the single surfactants being components of the 
mixture. This issue was the main objective of our studies. These studies were based on the 
surface tension measurements of the aqueous solutions of ELP and КН40 as well as their 
mixtures in the temperature range from 293 to 318 K. The obtained isotherms of surface 
tension were considered in terms of the concentration of RH40 and ELP in the mixed 
monolayer at the water-air interface and the composition of this monolayer. Different 
concepts were applied in this consideration. To explain the adsorption properties of the 
RH40 and ELP mixtures at the water-air interface, its thermodynamic parameters were 
determined and analyzed. 


2. Results and Discussion 
2.1. Surface Tension 


The surface tension of the aqueous solution of RH40 and ELP depends on that of water, 
RH40 and ELP. According to van Oss et al. [14-17] the surface tension of liquids and solids 
can be divided into two components. One, which is present in the surface tension of each 
liquid or solid, results from the Lifshitz-van der Waals intermolecular interactions (y DW) 
and the other one results from the acid-base intermolecular interactions (4/4). The acid- 
base component depends on the electron-acceptor (y™) and electron-donor (y~) parameters. 
However, van Oss et al. [14-17] suggest that the contribution of the dipole-dipole and 
induced dipole-dipole interactions to the Lifshitz-van der Waals component of the surface 
tension is smaller than 2%. This indicates that the у“ component of the solid and liquid 
surface tension depends only on the dispersion intermolecular interactions. 

Fowkes [18] stated that the dispersion intermolecular interactions at the interface can 
be deduced on the basis of the interactions of one element with twelve neighbors and the 
distance between particular elements. In the case of the aqueous solution of RH40 and 
ELP at the solution-air interface there can be highlighted such elements as follows: water 
molecule and -СНз, -CH»;-, -CH-, «CO and -OH groups. As the Lifshitz-van der Waals 
interactions at the first approximation are equal to the dispersion ones, they depend on the 
type of the groups being in the surfactants molecule, hence their surface tension is related to 
the orientation of surfactant molecules towards the interface. Taking into account this fact 
van Oss and Constanzo [15] stated that the surfactant surface tension is different depending 
on the orientation of its molecules at the surfactant-air interface. Hence, the surface tension 
of surfactant whose molecules are oriented by tail towards the air phase can be called the 
"surfactant tail surface tension" (yr). In the case of the surfactant molecules orientation 
by the hydrophilic group towards the air, its surface tension is called the "surfactant head 
surface tension” (үн). If there are no hydrophilic groups and/or a double bond between 
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the carbon atoms in the tail, then the yr results only from the Lifshitz-van der Waals 
intermolecular interactions. In the case of RH40 the -OH groups present in its tail and in 
the ELP tail apart from the -OH groups also a double bond occurs between the carbon 
atoms. For this reason, the hydrophobic properties of tail are reduced. Unfortunately, the 
literature lacks the data on the surface tension of RH40 and ELP tails and heads. However, 
it is possible that the surface tension of both surfactants tail is the same and close to the 
Lifshitz-van der Waals component of fatty acids. Depending on the type of fatty acid the 
values of "W are larger than the surface tension of hexadecane and smaller than 30 mN/m. 
The values of үн of RH40 and ELP can be approximately compared to that of Triton's, 
particularly Triton X-165 (TX165), which has 16 oxyethylene groups in one molecule. The 
гү!!! value for TX165 is equal to 27.7 mN/m at 293 K [19]. This value is insignificantly higher 
than that of 4!" for the water surface tension (26.85 mN/m) [20], but the electron-acceptor 
and electron-donor parameters of TX165 acid-base component of its head surface tension 
differ considerably from those of water surface tension. It seems that the components and 
parameters of RH40 and ELP surface tension can be close to those of TX165. If so, it can be 
stated that the adsorption of RH40 and ELP mixture at the water-air interface decreases 
only the acid-base component of the water surface tension. However, it is possible that in a 
studied range of mixture concentration the values "W of the surface tension of surfactant 
mixtures solution slightly increase as a result of their molecules tail orientation in the mixed 
monolayer towards the air phase, whereas the reduction in the water surface tension takes 
place due to the decrease in the acid-base component of its surface tension. In the other 
words, the number of hydrogen bonds between the water molecules decreases as a function 
of RH40 and ELP concentration. The greater changes of the number of hydrogen bonds 
take place in the range of surfactant and their mixtures concentration in which the linear 
dependence between the surface tension and log C/ log C12 is observed (Стэ is the mixture 
concentration) (Figure 1 and Ref. [21]). 

According to the Gibbs isotherm equation [4], the saturated mixed monolayer at the 
water-air interface is formed in this range of C12. The question arises why such a significant 
decrease in solution surface tension takes place in this concentration range. Indeed, this 
concentration range is smaller than CMC and it is difficult to take into account the influence 
of the micelles on the solution surface tension. It is more probable that in the saturated 
monolayer the surfactant molecules change the orientation and can change the part of tail 
being in the air phase. This fact causes the changes of the water molecules orientation 
preventing them from the hydrogen bond formation. On the other hand, the concentration 
gradient of surfactants at the interface changes with their concentration changes in the bulk 
phase. As a matter of fact, the destruction of hydrogen bonds is not complete because the 
minimal surface tension of the aqueous solution of RH40 and ELP mixtures is considerably 
higher than that resulting from the surfactants tail (Figure 16-е). The minimal values of the 
surface tension change as a function of the mixture composition (Figure 2). These changes 
are not linear and probably result from the synergetic effect of the micelle formation. It is 
known that only single molecules adsorb at the water-air interface and the surface tension 
of the surfactant aqueous solution depends largely on the concentration of surfactants 
in the monomeric form in the bulk phase. This concentration depends on the tendency 
of the surfactants and their mixture to form micelles. Unfortunately, based only on the 
surface tension isotherms it is difficult to explain more precisely changes in the minimum 
surface tension of the aqueous solutions of the RH40 and ELP mixture as a function of 
its composition. Additional measurements are required. However, we can state that the 
minimal values of the surface tension of the aqueous solution of RH40 and ELP mixtures do 
not differ significantly from the minimal values of the aqueous solution of TX165 at a given 
temperature (39.5 mN/m at 293 K) [22]. However, this value is achieved at a significantly 
smaller concentration than that of TX165. This indicates that the adsorption activity of 
КН40 and ELP as well as of their mixture is larger than that of TX165. 
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Figure 1. A plot of the surface tension (yzy) of the aqueous solutions of RH40 (a) as well as RH40 
and ELP mixtures at the mole fraction of ELP in the bulk phase, «, equal to 0.2 (b), 0.4 (c), 0.6 (d) and 
0.8 (e) vs. the logarithm of their concentration (logC and log C12) at different temperatures equal to 
293, 298, 303, 308, 313 and 318 K, respectively. 
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Figure 2. A plot of the minimal surface tension On) of the aqueous solutions of RH40 and ELP 
mixtures at the temperature equal to 293 (curve 1), 298 (curve 2), 303 (curve 3), 308 (curve 4), 313 
(curve 5) and 318K (curve 6) vs. the mole fraction of ELP in the bulk phase, a. 


The adsorption activity of a surfactant is closely related to the hydration of tail and 
head of surfactant molecules. Tail hydration is the driving force for the surfactant adsorp- 
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tion and the head hydration is the brake. It was proved that the volume of the surfactant 
molecule can be approximated by cubes in which a given part of the surfactant is inscribed. 
Hence, taking into account the length of the bonds and angle between them as well as the 
distance between molecules, it is possible to establish the contactable area of RH40 and 
ELP tail as equal ~967.7 А2. Since the contact plane of the water molecule with a molecule 
of another substance can be assumed to be equal to 10 À?, theoretically about 97 water 
molecules can be contacted with the RH40 or ELP tail simultaneously. On the other hand, 
it is known that one oxygen atom in the oxyethylene group can form a strong hydrogen 
bond with two molecules of water [23]. Assuming that the other oxygen atoms, except the 
one in the -OH group, can also form the hydrogen bond with two molecules of water, then 
the RH40 head can by hydrated by 92 molecules of water and head of ELP by 82 ones. This 
is in accordance with the earlier studies [24]. 

If, during the transfer of RH40 and ELP molecules from the bulk phase to the mixed 
surface layer, the head of their molecules does not exist or is poorly dehydrated, then 
the driving forces of adsorption are large, and therefore RH40 and ELP are more surface 
active than Triton’s [25]. It is possible that the dehydration of their head increases with 
the increasing temperature. On the other hand, the isotherms of the surface tension of the 
aqueous solutions of RH40 and ELP mixtures with different compositions obtained from 
the measurements are almost parallel (Figure 1b—e). 

The surface tension of the aqueous solution of studied mixtures at a given concentra- 
tion and composition as a function of temperature in most cases decreases almost linearly 
(Figure 3 as an example). This indicates that in the studied range of the temperature the 
dehydration of RH40 and ELP is insignificant. 


Surface tension, т, (mN/m) 


305 310 315 
T(K) 


Figure 3. A plot of the surface tension (yzy) of the aqueous solutions of RH40 and ELP mixtures at 
the mole fraction of ELP in the bulk phase equal to 0.2 (curve 1), 0.4 (curve 2), 0.6 (curve 3) and 0.8 
(curve 4) at Сэ = 5 x 1077 M vs. the temperature, T. 


On the other hand, there is not a linear dependence between the surface tension of 
the aqueous solution of RH40 and ELP mixture at a given concentration as a function of 
mixture composition (Figure 4 as an example). The positive and negative deviations of 
үгу values from the linear dependence were observed. This indicates that some changes 
in the hydration number of surfactants head can take place and/or the changes in the 
surfactant molecules packing as a result of their configuration variation. As follows from 
the literature [26,27] the possibility of forming different configuration by surfactants having 
large and branched molecules is greater than for the classical surfactants whose molecules 
are linear and not very long. It is possible that the oxyethylene and/or hydrocarbon 
chains are coiling. Such phenomenon was observed for other surfactants having many 
oxyethylene groups in their molecules [28,29]. 
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Figure 4. A plot of the surface tension (угу) of the aqueous solutions of RH40 and ELP mixtures at 
temperature T equal to 293 K and C or C4? equal to 1075 (curve 1), 5 x 1076 (curve 2), 1075 (curve 3) 
and 5 x 10^? M (curve 4) vs. the mole fraction of ELP in the bulk phase, «. 


From the theoretical and practical points of view it is interesting to describe and/or 
to predict the isotherm of the surface tension of the aqueous solution of the RH40 and 
ELP mixture on the basis of the isotherm of the individual component of the mixture. It 
appeared that ай isotherms of yzy can be described by the exponential function of the 
second order (Figures 51-56 in Supplementary Materials). This function has the form: 


-С е 
YLV = vo + Avexp( ==) 22:16:31 (1) 


where yzy is the surface tension of the aqueous solution of surfactant, C is the surfactant 
concentration, yọ, Ал, Аз, tj and t» are the constants. 

Unfortunately, no unequivocal dependence of the constants in Equation (1) on the 
physicochemical properties of surfactants has been reported so far. It was suggested that 
these constants are related to the components and parameters of the surfactants head and 
tail [30]. However, it is difficult to predict the constants in Equation (1) for the mixtures of 
surfactants based on the constants for individual components of the mixtures. In the case 
of RH40 and ELP mixtures the additional difficulties result from the fact that they can be 
treated not as binary but rather as multicomponent ones. However, it can be stated that 
the values of the yo constant are close to the minimal values of угу at each temperature 
(Figure 57). As the minimal values of yzy depend on the surface tension of surfactant tail, it 
can be stated that yo also is related to the tail of surfactant surface tension. Although the ур 
values are close to the minimum yzy values, the yo values for the RH40 with ELP mixture 
cannot be predicted based on the ур values for RH40 and ELP. This results from the fact that 
there is no linear relationship between the ур values and the composition of the mixture 
(Figure S8). Changes of the remaining constants from Equation (1) are also non-linear as a 
function of the mixture composition (Figures S9-S12). Hence, it can be concluded that the 
isotherms of the surface tension of the aqueous solution of RH40 and ELP mixtures can be 
satisfactorily described but not predicted. 

The changes of the surface tension of the aqueous solutions of many surfactants as a 
function of their concentration can be described by the Szyszkowski equation if the Gibbs 
surface excess concentration of the saturated monolayer at the water-air interface and the 
standard Gibbs free energy of adsorption are known. The Szyszkowski equation for the 
aqueous solution of RH40 and ELP mixtures can be expressed in the form [4]: 


yw — үу = л = —RIT™™ In (S T 1), (2) 
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where yw is the water surface tension, 7t is the monolayer at the water-air interface pressure, 
a is the constant, which depends on the standard Gibbs free energy of adsorption (AG? ds)» 
Ie" is the maximal Gibbs surface excess concentration of surfactant in the saturated 
monolayer at the water-air interface, T is the temperature and R is the gas constant. 

The a constant in Equation (2) fulfills the expression [4]: 


AG? 
а = w exp E | , (3) 


where: « is the number of the water moles in one dm?. 

To describe the surface tension isotherm Equation (2) can be solved numerically 
choosing proper values of Г” and a. For calculation of үгу from Equation (2) it was 
assumed that C changes from 0 to CMC. This resulted from the fact that not aggregated 
surfactant molecules are surface active. It appeared that at the first approximation it 
was possible to describe all isotherms of surface tension using the Szyszkowski equation 
(Equation (2)) (Figures 51-56). 

For the ideal mixture of surfactants, it is possible not only to describe but also to 
predict the isotherm of the aqueous solutions of this mixture at different compositions 
based on the data for individual surfactants. In the case of the ideal surfactant mixture 
the changes of the I"7*, a and CMC as a function of its composition should be linear. In 
addition, the surface tension of the aqueous solution of ideal surfactants mixture should 
satisfy the equation which for the binary mixture has the form: 


YW = YivX? + fiv X5, (4) 


where yty and ATW is the surface tension of the aqueous solution of surfactants 1 and 2 at a 
given concentration in the bulk phase and X and X3 are the mole fractions of surfactants 1 
and 2 in the monolayer. 

As it was suggested earlier, X? and X3 can be determined using the film pressure 
of surfactants 1 and 2 at their given concentration in their individual solutions. Hence, 
Are mau and © = m47; Un and л» are the layer of surfactants 1 and 2 pressure, 
respectively). The values of yzy calculated from Equation (4) for all studied systems are 
higher than the measured ones at the same RH40 and ELP concentration (Figures 52-55). 
However, the difference between the measured and calculated values of угу depends on 
the composition of this mixture. In the case of RH40 and ELP mixture with the mole fraction 
of ELP equal to 0.6, the calculated values of үгу are almost identical to those measured 
(Figure S4). For the aqueous solution of RH40 and ELP at a given composition the difference 
between the calculated from Equation (4) and measured values of угу practically does not 
change as a function of temperature (Figures 52-55). 

The difference between the measured and calculated values of угу indicates that the 
RH40 and ELP mixture does not behave as ideal one and the deviation from the ideal 
behaviour depends on the mixture composition. This conclusion is also confirmed by 
the I"?* values used in the Szyszkowski equation (Equation (2)) for the yzy calculations. 
The Г" values do not change linearly as a function of mixture surfactants composition 
(Figure 513), and they are in the range from 2.44 х 1076 to 3.12 x 1076 mol/m?. These 
values correspond to the minimal area occupied by one molecule which is in the range from 
53.2 to 68 À?. Thus, it can be concluded that there is a synergetic effect in the reduction in 
water surface tension by the RH40 and ELP mixture which depends on its composition. 
For this reason, for the RH40 and ELP mixture the surface tension of its aqueous solution 
cannot be precisely predicted based on the I"^7*, a and CMC for RH40 and ELP. However, 
this tension can be predicted using the Fainerman and Miller equation. This equation for 
binary mixtures has the form [31,32]: 


expll = expll; + expllI, — 1, (5) 
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where П = zto/ RT, II = 7104/RT and П, = 712@2/KT are the dimensionless pressure 
of the mixed monolayer at the water-air interface and individual surfactants 1 and 2, 
respectively, and @1, @2 and © are the areas occupied by one mole of surfactants 1 and 2 
and mixture at the water-air interface and л, лу and л» are the differences between the 
surface tension of the solvent and solution of the surfactants mixture and components 1 
and 2, respectively. 

The calculation made using Equation (5) confirmed the above mentioned conclusion 
that, on the basis of the Fainerman and Miller equation, the isotherm of surface tension of 
RH40 and ELP mixture can be predicted (Figures 52-55). However, at the concentration 
of the surfactant mixtures above their CMC the values calculated from Equation (5) are 
smaller than those measured. Some differences between the measured and calculated from 
Equation (5) угу values are observed for the mixture with the mole fraction of ELP equal 
to 0.6 at the constant temperatures 313 and 318 K (Figure 54). 


2.2. Concentration of Surfactants at the Water—Air Interface and Composition of the 
Mixed Monolayer 

In most cases the concentration of surfactants and/or their mixture is determined 
using the Gibbs adsorption isotherm Equation [4]. As a matter of fact, this equation allows 
us to determine only the surface excess concentration of surfactants at the interface, but 
it can be treated as the total dimensional concentration because of the small surfactant 
concentration in the bulk phase. For the studies on the adsorption properties of surfactants 
their very low concentration in the aqueous solution is used. Therefore, for such solution 
it can be assumed that the activity coefficient of surfactants is close to unity and its mole 
fraction is very close to € and w is practically constant in the range of studied surfactant 
concentrations. In such case the Gibbs isotherm equation for the nonionic surfactants and 
their mixtures has the form [4]: 


pe’ (22) = 1 ( OYLV ) (6) 
RTĄ OC jr 2.303RT V dlogC / т. 

The calculated values of I indicate that the adsorption of RH40 and ELP and their 
mixtures depends on the temperature and in the case of mixtures also on their composition 
(Figure 514). However, for the RH40 and ELP mixtures in which the mole fraction of ELP is 
equal to 0.2 and 0.8, respectively, the surface excess concentration at the water-air interface 
decreases and increases as a function of temperature. This indicates that for these mixtures 
the synergetic and/or antagonistic effect in the reduction in water surface tension takes 
place. It is also confirmed by the fact that there is no linear dependence between I"?* and 
composition of the RH40 and ELP mixtures. Thus, the negative and positive deviations 
from the linear dependence are observed. In many cases the Г"! values obtained from the 
Gibbs isotherm equation using the linear dependence between the surface tension of the 
solution and the logarithm from the surfactants concentration are similar to the Г” values 
obtained from the Szyszkowski equation (Figure $13). As mentioned above this equation 
was solved against the surface tension of the aqueous solution of surfactant numerically. 
The calculations of Г from Equation (6) also indicate that in many cases the values of 
are obtained for RH40 and ELP as well as their mixtures at the concentration equal to 
1 x 1076 mol/dm?. It was proved that in the range of the surfactants and their mixture 
concentrations from 0 to 1 x 107 mol/dm?, the independent adsorption at the water-air 
interface takes place (Figures 52-55). This results from the calculation of уру from the 
following expression: 


pmax 


YLV = Tw — 74 — 70, (7) 


In the mentioned range of surfactants concentration, the values of yzy calculated from 
Equation (7) are very close to those measured (Figures 52-55). This points out that in this 


1 б. S лә : ЯМ 
range of surfactants concentration Хү = ma and X5 = = frz: The relative composition 
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of the surfactants saturated mixed monolayer at the water-air interface can be determined 
using the Hua and Rosen equation which has the form [4,33]: 


(X$'In (aC42/ X$C1) 
(1 — X$)"In[(1 — а)С/ (1 - Xf) Co] 


where « is the mole fraction of surfactant 1 (ELP) in the mixture in the bulk phase, and 
indices 1, 2 and 12 are related to surfactants 1, 2 and their mixture, respectively. 

Taking into account the error in the calculation of mole fraction of surfactants in the 
mixed monolayer using Equation (8), it can be stated that in most cases the X? and X$ values 
determined from this equation are close to those determined from 7t, and л» (Figure 515). 
Hence, the values of x? and 6 obtained from 7t; and л», at the first approximation can 
be treated as a mole fraction of RH40 and ELP in the mixed monolayer at the water-air 
interface in the whole range of RH40 and ELP mixture concentration in the bulk phase. 

In the range of ELP mixture concentration corresponding to the unsaturated mixed 
monolayer, i.e., from zero to 1 x 1076 mol/dm? in the bulk phase the mutual effect of 
surfactants on their adsorption was not observed and for each composition of the surfactant 
mixtures the mole fraction of ELP in the mixed monolayer decreases in this concentration 
range of the mixture in the bulk phase. However, from the beginning of the saturated 
mixed monolayer formation at the water-air interface the X? increases to the value which 
depends on the composition of the mixture in the bulk phase. In any case, the mole fraction 
of X? is equal to a. However, the difference between these fractions does not yet prove the 
existence of synergy or antagonism in the reduction in the water surface tension by the 
mixture of RH40 and ELP. Hua and Rosen [33] stated that the synergy and antagonistic 
effects in the reduction in water surface tension by the adsorption of surfactants mixture 
can be determined on the basis of the parameter of intermolecular interactions (B^). This 
parameter fulfills the condition [4,34]: 


In (xC15/ X?C1) 
(1- xf) 


p = (9) 


From the calculations of the 6” parameter from Equation (9) it results that each value 
depends on the concentration and composition of the RH40 and ELP mixture as well as 
temperature (Figure S16). This is connected with the changes of the mixed monolayer 
composition at the water-air interface as a function of mixture concentration and temper- 
ature at a given constant mixture composition in the bulk phase, a. It should be noted 
that the values of the surfactant mole fraction in the mixed monolayer calculated from 
Equation (8), as mentioned above, are close to those determined on the basis of the surface 
tension isotherms of RH40 and ELP. The B^ parameter can assume both the positive and 
negative values (Table 51 in Supplementary Materials). For the mixture of RH40 and ELP 
with the ELP mole fraction equal to 0.6 the values of B^ are only positive and there is the 
almost linear dependence between В° and the surface tension of solutions and temperature 
(Figure S16d). This indicates the synergetic effect in the water surface tension does not 
occur for this mixture. For the other RH40 and ELP mixtures this effect was noted but not 
at each temperature applied in the studies. As mentioned above, in the case of the RH40 
and ELP mixture, the synergetic effect in the reduction in water surface tension results 
probably from the changes of the configuration of RH40 and ELP molecules in the mixed 
monolayer, in comparison to the individual surfactant monolayer and/or the changes of 
number of hydrogen bonds as a result of dehydration of surfactants tail and head. 


2.3. Thermodynamic Parameters of Adsorption 


The standard Gibbs free energy (AGO, ,), enthalpy (AH?) and entropy (AS? ds) of 
adsorption are powerful in understanding the adsorption process of the RH40 and ELP 


mixture at the water-air interface [4]. АСО, informs only about the surfactants tendency to 
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adsorb at the water-air interface. Based on this it is not possible to conclude what changes 
took place in the bulk phase as a result of surfactants adsorption. The AH? 4, Value informs 
us about chemical reactions which occurred during the adsorption process. In other words, 
based on AH9, , it can be stated whether more chemical bonds were broken or whether new 
ones were created. On the other hand, AS? ds informs about the structural and orientation 
changes of molecules in the bulk phase during the adsorption process. Л50, is the main 
force of this process. 

The literature reports numerous methods for determination of AG? ds" Among them, 
those that are based on the Langmuir equation modified by de Boer are very often used. 
This equation has the form [35,36]: 


A? AO С -А60, 
= j 10 
ААРА А0 cay RT v) 
where A and A? are the areas occupied by one molecule in the monolayer and limiting one. 


АС? can be also determined from the linear Langmuir equation [4,35]: 


C. © 8 
T [max ' pmax’ 


(11) 


where a fulfills Equation (4). 
The AG? js Values were determined from Equations (10) and (3). In Equation (3) there 
were used the values of a established from the Szyszkowski equation and obtained from 
the Langmuir linear equation. Hence, for each system three values of AG0,, were ob- 
tained. There are small differences between the AG. values obtained using these methods 
(Table S2 in Supplementary Materials). There is no linear dependence between АСО, and 
the composition of RH40 and ELP mixtures. The absolute values of AG? 4, for RH40, ELP 
and their mixtures are higher than those for other nonionic surfactants. This explains why, 
at small concentrations of RH40, ELP and their mixtures, the reduction in water surface 
tension by their adsorption is higher than that of other nonionic surfactants [4,22]. 
It is known that AG, fulfills the expression [4,35]: 
AG), = AH, = TAS? (12) 


айв ads’ 


Assuming that in the temperature range from 293 to 318 K the АН, is constant 


ads 
then [4,35]: 


dAG? 
a = – А504, (13) 
0 


Introducing the values of AS 
values of AH? were determined. 

It follows from Table 1 that the values of AH’, for the RH40 and ELP mixture with 
the mole fraction of 0.6 are close to zero. For this mixture no synergetic effect is observed 
in the reduction in water surface tension. For the mixture of different compositions, the 
positive and negative values of AH? , Were obtained. However, the absolute value of AH?, s 
is not large. 

The nonlinear dependence between the AG0,, and composition of the RH40 and ELP 
mixture can result from the fact that the Gibbs free energy of RH40 and ELP mixing is 
different from zero. In the other words, the mixed monolayer is not ideal and the coefficients 
of RH40 and ELP activity are different from zero (Table S1). The coefficient of surfactants 
activity and next the Gibbs free energy of surfactants mixing can be determined based on 
the parameter of the intermolecular interactions. The coefficients of surfactants activity 
satisfy the condition [4]: 


calculated from Equation (13) to Equation (12) the 


ads 


Inf, = g^ (1 — X$}, (14) 
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and 
Info = В" (Xi). (15) 


In turn the Gibbs free energy of surfactants mixing (GE,,) can be expressed by equa- 
tion [37]: 
Ст = RTQGIn fi + Хп), (16) 


It was proved that the absolute values of СЁ, are not great. This means that there is 


no great difference between the ideal and real behavior of the RH40 and ELP mixtures. 


Table 1. The values of the AH), Е апа AS? ds for the RH40 and ELP as well as their binary mixtures 


calculated from Equations (12) and (13), respectively. 


450 [kJ /molK] 


RH40 a =0.2 х= 04 a = 0.6 х= 0.8 ЕІР 
—0.172 —0.178 —0.191 —0.146 —0.159 —0.161 
АН0 [kJ/mol] 

T [K] RH40 а-02 x=0.4 x=0.6 x=0.8 ELP 
293 3.344 5.034 8.984 —3.262 —0.197 0.857 
298 3.403 5.123 9.140 —2.983 —0.201 0.871 
303 3.473 5.221 9.305 —3.034 —0.196 0.896 
308 3.292 5.069 9.101 —3.084 —0.200 0.910 
313 3.352 5.037 8.977 —3.125 —0.195 0.845 
318 3.411 5.116 9.123 —3.166 —0.199 0.869 


3. Materials and Methods 


For our studies the aqueous solution of Kolliphor® ELP (ELP) (Cremophor® ELP, 
Polyoxyl 35 Hydrogenated Castor oil, Polyoxyl-35 castor oil, CAS number 61791-12-6,) and 
Kolliphor® RH 40 (КН40) (Cremophor® RH 40, macrogolglycerol hydroxystearate, PEG-40 
castor oil, Polyoxyl 40 hydrogenated castor oil, CAS number 61788-85-0,) mixture were 
applied. RH40 and ELP were purchased form Sigma-Aldrich (St. Louis, MO, USA) and 
used without additional purification. The solution concentration was in the range from 0 
to the values significantly higher than the critical micelle concentration (CMC) of RH40 
and ELP that is from 0 to 0.01 M (mol/dm?) and the ELP mole fraction in the surfactant 
mixtures in the bulk phase, a, was equal to 0.2, 0.4, 0.6 and 0.8. The water used for the 
solutions preparation was doubly distilled and deionized (Destamat) the surface tension of 
which changed from 72.8 to 68.7 mN/m in the temperature range from 293 to 318 K. 

The surface tension (yy) measurements of the aqueous solution of the RH40 and 
ELP mixture were made in the temperature range from 293 to 318 K using the Krüss 
K100 tensiometer according to the platinum ring tensiometer method (du Nouy's method) 
calibrated before the measurements. The calibration was made only at 293 K using water 
and methanol whose surface tension at this temperature was equal to 72.8 and 22.5 mN/m, 
respectively. The surface tension measurements for each concentration and composition 
of the aqueous solution of RH40 and ELP mixtures were repeated at least ten times. The 
standard deviation of the results obtained from the measurements was + 0.1 mN/m, and 
the uncertainty was in the range from 0.3% to 0.9%. 


4. Conclusions 


Based on the results obtained from the surface tension measurements of the aqueous 
solution of RH40 and ELP mixtures and their analysis many conclusions can be drawn. 

In the small concentration range of RH40 and ELP as well as their mixture in the 
aqueous solution the greater reduction in water surface tension takes place as a result 
of the adsorption at the water-air interface that due to the adsorption of other nonionic 
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surfactants. However, the minimal surface tension of the aqueous solution of RH40 and 
ELP and their mixtures which can be obtained is close to that of TX165. 

There is no linear dependence between the surface tension of RH40, ELP mixtures and 
their composition. The isotherm of the surface tension of RH40 and ELP, as well as their 
mixture, can be precisely described by the exponential function of the second order. It is 
possible to describe the isotherm of the surface tension by the numerical solution of the 
Szyszkowski equation. The maximal Gibbs surface excess concentration obtained solving 
the Szyszkowski equation against the surface tension is comparable to that determined 
from the equation of the Gibbs isotherm of adsorption. 

The surface tension of the aqueous solution of RH40 and ELP mixtures can be predicted 
using the Fainerman and Miller equation in the range of mixtures concentration from zero 
to CMC. In some cases, it is possible to predict the changes of the surface tension of the 
aqueous solution of RH40 and ELP mixtures in the whole studied range of concentration. 

The synergetic effect of the RH40 and ELP mixtures in the reduction in the water 
surface tension but not at all their composition was found based on the intermolecular 
interactions parameter. 

The composition of the RH40 and ELP mixed monolayer at the water-air interface 
can be predicted from the contribution of a particular surfactant to the reduction in water 
surface tension. 

The calculated values of the mole fraction of RH40 and ELP in the mixed monolayer 
at the water-air interface, on the basis of their aqueous solution surface tension isotherms, 
are in most cases close to those determined based on the Rosen and Hua concept. 

The values of the Gibbs free energy of adsorption for a given system calculated from 
the constant in the Szyszkowski equation are close to those calculated from the Langmuir 
equation modified by de Boer and the linear Langmuir equation. 

The absolute values of the Gibbs free energy of adsorption for the studied systems are 
higher than those for Triton’s. 

The standard enthalpy of adsorption assumes negative and positive values depending 
on the composition of the RH40 and ELP mixture, but the absolute values of enthalpy do 
not differ significantly from zero. 


Supplementary Materials: The following are available online, Scheme S1: The structure of the 
main components of ELP (a) and КН40 (b); Table 51: The values of the mole fractions of ELP (X?) 
and RH40 (X3) in the mixed monolayer at the water-air interface, parameter of intermolecular 
interactions (80), activity coefficient of ELP (f,) and RH40 (f2) as well as and Gibbs free energy 
of surfactants mixing (GF) (kJ/mol) for the mixtures at the mole fraction of ELP in the bulk 
phase, « equal to 0.2 (a), 0.4 (b), 0.6 (c) and 0.8 (d). There is also the condition for synergism or 
antagonism existence (In(C1,C2)) [4]; Table S2: The values of the standard Gibbs free energy of 
adsorption (AG?) calculated from Equations (3), (10) and (11); Figure 51: A plot of the surface 
tension (Үд) of the aqueous solutions of RH40 vs. the logarithm of its concentration (log C) at 
different temperatures equal to 293 (a), 298 (b), 303 (c), 308 (d), 313 (e) and 318 K (f). Points 1 
correspond to the measured values, curves 2 and 3 correspond to the values calculated from the 
Szyszkowski equation (Equation (2)) and the exponential function of the second order (Equation (1)), 
respectively; Figure 52: A plot of the surface tension (yzy) of the aqueous solutions of RH40 and 
ELP mixtures at the mole fraction of ELP in the bulk phase equal to 0.2 vs. the logarithm of 
their concentration (log C12) at different temperatures equal to 293 (a), 298 (b), 303 (c), 308 (d), 
313 (e) and 318 K (f). Points 1 correspond to the measured values, curves 2-6 correspond to the 
values calculated from the Szyszkowski equation (Equation (2)), exponential function of the second 
order (Equation (1)), Fainerman and Miller equation (Equation (5)), Equation (7) and Equation (4), 
respectively; Figure 53: A plot of the surface tension (yzy) of the aqueous solutions of RH40 and 
ELP mixtures at the mole fraction of ELP in the bulk phase equal to 0.4 vs. the logarithm of 
their concentration (log C12) at different temperatures equal to 293 (a), 298 (b), 303 (c), 308 (d), 
313 (e) and 318 К (f). Points 1 correspond to the measured values, curves 2-6 correspond to the 
values calculated from the Szyszkowski equation (Equation (2)), exponential function of the second 
order (Equation (1)), Fainerman and Miller equation (Equation (5)), Equation (7) and Equation (4), 
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respectively; Figure 54: A plot of the surface tension (yzy) of the aqueous solutions of RH40 and 
ELP mixtures at the mole fraction of ELP in the bulk phase equal to 0.6 vs. the logarithm of 
their concentration (log C12) at different temperatures equal to 293 (a), 298 (b), 303 (c), 308 (d), 
313 (e) and 318 K (f). Points 1 correspond to the measured values, curves 2-6 correspond to the 
values calculated from the Szyszkowski equation (Equation (2)), exponential function of the second 
order (Equation (1)), Fainerman and Miller equation (Equation (5)), Equation (7) and Equation (4), 
respectively; Figure 55: A plot of the surface tension (yzy) of the aqueous solutions of RH40 and 
ELP mixtures at the mole fraction of ELP in the bulk phase equal to 0.8 vs. the logarithm of 
their concentration (log C12) at different temperatures equal to 293 (a), 298 (b), 303 (c), 308 (d), 
313 (e) and 318 K (f). Points 1 correspond to the measured values, curves 2-6 correspond to the 
values calculated from the Szyszkowski equation (Equation (2)), exponential function of the second 
order (Equation (1)), Fainerman and Miller equation (Equation (5)), Equation (7) and Equation (4), 
respectively; Figure 56: A plot of the surface tension (yzy) of the aqueous solutions of ELP vs. the 
logarithm of its concentration (log C) at different temperatures equal to 293 (a), 298 (b), 303 (c), 308 
(d), 313 (e) and 318 K (f). Points 1 correspond to the measured values, curves 2 and 3 correspond to 
the values calculated from the Szyszkowski equation (Equation (2)) and the exponential function 
of the second order (Equation (1)), respectively; Figure 57: A plot of the values of constant уо in 
Equation (1) (curve 1) and the minimal surface tension of aqueous solution (үгэн) (curve 2) vs. the 
temperature (T) for the RH40 and ELP mixtures at the mole fraction of ELP in the bulk phase equal 
to 0 (RH40 (a)), 0.2 (b), 0.4 (c), 0.6 (4), 0.8 (е) and 1 (ELP (f)), respectively; Figure 58: A plot of the 
values of constant ур in Equation (1) for studied surfactant mixtures vs. the mole fraction of ELP in 
the bulk phase (а) at the temperatures equal to 293 (curve 1), 298 (curve 2), 303 (curve 3), 308 (curve 
4), 313 (curve 5) and 318 K (curve 6), respectively; Figure 59: A plot of the values of constant A, in 
Equation (1) for studied surfactant mixtures vs. the mole fraction of ELP in the bulk phase (а) at the 
temperatures equal to 293 (curve 1), 298 (curve 2), 303 (curve 3), 308 (curve 4), 313 (curve 5) and 318 
K (curve 6), respectively; Figure 510: A plot of the values of constant A> in Equation (1) for studied 
surfactant mixtures vs. the mole fraction of ELP in the bulk phase («) at the temperatures equal to 293 
(curve 1), 298 (curve 2), 303 (curve 3), 308 (curve 4), 313 (curve 5) and 318 K (curve 6), respectively; 
Figure 511: A plot of the values of constant tı in Equation (1) for studied surfactant mixtures vs. the 
mole fraction of ELP in the bulk phase (4) at the temperatures equal to 293 (curve 1), 298 (curve 
2), 303 (curve 3), 308 (curve 4), 313 (curve 5) and 318 K (curve 6), respectively; Figure 512: A plot 
of the values of constant t» in Equation (1) for studied surfactant mixtures vs. the mole fraction of 
ELP in the bulk phase (x) at the temperatures equal to 293 (curve 1), 298 (curve 2), 303 (curve 3), 308 
(curve 4), 313 (curve 5) and 318 K (curve 6), respectively; Figure 513: The values of Г"! calculated 
from Equations (6) and (2) for studied surfactant mixtures at the mole fraction of ELP in the bulk 
phase (д) equal to 0, 0.2, 0.4, 0.6, 0.8 and 1 at different temperatures equal to 293 (a), 298 (b), 303 (c), 
308 (4), 313 (e) and 318 К (f), respectively; Figure 514: A plot of the values of Gibbs surface excess 
concentration (I) vs. logarithm of the concentration in the bulk phase (C or C12) for the RH40 and 
ELP mixtures at the mole fraction of ELP in the bulk phase equal to 0 (КН40 (a)), 0.2 (b), 0.4 (с), 0.6 (4), 
0.8 (e) and 1 (ELP (f)), respectively; Figure S15: The values of the mole fraction of ELP in the mixed 
monolayer at the water air interface (X ) calculated from the relationship: x$ = Rom (bars 1) and 
from Equation (8) (bars 2) and temperature range 293-318K for mixtures at the mole fraction of ELP 
in the bulk phase (4) equal to 0.2 (a), 0.4 (b), 0.6 (c) and 0.8 (d). Figure 516: A plot of the values of 
parameter of intermolecular interactions, £^, for studied binary surfactant mixtures (a) at the mole 
fraction of ELP in the bulk phase л, equal to 0.2, 0.4, 0.6, 0.8 and T = 293K vs. surface tension, (үгу), 
(b) at the surface tension, (yry), equal to 65, 55, 45 mN/m and T = 293K vs. mole fraction of ELP in 
the bulk phase a, (c) at the mole fraction of ELP in the bulk phase «, equal to 0.2, 0.4, 0.6, 0.8 and 
YLv = 65 mN/m vs. temperature, T, (d) at the mole fraction of ELP in the bulk phase л, equal to 0.6 
and үгу = 65, 55, 45 mN/m vs. temperature, T. 
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Scheme S1. The structure of the main components of ELP (a) and RH40 (b). 


Table S1. The values of the mole fractions of ELP (X7) and RH40 (X$) in the mixed monolayer at the water-air interface, 
parameter of intermolecular interactions (8°), activity coefficient of ELP (f1) and RH40 (№) as well as and Gibbs free 
energy of surfactants mixing (С2,,) (kJ/mol) for the mixtures at the mole fraction of ELP in the bulk phase, a equal to 


0.2 (a), 0.4 (b), 0.6 (c) and 0.8 (d). There is also the condition for synergism or antagonism existence (In( C1,C7)) [4]. 


a) 

a ELP = 0.2, үу = 65 mN/m 
T [K] xi Xj B^ In( 6,,C) fi fz Gmix 
293 0.2708 0.7292 -1.7006 0.3838 0.4048 0.8828 -0.8180 
298 0.2603 0.7397 -1.6267 0.4380 0.4106 0.8956 -0.7760 
303 0.2365 0.7635 -1.4850 0.4960 0.4208 0.9203 -0.6755 
308 0.2405 0.7595 -1.4010 0.5388 0.4456 0.9222 -0.6552 
313 0.2518 0.7482 -1.2870 0.3152 0.4865 0.9216 -0.6310 
318 0.2673 0.7327 -1.2858 0.2217 0.5014 0.9122 -0.6658 

a ELP = 02, y,y = 55 mN/m 

293 0.2722 0.7278 -0.8054 0.2215 0.6527 0.9421 -0.3887 
298 0.2610 0.7390 -0.2582 0.3000 0.8685 0.9826 -0.1234 
303 0.2400 0.7600 0.3737 0.3797 1.2409 1.0218 0.1717 
308 0.2510 0.7490 0.0917 0.4241 1.0528 1.0058 0.0441 
313 0.2568 0.7432 -1.2022 0.2597 0.5148 0.9238 -0.5971 
318 0.2856 0.7144 -1.4203 0.1393 0.4844 0.8906 -0.7662 

aELP=0.2,y,y = 45 mN/m 
293 0.27682 0.72318 0.27178 0.5473 1.1527 1.0210 0.1325 
298 0.2680 0.7320 0.72404 0.3283 1.4740 1.0534 0.3519 
303 0.2500 0.7500 2.42876 0.1432 3.9203 1.1639 1.1472 
308 0.2590 0.7410 -0.20318 0.0836 0.8944 0.9865 -0.0999 
313 0.3200 0.6800 -1.82648 0.1750 0.4297 0.8294 -1.0343 
318 0.34357 0.65643 -1.77151 0.1837 0.4661 0.8113 -1.0563 

b) 

a ELP = 0.4, Y;y = 65 mN/m 
T [K] xj x5 В? In(C1,Cz) f1 f2 Gmix 
293 0.3267 0.6733 -0.1910 0.3838 0.9171 0.9798 -0.1023 
298 0.2837 0.7163 0.1910 0.4380 1.1029 1.0155 0.0962 
303 0.2935 0.7065 -0.0529 0.4960 0.9739 0.9955 -0.0276 
308 0.3312 0.6689 -0.7090 0.5388 0.7282 0.9252 -0.4021 
313 0.3854 0.6146 -1.7359 0.3152 0.5191 0.7727 -1.0700 
318 0.4333 0.5667 -2.6775 0.2217 0.4232 0.6050 -1.7382 

a ELP = 0.4, Yw = 55 mN/m 

293 0.3801 0.6199 -0.5749 0.2215 0.8018 0.9203 -0.3300 
298 0.3532 0.6468 -0.3420 0.3000 0.8667 0.9582 -0.1936 
303 0.3334 0.6666 -0.2679 0.3797 0.8878 0.9707 -0.1500 
308 0.3397 0.6603 -0.5233 0.4241 0.7960 0.9414 -0.3006 
313 0.3950 0.6050 -1.1453 0.2597 0.6576 0.8363 -0.7122 
318 0.4287 0.5713 -1.8049 0.1393 0.5548 0.7177 -1.1687 

aELP=0.4,y,, = 45 mN/m 
293 0.5260 0.4740 -0.7231 0.5473 0.8501 0.8187 -0.4392 
298 0.4873 0.5127 -0.7825 0.3283 0.8141 0.8304 -0.4844 
303 0.4536 0.5464 -0.8350 0.1432 0.7793 0.8422 -0.5213 
308 0.4437 0.5563 -0.8497 0.0836 0.7688 0.8460 -0.5370 
313 0.4630 0.5370 -1.1080 0.1750 0.7265 0.7886 -0.7169 
318 0.4664 0.5336 -1.2825 0.1837 0.6941 0.7566 -0.8438 


a ELP = 0.6, y, = 65 mN/m 


4) 


TIK] xs xs p^ In( C1,C;) fi fa бин 
293 0.5060 0.4940 0.2101 0.3838 1.0526 1.0553 0.1279 
298 0.4908 0.5092 0.2404 0.438 1.0643 1.0596 0.1489 
303 0.4733 0.5267 0.3271 0.4960 1.0950 1.0760 0.2054 
308 0.4679 0.5321 0.3505 0.5388 1.1043 1.0798 0.2235 
313 0.5201 0.4799 0.3944 0.3152 1.0951 1.1126 0.2562 
318 0.5583 0.4417 0.4225 0.2217 1.0859 1.1407 0.2755 
a ELP = 0.6, y,y = 55 mN/m 
293 0.6090 0.3910 1.1888 0.2215 1.1993 1.5541 0.6896 
298 0.5933 0.4067 1.4586 0.3000 1.2729 1.6710 0.8720 
303 0.5660 0.4340 1.6555 0.3797 1.3659 1.6995 1.0244 
308 0.5580 0.4420 1.7510 0.4241 1.4079 1.7250 1.1059 
313 0.6401 0.3599 1.8977 0.2597 1.2787 2.1761 1.1377 
318 0.7400 0.2600 2.0000 0.1393 1.1448 2.9898 1.0174 
a ELP = 0.6,y,y = 45 mN/m 
293 0.6400 0.3600 1.5235 0.5473 1.2183 1.8664 0.8551 
298 0.6300 0.3700 1.8585 0.3283 1.2897 2.0910 1.0553 
303 0.6400 0.3600 1.9591 0.1432 1.2890 2.2310 1.1371 
308 0.6500 0.3500 2.0930 0.0836 1.2923 2.4213 1.2193 
313 0.7000 0.3000 2.3117 0.1750 1.2313 3.1041 1.2633 
318 0.8000 0.2000 2.3452 0.1837 1.0983 4.4859 0.9921 
a ELP = 0.8, үүү = 65 mN/m 
T [K] X] Х3 p^ In( C,,C;) fi fz Ghix 
293 0.6783 0.3217 -0.7193 0.3838 0.9283 0.7182 -0.3823 
298 0.6926 0.3074 -0.3532 0.4380 0.9672 0.8442 -0.1863 
303 0.7063 0.2937 -0.0337 0.4960 0.9971 0.9833 -0.0176 
308 0.6807 0.3193 0.2259 0.5388 1.0233 1.1103 0.1257 
313 0.7098 0.2902 0.5465 0.3152 1.0471 1.3170 0.2929 
318 0.8291 0.1709 0.6282 0.2217 1.0185 1.5400 0.2353 
a ELP = 0.8, y;y = 55 mN/m 
293 0.6831 0.3169 -1.0828 0.2215 0.8970 0.6033 -0.5710 
298 0.7016 0.2984 -0.8880 0.3000 0.9240 0.6459 -0.4606 
303 0.6893 0.3107 -0.5618 0.3797 0.9472 0.7657 -0.3031 
308 0.6947 0.3053 -0.3522 0.4241 0.9677 0.8437 -0.1912 
313 0.7883 0.2117 0.3289 0.2597 1.0148 1.2267 0.1428 
318 0.9298 0.0702 1.5557 0.1393 1.0077 3.8385 0.2683 
a ELP = 08, y,y = 45 mN/m 
293 0.7063 0.2938 -2.5609 0.5473 0.8017 0.2788 -1.2942 
298 0.6929 0.3071 -2.3356 0.3283 0.8023 0.3258 -1.2314 
303 0.6925 0.3075 -1.8596 0.1432 0.8388 0.4099 -0.9975 
308 0.7251 0.2749 -1.1108 0.0836 0.9195 0.5577 -0.5670 
313 0.7627 0.2374 -0.7509 0.175 0.9586 0.6461 -0.3537 
318 0.8254 0.1746 -0.0268 0.1837 0.9992 0.9819 -0.0102 


Table S2. The values of the standard Gibbs free energy of adsorption (AGQ,,) calculated from Eqs. (3), (10) and (11). 


Temperature 0 
AGads [kJ/mol] 
[KI RH40 a=0.2 а = 0.4 а = 0.6 а = 0.8 ELP 
Eq. (3) 
293 -47.02 -47.02 -47.02 -45.99 -46.81 -46.29 
298 -47.82 -47.82 -47.82 -46.44 -47.61 -47.08 
303 -48.61 -48.61 -48.61 -47.22 -48.40 -47.86 
308 -49.65 -49.65 -49.77 -48.00 -49.20 -48.65 
313 -50.45 -50.57 -50.85 -48.77 -49.99 -49.52 
318 -51.25 -51.38 -51.66 -49.54 -50.79 -50.30 
Eq. (10) 
293 -46.12 -45.41 -45.90 -45.03 -44.86 -45.00 
298 -46.86 -47.08 -46.68 -45.81 -45.81 -45.72 
303 -47.75 -48.06 -47.36 -46.38 -46.66 -46.53 
308 -48.46 -49.19 -48.50 -47.23 -47.41 -47.05 
313 -49.03 -49.76 -50.04 -47.68 -47.89 -48.33 
318 -49.49 -50.90 -51.19 -48.32 -48.39 -49.14 
Eq. (11) 
293 -49.20 -48.14 -48.36 -48.21 -46.88 -47.68 
298 -49.96 -51.11 -49.18 -49.15 -48.03 -48.62 
303 -50.89 -52.07 -49.99 -49.85 -48.88 -49.55 
308 -51.65 -53.23 -51.36 -50.85 -49.73 -50.03 
313 -52.31 -53.39 -53.15 -51.38 -50.27 -51.12 


318 -52.82 -54.72 -54.57 -52.01 -50.82 -52.37 
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Figure 51. A plot of the surface tension (ууу) of the aqueous solutions of RH40 vs. the logarithm of its concentration 
(logC) at different temperatures equal to 293 (a), 298 (b), 303 (c), 308 (d), 313 (e) and 318 K (f). Points 1 correspond 
to the measured values, curves 2 and 3 correspond to the values calculated from the Szyszkowski equation (Eq. (2)) and 
the exponential function of the second order (Eq. (1)), respectively. 
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Figure S2. A plot of the surface tension (у;у) of the aqueous solutions of RH40 and ELP mixtures at the mole fraction of 
ELP in the bulk phase equal to 0.2 vs. the logarithm of their concentration (logC,,) at different temperatures equal to 
293 (a), 298 (b), 303 (c), 308 (d), 313 (e) and 318 K (f). Points 1 correspond to the measured values, curves 2 - 6 
correspond to the values calculated from the Szyszkowski equation (Eq. (2)), exponential function of the second order 
(Eq. (1)), Fainerman and Miller equation (Eq. (5)), Eq. (7) and Eq. (4), respectively. 
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Figure S3. A plot of the surface tension (у;у) of the aqueous solutions of RH40 and ELP mixtures at the mole fraction of 
ELP in the bulk phase equal to 0.4 vs. the logarithm of their concentration (logC;;) at different temperatures equal to 
293 (a), 298 (b), 303 (c), 308 (d), 313 (e) and 318 К (f). Points 1 correspond to the measured values, curves 2 - 6 
correspond to the values calculated from the Szyszkowski equation (Eq. (2)), exponential function of the second order 
(Eq. (1)), Fainerman and Miller equation (Eq. (5)), Eq. (7) and Eq. (4), respectively. 
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Figure S4. A plot of the surface tension (у;у) of the aqueous solutions of RH40 and ELP mixtures at the mole fraction of 
ELP in the bulk phase equal to 0.6 vs. the logarithm of their concentration (logC,,) at different temperatures equal to 
293 (a), 298 (b), 303 (c), 308 (d), 313 (e) and 318 K (f). Points 1 correspond to the measured values, curves 2 - 6 
correspond to the values calculated from the Szyszkowski equation (Eq. (2)), exponential function of the second order 
(Eq. (1)), Fainerman and Miller equation (Eq. (5)), Eq. (7) and Eq. (4), respectively. 
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Figure S5. A plot of the surface tension (у;у) of the aqueous solutions of RH40 and ELP mixtures at the mole fraction of 
ELP in the bulk phase equal to 0.8 vs. the logarithm of their concentration (logC,,) at different temperatures equal to 
293 (a), 298 (b), 303 (c), 308 (d), 313 (e) and 318 K (f). Points 1 correspond to the measured values, curves 2 - 6 
correspond to the values calculated from the Szyszkowski equation (Eq. (2)), exponential function of the second order 
(Eq. (1)), Fainerman and Miller equation (Eq. (5)), Eq. (7) and Eq. (4), respectively. 
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Figure S6. A plot of the surface tension (y;y) of the aqueous solutions of ELP vs. the logarithm of its concentration 
(logC) at different temperatures equal to 293 (a), 298 (b), 303 (c), 308 (d), 313 (e) and 318 K (f). Points 1 correspond 
to the measured values, curves 2 and 3 correspond to the values calculated from the Szyszkowski equation (Eq. (2)) and 


the exponential function of the second order (Eq. (1)), respectively. 
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Figure 57. A plot of the values of constant yọ in Eq. (1) (curve 1) and the minimal surface tension of aqueous solution 


(vi) (curve 2) vs. the temperature (T) for the RH40 and ELP mixtures at the mole fraction of ELP in the bulk phase 
equal to 0 (RH40 (a)), 0.2 (b), 0.4 (c), 0.6 (d), 0.8 (e) and 1 (ELP (f)), respectively. 
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Figure S8. A plot of the values of constant y, in Eq. (1) for studied surfactant mixtures vs. the mole fraction of ELP in the 
bulk phase (æ) at the temperatures equal to 293 (curve 1), 298 (curve 2), 303 (curve 3), 308 (curve 4), 313 (curve 5) 
and 318 K (curve 6), respectively. 


1 —1— 293K 
2 —*— 298K 
3 —4— 303K 
4 —v— 308K 
5m$9—313K 
6 —4— 318K 


a EEP 


Figure S9. A plot of the values of constant A, in Eq. (1) for studied surfactant mixtures vs. the mole fraction of ELP in the 
bulk phase (æ) at the temperatures equal to 293 (curve 1), 298 (curve 2), 303 (curve 3), 308 (curve 4), 313 (curve 5) 
and 318 K (curve 6), respectively. 
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Figure S10. A plot of the values of constant A, in Eq. (1) for studied surfactant mixtures vs. the mole fraction of ELP in 
the bulk phase (a) at the temperatures equal to 293 (curve 1), 298 (curve 2), 303 (curve 3), 308 (curve 4), 313 (curve 
5) and 318 K (curve 6), respectively. 
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Figure S11. A plot of the values of constant t, in Eq. (1) for studied surfactant mixtures vs. the mole fraction of ELP in 
the bulk phase (a) at the temperatures equal to 293 (curve 1), 298 (curve 2), 303 (curve 3), 308 (curve 4), 313 (curve 
5) and 318 K (curve 6), respectively. 
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Figure S12. A plot of the values of constant t; in Eq. (1) for studied surfactant mixtures vs. the mole fraction of ELP in 
the bulk phase (a) at the temperatures equal to 293 (curve 1), 298 (curve 2), 303 (curve 3), 308 (curve 4), 313 (curve 
5) and 318 K (curve 6), respectively. 
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Figure 513. The values оГГ” calculated from Eq. (6) and (2) for studied surfactant mixtures at the mole fraction of 


ELP in the bulk phase (æ) equal to 0, 0.2, 0.4, 0.6, 0.8 and 1 at different temperatures equal to 293 (a), 298 (b), 303 (c), 
308 (d), 313 (e) and 318 K (f), respectively. 
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Figure 514. A plot of the values of Gibbs surface excess concentration (T) vs. logarithm of the concentration in the bulk 
phase (Cor (12) for the RH40 and ELP mixtures at the mole fraction of ELP in the bulk phase equal to 0 (RH40 (a)), 0.2 
(b), 0.4 (c), 0.6 (d), 0.8 (e) and 1 (ELP (f)), respectively. 
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Figure S15. The values of the mole fraction of ELP in the mixed monolayer at the water air interface (X;) calculated 


from the relationship: X? = m (bars 1) and from Eq. (8) (bars 2) and temperature range 293-318K for mixtures at 


the mole fraction of ELP in the bulk phase (a) equal to 0.2 (a), 0.4 (b), 0.6 (c) and 0.8 (d). 
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Figure S16. A plot of the values of parameter of intermolecular interactions, 3°, for studied binary surfactant mixtures 
a) at the mole fraction of ELP in the bulk phase a, equal to 0.2, 0.4, 0.6, 0.8 and T = 293K vs. surface tension, (y,y), b) at 
the surface tension, (y; y), equal to 65, 55, 45 mN/m and T = 293K vs. mole fraction of ELP in the bulk phase a, c) at the 


mole fraction of ELP in the bulk phase a; equal to 0.2, 0.4, 0.6, 0.8 and y,y = 65 mN/m vs. temperature, T, d) at the mole 


fraction of ELP in the bulk phase a, equal to 0.6 and y; = 65, 55, 45 mN/m vs. temperature, T. 
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Abstract: The poor solubility of berberine (Ber) in water limits its practical use. Its solubility can be 
increased, among other ways, by the addition of surfactants. Of the surfactants, Kolliphor® ELP (ELP) 
and Kolliphor® RH 40 (RH40) can be very useful in this respect. The increase of Ber’s solubility in 
water in the presence of ELP and RH40 should be reflected in the composition of the surface layers at 
the water-air interface and the micelles. The determined composition is reflected in the Gibbs energy 
of interactions of berberine with ELP and RH40 through the water phase and the standard Gibbs free 
energy, enthalpy, and entropy of adsorption and micellization. These energies were determined from 
the equations proposed by us, based on the Gibbs surface excess concentration of the Ber mixture 
with ELP and КН40, the activity of these compounds in the surface layer at the water-air interface and 
in the micelles obtained by the Hua and Rosen method, and the contributions of Ber, ELP, and RH40 
to the reduction in the water surface tension. For this determination, the measurements of the surface 
tension of the aqueous solution of the Ber mixture with ELP or RH40 and that of the Ber mixture 
with these two surfactants, as well as the density and conductivity were performed. Moreover, the 
fluorescence emission spectra for the Ber + surfactant mixtures were recorded. 


Keywords: berberine; Kolliphor® ELP; Kolliphor® RH 40; adsorption; micellization; thermodynamic 
parameters of adsorption and micellization 


1. Introduction 


Berberine (Ber), which occurs in the roots, rhizomes, stems and bark of medicinal 
plants is widely used in medicine [1—10]. Its molecule consists of condensed hydrocarbon 
rings, a ring containing oxygen atoms, and a ring with a quaternary nitrogen atom, as 
well as the -CH3O groups (Figure 1). The predominance of hydrophobic groups over 
hydrophilic ones in the berberine molecule causes its poor solubility in water [11,12], which, 
among others, is the reason for the limited use of berberine in practice. The solubility of 
Ber in water can be increased, among other ways, by mixing it with suitable surfactants. 
Such surfactants include Kolliphor® ELP (ELP) and Kolliphor® RH 40 (RH40). These 
surfactants, which are well soluble in water, are nontoxic, allowing them to be used, for 
example, in cosmetics and pharmaceutical products [13-16]. The good solubility of ELP 
and RH40 results from the presence of many oxyethylene groups, as well as -C=O and 
=C-O groups in their molecules. However, in the ELP and RH40 molecules, hydrocarbon 
chains containing one -OH group are also present. This group decreases the hydrophobic 
properties of these chains to some extent. 

The type and amount of the chemical groups present in the molecule determine the 
values of the compound surface tension and its components, as well as the number of the 
water molecules that can be contacted with this compound in aqueous media. According to 
van Oss and Constanzo [17], the surface tension of surfactants depends on their molecules’ 
orientation toward the air phase. If they are oriented toward this phase by the hydrophobic 
parts of their molecules, then the surfactant surface tension is treated as its tail surface 
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tension. Otherwise, when the surfactant molecules are oriented by their hydrophilic 
part toward the air phase, then the surface tension is called the head surface tension. In 
turn, van Oss et al. divided the surface tension of all substances into the Lifshitz-van 
der Waals (LW) and acid-base (AB) components. The AB component associated with 
hydrogen bond formation is a function of electron-acceptor and electron-donor parameters 
of the surface tension [17—20]. For the condensed phases, van Oss et al. proved that 
the LW component results from the dispersion intermolecular interactions because the 
contributions of dipole-dipole and induced dipole-dipole interactions to this component 
are smaller than 2% [18-20]. 
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Figure 1. A plot of the surface tension (угу) of berberine (Ber) aqueous solutions vs. the logarithm of 
Ber concentration (log Свег) at a constant temperature equal to 293 K (curve 1), 308 K (curve 2), and 
313 K (curve 3). 


Knowledge of the values of the components and parameters of the surfactants' sur- 
face tension allows for the determination of the Gibbs free energy of interactions of the 
surfactants through the water phase. Additionally, knowing the contactable area of the 
hydrophobic and hydrophilic parts of the surfactants' molecules, it is possible to predict 
their tendency to adsorb and aggregate [21,22]. In the other words, based on the compo- 
nents and parameters of the surfactants' tail and head surface tension, as well as their 
contactable area, it is possible to predict the standard Gibbs free energy of adsorption 
and micellization. This process is possible for single surfactants. In the case of surfactant 
mixtures or their mixture with additives, it is more complicated. Moreover, such considera- 
tions are rarely found in the literature. In the case of such mixtures, the standard Gibbs 
energy of adsorption at the water-air interface and aggregation depends on the activity of 
all mixture components in the surface layer and in the micelles, which has not often been 
considered by investigators [23,24]. In fact, the activity depends on the composition and 
activity coefficient. 

The knowledge of the interactions between Ber and ELP, as well as RH40 molecules in 
aqueous media associated with the standard Gibbs free energy of adsorption and micel- 
lization, should enable understanding of the increase in Ber solubility in water, which is 
important from the practical point of view. As it is difficult to find such studies, the aim of 
our paper was to determine the relationship among the intermolecular interactions of the 
Ber, ELP, and RH40 through the water phase and the relationship between the composition 
of the mixed surface layer and mixed micelles, as well as the standard Gibbs free energy of 
adsorption and micellization. To solve this problem, measurements of the surface tension, 
conductivity, and density at temperatures of 293, 303 and 313 K of the following solutions 
were performed: ELP + Ber (M1), RH40 + Ber (M2), and ELP + RH40 + Ber (M3). In each 
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solution, the concentration of Ber was constant and equal to 1 х 1074 mol/dm*, and the 
concentration of surfactants was in the range of 1 x 1076-1 x 1072 mol/dm?. In the МЗ 
mixture, the mole fraction of ELP in the bulk phase was equal to 0.8. Moreover, the fluores- 
cence emission spectra for the Ber + surfactant mixtures were recorded. Additionally, for the 
determination of the components and parameters of berberine's surface tension, the contact 
angles of water, diiodomethane, and formamide on pressed Ber tablets were measured. 


2. Theory 


The amount of the surfactants adsorbed at the water-air (W-A) interface as well as 
the standard Gibbs free energy of adsorption and micellization, is directly associated with 
the chemical potential of the surfactants in the surface layer and/or micelles, as well as in 
the bulk phase. It should be remembered that, at a constant temperature (7) and pressure 
(p), the chemical potential of i-th solution component (1) is equal to the contribution of 
the Gibbs free energy (G;) of this component. The adsorption and micellization processes 
of a single surfactant are related to the standard Gibbs free energy of these processes. In 
turn, the standard Gibbs free energy of adsorption (АСО) and micellization (AGV, с) is 
equal to the difference between the standard potential, which is symmetrically (иб) апа 
asymmetrically (u;) defined. In the case of solutions of multicomponent mixture, the 
changes of the Gibbs free energy of the solution during the adsorption and micellization 
processes can be also related to the Gibbs free energy of surfactants mixing in the monolayer 
at the W-A interface and in the micelles. 

In the case of aqueous solutions, the amount of surfactant mixture adsorbed at the 
W-A interface is connected to the changes in the solution surface tension (yzy), which is 
expressed by the Gibbs-Duhem equation [25,26]: 


i=1 
Adyry + a nidu; = 0, (1) 


where n; is the number of moles of the i-th component, A is the area of the interface, and j 
is the number of the components. 

Knowing that, in the equilibrium state, the chemical potential of i-th component in 
the surface layer is equal to that in the bulk phase, assuming that the activity of j £ i 
components is constant and that the Gibbs surface excess concentration of i-th component, 
T;, is equal to 4, the following can be written: 


r, =i (9mv — 1 (dyw 0) 
! тКТ \ da; JT mRT \ dIna; md 


where the value of m depends on the type of surfactant; for example, for the ionic surfactant 
type 1:1 electrolyte is equal to 2, a; is the activity, and R is the gas constant. 

If the concentration of surfactants is small, then according to the asymmetrical defini- 
tion of the chemical potential, its activity coefficient (f;) is close to unity, and a; = х; (х; is 
the molar fraction of i-th component in the bulk phase). In such cases, x; = C;/w, where Cj 
is the mole concentration of the i-th component in the bulk phase, and w is the number of 
the water moles in 1 dm? of solution. Taking these assumptions into account, Equation (2) 
can be written as: 


r= C; OYLV = 1 дүгү (3) 
1 mRT \ OC; Cast 2.303mRT \ 0 log C; булт 


In the case of the constant composition of the multicomponent mixtures of substrates 
in the bulk phase, it was proved that the sum of the mixture’s Gibbs surface excess concen- 
tration at the W-A interface (Tsum) can be determined from the following equation: 


re Oum dyw Y — 1 дүгү (4) 
цэн ТАТ ОС r 203тШАТ \ dlog Csum / r 
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where Cs, is the sum of the mole concentration of all components of substrate in the solution. 

Unfortunately, in the case of the mixtures in which the concentrations of all compo- 
nents except one in the bulk phase are constant, the Gibbs surface excess concentration 
at the W-A interface of the components at a constant concentration in the bulk phase 
changes as a function of component at variable concentration. In such cases, it is possible 
to determine the Gibbs surface concentration using the Frumkin equation, which for i-th 
component of the mixture has this form: 


Т, 
Ti = – КТГ!" In (1 - pe) (5) 
1 


where 7t; is the contribution of i-th component to the reduction in water surface tension 
(Tw), Г; іѕ the concentration of i-th component of the mixture at the W-A interface, and I7*** 
is the maximal concentration of i-th mixture component at the W-A interface. 

It is proved that the fraction of area occupied by the i-th component in the mixture 


at the W-A interface (x?) is equal to where n j is the film pressure at the W-A 


ni" 
xis gn 

i=j i 
interface of the individual component of the mixture in the absence of others. Hence 7t; can 
be determined from the following expression [21]: 


7t = зал = Xi ? (Tw — yiv), (6) 


The studies of the composition of binary and ternary mixtures of surfactants have 
shown that the x? values determined based on the contributions of particular mixture 
components to the reduction in water surface tension are comparable to those determined 
from the Hua and Rosen concept [27]. The Hua and Rosen concept, based on the non- 
ideal solution theory, allows for determining the composition of the binary mixture of the 
surfactants in the interface region [27]. However, our studies showed that it is possible to 
apply the Hua and Rosen concept successfully to determine the composition of the mixed 
monolayer of the ternary surfactants’ mixture [21]. For the i-th, k-th, and l-th components 
of the ternary mixture, the Rosen and Hua equation can be written in this form: 


(5) In(a кСы/ХдС1) 21 (7) 
(oai) mA- aC 60] 
and 
Gd In (a Cj хиси) =i; (8) 
d= xi) In{(1 — ad)Cjy/ (1 — x) CP] 
and a SY In (a C9,/x$C9 
иСи/хдС4) E (9) 


(1 — 3$) In[( — a)C9,/ (1 — x$) C] 

where « is the mole fraction of a given component in the surfactant mixture in the bulk 
phase of the solution, C? is the concentration of the individual component of the mixture or 
the sum of components' concentrations at which the water surface ION is шэн to 
the 2. and ар = a; + Oy, р = Ri + Oy, Oy = Ry + oy, XS = х + xP, х5 = xP +27, 
x euwqe 0 - C? +C? + CO, Co = = со C9, CO = C? +С? and Со, = + ce 

The composition of the mixed micelles of binary mixtures of the ӨРЖЛЭНГ in the 
aqueous media can be determined using the Rubingh equation [28,29]. This equation can 
be successfully used for the ternary surfactant mixture in a modified way [21]. For the 
ternary mixtures, it can be written as: 


CR n(a Wining / Xin Cię) 
(1—xM)*In{(1 ПСИ — хМ)сМ] 


= 1, (10) 
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(xM) In(au C /xMCM) 21 (11) 
2 r 
(1-3 ) n[(— aa)Czj/ (1 — x) Cr] 
and 
M2 мем 
(xi ) In n(a aChy x; C; 1) -1 (12) 


(1 хм) In[(1 — a)CM/ (1 — xM)cM] 


where M refers to the micelles, and the other symbols have the same meaning as above. 

It is known that, for the ideal mixtures, there are direct relationships between the 
standard Gibbs free energy of adsorption and the mole fraction of components in the 
mixed monolayer and between the standard Gibbs free energy of micellization and the 
mole fraction of components of the mixed micelle. In the case of the non-ideal mixtures of 
surfactants, the activity of particular components of the mixture in the surface region and 
micelle should be known. The literature reports different methods for the determination 
of standard Gibbs free energy of adsorption (AG) ds) and micellization (АСО, с) [29]. These 
methods are useful for aqueous solutions of single surfactants. However, they are not quite 
suitable for aqueous solutions of the multicomponent surfactant mixtures [21]. 

In the aqueous solution of the multicomponent surfactant mixture, three regions can 
be distinguished: interface (5), bulk (B), and micellar (M). These regions can be treated as 
separate phases. The chemical potential of the i-th component of the surfactant mixture in 
these phases can be expressed as follows [30]: 


и? = ub + RT Ina? + 7;00;, (13) 

and 
uM = y + RTInaM, (14) 

and 
u? = pi? + RT Ina?, (15) 


In the equilibrium state, the chemical potential of i-th component is the same. Hence: 


B 


a: 
л” = ję * = RT In S — iOi, (16) 
i 
and > 
„OM — y = RTIn a (17) 


i 
At constant T and p, юэ m p. = = AG; and вэ" — už? = AG jci: In the case of 
surfactants, their concentrations in aqueous solution are small, and it can be assumed that 


ав = x? and x? = C;/w. Based on these assumptions the following can be written: 
Ci 
AG as i ; = RTIn SBE 7j Qj, (18) 
and : 
CMCa; 
AG) ic; = RTIn 5 19 
mic,i ig MEM 2 


where ©; is the area occupied by one mole of surfactant, CMC is the critical micelle 
concentration, and f; is the activity coefficient. 

According to the thermodynamic rule, the standard Gibbs free energy of adsorption or 
micellization of the i-th component of the surfactant mixture (AG?) satisfies the following 
equation [25,29]: 

AG? = AH? — TAS), (20) 
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where AH? is the standard enthalpy of adsorption or micellization, and AS? is the standard 
entropy of adsorption or micellization. 

Knowing the values of AG? at different temperatures and assuming that, in the given 
temperature range, AH? is constant, it is possible to determine As? from the following 
equation [25,29]: 

d(AG°) 
dT 


Indeed, knowing the AG? and AS? values, the AH? values can be calculated from 
Equation (20). The activity coefficients of the i-th component of the surfactant mixture 
needed for the standard Gibbs free energy of adsorption and micellization determination 
can be obtained, among other ways, on the basis of intermolecular interactions parameter 
for the mixed monolayer (f^) and the micelles (ВМ) from the Rosen and Hua equations, 
which for the two-component mixtures of surfactants have the following forms [27,29]: 


=—A9, (21) 


In ff eg oam (22) 
and 
In f? = B^ (3y, (23) 
and 2 
In £M = gM(1—3M), Q4) 
and 5 
юу =p” (4 ). 05) 


For calculations of the B^ and BM parameters, Rosen and Hua derived the following 
equations [27,29]: 


B = In(aiCh, /x1 C0) (26) 
2 
(1 — xi) 
and ГЭРЭЭ” 
In (1C: / x^ C 
pM = (41С/5 rer) (27) 


(1 — xy”) 


Janczuk et al. [31,32] showed that there is a relation between the Gibbs free energy 
of adsorption and/or micellization and the components and parameters of the surface 
tension of water, as well as the tail and head of surfactant surface tension. They proposed 
the following equation: 


AGO, = (yr — Yrw)ATN + (тині — YwH) AHN, (28) 


where “үт is the tail surface tension; “түү, YwH, and Түүнд are the tail-water, head-water 
and dehydrated head-water interface tensions, respectively; Ar and Ан are the contactable 
areas of the tail and head; N is the Avogadro number. If during the transfer of the surfactant 
molecule from the bulk phase to the surface region, its head does not dehydrate, then: 


АС = (yr — yrw) ArN, (29) 


0 


The proposed equation for AG; ;. has the following form: 


AG), = —N |27wrAr = (2ywn - АС )3 н], (30) 


where AGF" is the Gibbs free energy of electrostatic interactions. 


For nonionic surfactants, Equation (30) assumes the following form: 


АСО ie = —N(2YwrAr — 2Ywu An), (31) 


mic 


Molecules 2023, 28, 3115 


7 of 24 


3. Results and Discussion 
3.1. Some Physicochemical Properties of Ber, ELP and RH40 


The behavior of surfactants and their mixtures with organic additives in aqueous me- 
dia depends not only on the types of the chemical groups present in the surfactant molecules 
and additives but also on their number, as well as their polarity. The sizes of the contactable 
area of the apolar and polar parts of the surfactants and additives with other molecules 
and their surface tension components and parameters are responsible for the adsorption 
and aggregation behaviors of surfactants and their mixtures with organic additives. 

The sizes of Ber, ELP, and RH40 molecules were calculated based on the bonds’ lengths, 
the angles between the bonds, and the distance between the molecules. It was proved 
that the volume of one molecule each of Ber, ELP, and RH40 can be deduced based on 
the volumes of cubes into which particular fragments of the molecule can enter. The 
volumes of Ber, ELP, and RH40 molecules obtained in this way are equal 428.25, 4378.64, 
and 4825.21 АЗ, respectively. As follows from these volumes and the Ber, ELP, and RH40 
mole weights, the density of these compounds is equal to 1.3042, 0.91947, and 0.9979 g/cm, 
respectively. Taking into account the molecule size of a given compound obtained in 
this way, it was possible to establish the contactable area of this molecule with those of 
water. This area is close to 277.45 A? for Ber, 3779.9 A? for ELP, and 4210.2 A? for RH40. 
Considering the contactable area of water molecules, which is equal to 10 A? [33], it can be 
stated that there can theoretically be contacted 28 water molecules with Ber, 378 with ELP, 
and 421 with КН40. In turn, one oxyethylene group in the surfactant molecule can be joined 
with two water molecules by strong hydrogen bonds and three by the weak ones [34,35]. 
Hence, it can be concluded that, depending on the configuration, the oxyethylene chains in 
the RH40 molecule can be surrounded at most by about 280 water molecules and ELP by 
245. As the contactable area of ELP and RH40 tails is twice as small as the head of their 
molecules, approximately two times fewer water molecules surround the tails of ELP and 
RH40 molecules than the heads. The numbers of water molecules surrounding the heads 
of surfactant molecules make them largely soluble in water. On the other hand, the number 
of water molecules surrounding the tails of RH40 and ELP molecules is the driving force in 
their adsorption at the W-A interface and micellization. This driving force depends on the 
components and parameters of water and surfactants’ surface tension, as well as on the 
size of the contactable area. The contactable area of the surfactant molecules through the 
water phase is smaller than that of the whole one. These areas for ELP and RH40 are equal 
to 951.616 and 1054.582 А?, respectively. In the case of Ber, its contactable area between 
two molecules is equal to 108.96 À?. This area is smaller than the maximal contactable area 
between the tails of two ELP and RH40 molecules, which is close to 324.714 A2. 

The surface tension of Ber was determined based on the contact angle (0) measured 
for water (56.2°), formamide (38.4°), and diiodomethane (46.1°) on pressed berberine using 
the van Oss et al. concept [18-20]. According to this concept and the Young equation, the 
following can be written: 


yLv(cos@ + 1) = 2 (v W + ty ty "es ). (32) 


Knowing the Lifshitz-van der Waals component (7) of the water, formamide and 
diiodomethane surface tension, and the electron-acceptor (р) and electron-donor (yry) 
parameters of this tension [36], it was possible to obtain from Equation (32) the Lifshitz-van 
der Waals component (927 = 36.42 mN/m) of the Ber surface tension (ysy = 46.52 mN/m), 
as wellas its electron-acceptor ay = 1.56 mN/m) and electron-donor (усу = 16.38 mN/m) 
parameters. From the calculations based on Equation (32), it resulted that the Ber surface 
tension is not much different from that of Triton X-165 (TX165) if its molecules are positioned 
with the hydrophilic part toward the air phase [37]. However, the contribution of the 
LW component to the Ber surface tension is smaller, and that of AB larger than those 
of TX165 surface tension. The contribution of the LW component to Ber ys is equal 
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to 78%, and it results in the poor solubility of Ber in water. Indeed, the solubility of 
organic substrates in water depends not only on the LWtoAB ratio and their surface 
tension but also on the contactable area of the hydrophobic and hydrophilic groups in 
their molecules. To increase the solubility of Ber in water, it is added to water in the 
form of an alcohol solution. Probably due to strong interactions between the hydrocarbon 
groups in the alcohol and in the Ber molecules, the displacement of water molecules 
surrounding the hydrophobic groups in the Ber molecule, as well as the orientation of the 
alcohol molecules with the -OH group toward the water phase, enhances the hydrophilic 
character of the Ber molecules. The constant concentration of ethanol used to dissolve Ber 
in water, equal to 1 x 1073 mol/dm?, practically changes the water surface tension by a 
value comparable to the accuracy of surface tension measurements. It seems that, in such 
cases, all thermodynamic considerations dealing with the adsorption and micellization 
properties of Ber, ELP, and RH40, as well as their mixtures, are reasonable. 

In the Ber molecule, in contrast to ELP and КН40, it is difficult to distinguish the 
hydrophobic (tail) and hydrophilic (head) parts. Therefore, it should be expected that 
the preferential orientation of Ber molecules at the W-A interface is parallel to this inter- 
face. On the other hand, it is possible that the hydrophilic group in the Ber molecule 
causes its orientation at the W-A interface in the form of an inclined plane. A greater 
tendency of hydrophobic chains composed of -CH3, -CH»- or =CH- groups for horizontal 
orientation than toward reorientation is suggested by the investigators on the basis of ther- 
modynamic considerations of the adsorption processes of different compounds at the W-A 
interface [38—40]. Indeed, this orientation depends largely on the packing of the compound 
molecules being influenced by the strong intermolecular interactions in the surface region. 

Berberine reduces the water surface tension to a small extent (Figure 1), and its 
maximal reduction corresponds to a Ber concentration close to 1 х 1074 mol/dm?. At 
higher Ber concentrations in aqueous solutions, the surface tension of solution is constant. 
However, the inflection point on the surface tension isotherm of the berberine aqueous 
solution cannot be treated as the CMC. The Gibbs surface excess concentration of Ber at 
the W-A interface calculated from Equation (4) is equal to 1.5 х 1075, 1.48 х 1075 and 
1.46 x 1076 mol/dm? at 293, 303, and 313 K, respectively. The minimal area occupied by 
one Ber molecule is close to the contactable area of berberine at its parallel orientation at 
the W-A interface. In contrast to berberine, the ELP and RH40 molecules in the surface 
layer can be oriented perpendicularly and/or at an angle to the W-A interface because their 
molecules are of an amphiphilic nature. 

It seems that, in the case of ELP and RH40 the directly measured surface tension can 
be treated as the head surface tension because of very strong hydrophobic interactions 
between the tails of ELP and RH40 molecules. For this reason, the probability of the 
orientation of ELP and RH40 molecules at the surfactant-air interface with their heads 
directed toward the air phase is very high. If so, it is possible to determine the components 
and parameters of the heads of ELP and RH40 surface tension by measuring their contact 
angles on the PTFE and PMMA surfaces. It is known that PTFE is a modeling apolar solid, 
the surface tension of which results from only the Lifshitz-van der Waals intermolecular 
interactions [18-20]. According to Fowkes [41] and van Oss et al. [18-20], the equilibrium 
state of the PTFE-liquid drop-air system fulfills this equation: 


Yuv (cos 8 + 1) = TT (33) 


On the basis of ^77! and the contact angle on the monopolar PMMA, the electron- 
acceptor parameter of ELP and RH40 head surface tension can be determined. It should 
be noted that PMMA is treated as a monopolar solid because its surface tension results 
from the LW intermolecular interactions [32]. On the other hand, PMMA can also interact 
with polar liquids by acid-base (AB) forces because the y~ parameter of the PMMA surface 
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tension is greater than zero. The contact angle of liquids on the PMMA surface can be 
expressed by the following equation [18]: 


rev (cos 8+ 1) = 2( YAE + nsns (34) 


Knowing that yA? = «ry — YH, it is possible to determine all components and 
parameters of the ELP and RH40 head surface tension based on the measured values of yry 
and 0 on the PTFE and PMMA surfaces. These components and parameters determined in 
this way are similar to those of the head of TX165 surface tension (Table 1). This outcome 
confirms the assumption that ELP and RH40 molecules are oriented by their heads toward 
the air phase. In turn, as suggested earlier, the surface tension of the ELP and RH40 heads 
should be close to those of unsaturated fatty acids [22]. One -OH group in the hydrophobic 
single chain of ELP and RH40 molecules does not seem to have a significant effect on the 
value of the tail surface tension. The large size of the head of ELP and RH40 molecules 
guarantees good solubility of these surfactants in water, and the large size of the tail causes 
good surface activity and the tendency to form micelles in aqueous media. The maximal 
Gibbs concentration of ELP and RHAO is greater than for Tritons, and micelle formation 
takes place at concentrations smaller than that in the case of Tritons [42]. 


Table 1. The values of the contact angles of water (0w), formamide (0r), and diidomethane (0р) 
on the pressed Ber; the contact angles (0) of ELP and RH40 on PTFE and PMMA surfaces; and the 
components and parameters of the surface tension of water, Ber, ELP, ELP/RH40/TX165 head, PTFE 
and PMMA. 


0w 0r 0p 
Ber 56.2 38.4 46.1 
0 PTFE 0 PMMA 
ELP 74 23 
RH40 75 25 
TW Yiv TW YLV 
Ber 36.42 1.556 16.38 46.52 
ELP head 27.51 0.460 48.95 37.00 
КН40 head 27.38 0.565 44.42 37.40 
Water 26.85 22.975 22.975 72.80 
TX165 head 27.70 0.33 50.20 35.84 
PTFE 20.24 - - 20.24 
PMMA 41.28 - 7.28 41.28 


3.2. Surface Tension of the Aqueous Solution of ELP, RH40, and ELP + RH40 Mixtures with Ber; 
Concentration and Composition of Mixed Monolayers at the Water-Air Interface 

Taking into account the contribution of the Lifshitz-van der Waals intermolecular 
interactions in the surface tension of water, Ber, ELP, and RH40, it can be stated that 
reduction in water surface tension due to their adsorption at the W-A interface depend 
on only the decrease in the acid-base component of the water surface tension. The LW 
component of Ber, ELP, and RH40 surface tension is larger than that of water’s surface 
tension. In the case of the ELP and КН40 tail surface tension, the LW component is 
only insignificantly larger than that for water (Table 1) [36]. As mentioned above, in Ber 
molecules, it is difficult to distinguish the tail and head. Therefore, Ber molecules adsorbed 
in the monolayer at the W-A interface should be oriented parallel to the interface. Hence, 
the minimal surface tension of the aqueous solution of Ber cannot be smaller than its surface 
tension (46.52 mN/m) (Table 1). 
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Assuming that only the surface tension of the ELP and RH40 tail determines the 
reduction in water surface tension, it can be stated that the minimal surface tension of the 
aqueous solution should be considerably smaller than that of the Ber aqueous solution. 
From the differences between the surface tension of Ber and the tail of ELP and RH40, 
it results that, if the mixed monolayer is formed at the W-A interface from the solution 
including ELP, RH40, and ELP + RH40 with Ber, the surface tension of the aqueous solution 
of these mixtures should be greater than that for the solution in the absence of Ber at 
the same concentrations of ELP, RH40 and ELP + RH40. Removal of the ELP or RH40 
molecules from the mixed monolayer at the W-A interface by Ber molecules increases the 
solution surface tension. This suggestion is confirmed by the isotherms of the surface 
tension of the aqueous solution of ELP, RH40, and ELP + RH40 mixture in the presence and 
absence of Ber [22] (Figures 2—4). These isotherms have almost the same shape, and at a 
given temperature, there is almost a linear dependence between the surface tension of the 
solution and the logarithm of the surfactants or their mixture concentration. The isotherms 
of the surface tension of the aqueous solution of ELP, RH40, and ELP + RH40 mixture in 
the presence of Ber, similar to those in it absence [22], can be described by the exponential 
function of the second order (Figure S1): 


26 26 
Tw -y 22163 222163) 89) 


where A1, А», ty, to, and y? are constants. 


60 l 3 
55 |- 


50 H 


y, (mN/m) 


45 | 


40 - 


35 LL А 1 А 1 А 1 А 1 


Figure 2. A plot of the surface tension (yzy) of M1 aqueous solutions (curves 1, 1” and 1") and ELP 
aqueous solutions [22] (curves 2, 2' and 2") vs. the logarithm of ELP concentration (log Cg; p) at a constant 
temperature equal to 293 K (curves 1 and 2), 303 K (curves 1’ and 2^), and 313 K (curves 1” and 27). 


It was suggested earlier that the constants in Equation (35) are related to the compo- 
nents and parameters of the surface tension of all components of the solution, which are 
decisive for the interactions between these components [22]. It seems that the constant 
y? is related to the LW component of the surface tension of all compounds present in the 
solution and is close to the minimal surface tension of a given solution. This constant 
changes as a function of temperature almost linearly (Figure 52а). It is likely that the 
other constants in Equation (35) are closely related to the AB component of the surface 
tension of water and other compounds present in the aqueous solution. As a matter of 
fact, the changes in the hydration degree of the Ber, ELP, and RH40 molecules and the 
distance between these molecules and water can influence on the values of the Aj, A>, 
tı, and to constants. Presumably for these reasons, the changes in А], Аз, t1, and tę asa 
function of the temperature are more complicated than the y? constant (Figure $2). The LW 
components of the Ber, ELP, and RHAO tail surface tension as well as the degree of Ber, tail, 
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and head of ELP and RH40 molecule hydration, are decisive regarding the concentration 
and composition of the mixed monolayer at the W-A interface. 


60 T = 293K 
1—»—M2 
2—e—RH40 

ES T = 303K 
1'—>— М2 
= 2' —0—RH40 

Е T = 313K 

Ex 1" —r—M2 

= 2" —9— RH40 

= 


45 


40 


Figure 3. A plot of the surface tension (угу) of M2 mixture aqueous solutions (curves 1, 1’ and 1") 
and RH40 aqueous solutions [22] (curves 2, 2' and 2”) vs. the logarithm of RH40 concentration (log 


CRH40) at a constant temperature equal to 293 K (curves 1 and 2), 303 K (curves 1' and 2’), and 313 К 
(curves 1" and 2”). 


65 + 
L 4% T = 293K 
"D 1—e«— M3 
2 —x— ELP+RH40 
Г Т = 303K 
Be 1" —0— M3 
e i 2' —+— ELP+RH40 
5 5L T = 313K 
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Figure 4. A plot of the surface tension (угу) of M3 mixture aqueous solutions (curves 1, 1' and 1") and 
ELP + RH40 mixture aqueous solutions at a mole fraction of ELP equal to 0.8 [22] (curves 2, 2' and 
2") vs. the logarithm of the sum of ELP and RH40 concentration (log ((CgLp + Скнао)) at a constant 
temperature equal to 293 K (curves 1 and 2), 303 K (curves 1' and 2’), and 313 K (curves 1" and 2"). 


The concentration of a given component of the surfactant mixture and surfactant mixtures 
with additives can be indirectly determined, among other ways, using the Gibbs isotherm 
(Equation (4)) and Frumkin equation (Equation (5)). Unfortunately, in the studied concen- 
tration range, which corresponds to the saturated mixed monolayer at the W-A interface, it 
was possible to determine the maximal concentration of ELP, RH40, and ELP + RH40 mixture 
in the presence of Ber. It should be noted that the poor solubility of berberine in water will 
make it impossible to measure the surface tension of the aqueous solutions of the mixture 
of surfactant and Ber in a broad range of constant surfactant concentration. For this reason, 
the concentration of the particular surfactant in the single and mixed monolayers at the W-A 
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interface was determined using the Frumkin equation (Equation (5)). To solve this equation in 
the case of the mixtures, the contribution of their particular components to the reduction in 
the water surface tension must be known. It was found earlier that the composition of the 
mixed monolayer can be determined based on the yzy isotherm of the aqueous solution of 
the single components of the mixture [21,22]. If so, the reduction in the water surface tension 
by a given mixture component (7t;) can be calculated from Equation (6). 

In many cases, it is possible to determine the composition of the binary and ternary mixtures 
of surfactants using the modified Rosen and Rubingh equations (Equations (8)-(10)) [27,29]. It 
is known that these equations can be used for the determination of the mole fraction of 
the mixed monolayer only in limited mixture compositions in the bulk phase. For this 
reason, it was possible to determine the composition of the mixed monolayer only at small 
concentrations of surfactants in the М1-М3 mixtures. It was proved that the composition of 
the mixed monolayer established in this way is close to that determined based on the Try 
isotherm of the aqueous solutions of the individual mixture components. The mole fraction 
of particular components of the M1, M2, and M3 mixtures in the mixed monolayer at the 
W-A interface is considerably different from that in the bulk phase and depends on the 
temperature (Figures 53-55). It is interesting that, at large concentrations of ELP, RH40, or 
the mixture of ELP + RH40, the mole fraction of Ber in the mixed monolayer is larger than 
in the bulk phase. This finding indicates that the Ber molecules likely adsorb at the W-A 
interface together with the surfactant ones. In this case, as a result of strong hydrophobic 
interactions between Ber and surfactant molecules, the Ber molecules in the mixed mono- 
layer are not oriented parallel and toward the interface but perpendicularly or at a given 
angle to the interface. This orientation results in an increase in the Ber concentration in the 
mixed monolayer. In fact, if the berberine concentration in the mixed monolayer is larger 
than in the bulk phase, then the relation of the mole fraction of surfactants in these phase is 
reversed. The relation of the composition of the mixed monolayer, including ELP, RH40, 
and Ber, to that in the bulk phase is slightly different from that for the solutions of mixture 
M3. In the whole concentration range, the mole fraction of RHAO in the mixed monolayer is 
larger than in the bulk phase. This relation is consistent with the results for the ELP and 
RH40 mixture in the absence of Ber [22]. 

Taking into account the mole fraction of particular mixture components in the mixed 
monolayer at the W-A interface, the contribution of these components to the reduction 
in water surface tension was determined using Equation (6). Next, using the Frumkin 
equation, the concentration of the mixture components in the surface region was calculated. 
The concentration of the particular components in the mixed monolayer at the W-A interface 
for all studied mixtures is smaller than those of the individual components in the absence 
of others (Figures 56-58). Moreover, the sum of the concentrations of Ber and the surfactant 
or of Ber and the ELP + RH40 mixture in the mixed surface monolayer is smaller than that 
of the surfactant or surfactants mixture in the absence of Ber (Figures 56-58). What causes 
these phenomena? It is possible, as mentioned above, that the Ber molecules adsorb together 
with the surfactant ones. The berberine molecules can join the tail of surfactant molecules, 
decreasing the extent of hydration. As a consequence, a complex of surfactant-Ber is formed 
with a smaller tendency to adsorb at the W-A interface than the single molecule of ELP or 
RH40. On the other hand, the adsorbed berberine molecule + surfactant tail complex in 
the mixed monolayer can be directed toward the water phase, changing its surface tension. 
Due to the Ber molecule's structure, it has a positive charge; hence, repulsive electrostatic 
interactions between the Ber + surfactant tail complexes can appear. These interactions can 
reduce the packing of the mixed monolayer, decreasing the concentration of surfactant in 
comparison to its concentration in the absence of Ber. 


3.3. CMC of ELP + Ber, RH40 + Ber and ELP + RH40 + Ber Mixtures 


The ability of surfactants to form aggregates in aqueous solution is one of their char- 
acteristic and a very important property. The micellar phase plays a very important role, 
among others, because of its solubilization properties. We were interested in whether Ber is 
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present in the micelles formed by ELP, RH40, and their mixture. To consider this problem, 
it was first necessary to determine the CMC and the composition of the mixed micelles. 
In the literature, there are many methods for CMC determination, for example, based on 
the isotherms of the surface tension, density, and conductivity and using spectroscopic 
methods. Formation of micelles in aqueous media is reflected by the inflection point ap- 
pearing on these isotherms [29]. However, it should be remembered that each method 
can be sensitive to a different size and shape of aggregates formed by surfactants. The 
CMC values of surfactant + Ber mixtures (М1-М3) determined from the yzy isotherms 
(Figures 2—4), as well as conductivity (Figure 5), density (Figure 6) measurements, and 
fluorescence emission spectra (Figures 7 and 8 as an example), are different even for the 
same mixture (Table 2). The CMC values of the surfactants or their mixture in the presence 
of Ber are higher than those in its absence. 
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Figure 5. A plot of the specific conductivity (к) of mixtures M1 (a), M2 (b), and МЗ (с) vs. the 
logarithm of the surfactant or their mixture concentrations at a constant temperature equal to 293 K 
(curve 1), 303 K (curve 2) and 313 K (curve 3). 
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3—+—313K 


log (Corp * Care) 


Figure 6. A plot of the density (0) of mixtures M1 (a), M2 (b) and M3 (c) vs. the logarithm of the 
surfactant or their mixture concentration at a constant temperature equal to 293 K (curve 1), 303 K 
(curve 2), and 313 K (curve 3). 
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3—4— 313K 


log (Corp + Caro) 


Figure 7. A plot of the fluorescence intensity (I) of mixtures M1 (a), M2 (b), and M3 (c) vs. the 
logarithm of the surfactant or their mixture concentration at a constant temperature equal to 293 K 
(curve 1), 303 K (curve 2), and 313 K (curve 3). 
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Ix 10? 


2 (nm) 


Figure 8. Fluorescence emission spectra of mixture M1 at T = 293 K. 


Table 2. The values of the CMC (mol/ dm?) of the aqueous solutions M1, M2, and M3 mixtures as, 
well as ELP, RH40 and their mixture at T = 293, 303, and 313 K, determined from the surface tension 
isotherms (угу = f(log C)), conductivity (x = f(1/C)), and density (p = f(1/C)) measurements, as well 
as fluorescence emission spectra (I = f(1/C)). The values of CMC for ELP, RH40, and their mixture, 
determined from surface tension measurements, are taken from Ref. [22]. 


T(K) CMC from CMC from CMC from CMC from I 
YLV K р 
293 2.11 х 105 3.87 x 107% 5.92 x 10-4 
a 303 2.78 x 1075 4.02 x 1074 5.40 х 10-4 
313 2.72 x 107? 4.67 x 10-4 5.05 x 10-4 
293 2.14 x 1075 5.13 x 1074 6.03 x 10-4 
ELP 303 2.03 x 1075 7.78 x 1074 4.16 x 1074 
313 1.91 x 105 5.32 x 10-4 3.71 х 1074 
293 4.24 x 1075 7.31 x 1074 4.97 x 107% 
M 303 5.62 x 1075 7.70 x 1074 4.67 x 1074 
313 5.05 x 1075 4.46 x 1074 4.16 x 1074 
293 6.64 x 1075 6.03 x 107% 4.97 x 1074 
RH40 303 2.50 х 1075 5.77 х 107% 4.67 х 10-4 
313 2.12 x 1075 5.22 x 107% 4.42 x 10-4 
293 3.13 x 1075 3.81 x 1074 4.74 x 107% 
ELEA ав + Ber 303 2.81 x 10-5 4.23 x 1074 4.33 x 1074 
313 2.80 x 1075 4.74 x 107% 4.16 x 107% 
293 1.92 x 105 8.76 x 1074 4.89 x 1074 
ELP + RH40 303 1.84 x 1075 7.21 x 107% 4.40 x 10-4 
313 1.70 x 1075 7.03 x 107% 3.90 x 107% 


Formation of micelles in aqueous solution by surfactants at their given concentrations 
is due to hydrophobic interactions between the tails of surfactants through the water phase. 
The values of these interactions are positive, in contrast to the interactions of the surfactant 
head, which are negative independent of whether micelles are formed by nonionic or ionic 
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surfactants. According to Equations (30) and (31) [31,32], the power of the hydrophobic 
interactions depends on the tail-water interface tension and the contactable area of the 
tail. The interface tension can be calculated from, among others, the van Oss et al. [18—20] 
concept. From this concept it follows: 


Wyant -2( тре + af + sf) (36) 


If the surface tension of at least one of the phases being in contact results from only 
the LW intermolecular interactions, then Equation (36) takes the following form: 


qij = Vi + Yj = 2y УР", (37) 


Taking into account Equations (36) and (37) the components and parameters of berber- 
ine, water, and surfactant tail surface tension (Table 1), the values of the water tails of 
surfactants, water-Ber, and Ber-surfactant tail interface tension were calculated. The ob- 
tained values of the water-surfactant tail, water-Ber, and Ber-surfactant tail interface tension 
are close to 46.0, 6.0, and 10.7 mN/m, respectively. Based on these values, it was possible 
to determine the adhesion work of the surfactant tail to the surfactant tail, the surfactant 
tail to Ber, and Ber to the surfactant tail through the water phase. The adhesion work of the 
surfactant tail to the surfactant tail through the water phase is equal to 2ywr = 92 mJ/m? 
and that of Ber to Ber through the water phase to ywg = 12 mJ/ m^. In the case of the Ber- 
water-surfactant tail system, the adhesion work (WPT) was calculated from the following 
equation [25]: 

WE! = ywr + Yws — YBT, (38) 


The value of WPT calculated from Equation (38) is equal to 41.3 mJ/m?. 

As the surfactant head-water interface tension of both surfactants is negative and close 
to —18 mJ /m?, their total adhesion work through the water phase is close to 56 mJ /m?. 
Thus, there is not a large difference in the tendency to contact the surfactant with the 
surfactant molecule and the surfactant with the Ber molecule through the water phase. 
However, the probability of the binding of Ber+ surfactant complexes through the water 
phase is smaller than that of surfactant molecules. Presumably for this reason, the CMC of 
ELP, RH40, and ELP + RH40 mixtures is higher in the presence of Ber than in its absence 
(Table 2). Based on these facts, it can be assumed that the micellization process of the 
surfactant mixture occurs as a result of connecting not only surfactant molecules but also 
surfactants with Ber molecules. Thus, it seems that the presence of Ber in the micelles of 
ELP, RH40, and their mixture is not due to the adsorption of Ber molecules on the micelles 
and its penetration into the micelles but rather by common aggregation. For mixed micelles 
of ELP + Ber (M1) and RH40 + Ber (M2), it is difficult to determine the mole fraction of 
surfactant and Ber in the micelles based on their CMC values. This possibility was based 
on the modified concept of Hua and Rosen (Equations (10)-(12)) [27,29]. The calculated 
mole fractions of ELP, RH40, and Ber showed that the mole fraction of Ber does not differ 
significantly from 0.5. This outcome indicated that the interactions of berberine + surfactant 
complexes through the water phase play a major role in the aggregation process. Indeed, 
the mole fraction of particular components of ELP + RH40 + Ber mixtures (M3) changes to a 
small extent as a function of temperature (Table 3). The change in the mole fractions of the 
mixture components, as well as the CMC value itself as a function of temperature, results, 
on the one hand, from the change in kinetic energy and, on the other hand, from the change 
in hydration degree, especially of the surfactant heads. A change in the configuration of 
the surfactant molecules can also have an effect. 
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Table 3. The values of the fraction of area occupied by a component in the micelles in mixture M3 
(xM), intermolecular interactions parameters for mixed micelle (BM ), activity coefficient (f M), and 
the excess Gibbs energy of micelle formation per mole of the surfactant mixture (GM) calculated 
from the Rosen and Hua concept (Equations (10)-(12), (24), (25) and (27)) at T = 293, 303, and 313 K 
(1—ELP, 2—RH40, 3—Ber). 


T = 293 К Т = 303 K Т = 313 К 
CMC 3.13 x 105 2.81 x 105 2.80 x 1075 
Cy 2.50 x 105 2.25 x 1075 2.24 x 105 
Co 6.26 x 1075 5.63 x 1076 5.60 x 10-6 
C3 1 x 10-4 1 x 10-4 1 x 1074 
y CMC 40.27 39.17 37.92 
xM 0.3954 0.4568 0.4508 
xM 0.2313 0.0418 0.0024 
сон 0.3733 0.5014 0.5467 
(8)-2 —5.0834 0.3762 1.0993 
Ві) 1 —0.4159 — 1.5256 —1.1598 
“m 0.7619 1.0007 1.0007 
M 0.0496 1.4125 2.8469 
2 0.9371 0.7273 0.7900 
M 0.8590 0.6376 0.7048 
Gls) 2 -2.2018 0.0380 0.0682 
GO 1 —0.2422 —0.9536 —0.7472 


The influence of these factors on the composition and size of micelles should be 
reflected in the parameter of intermolecular interactions. This parameter can be determined, 
among other ways, from Equation (27) [27,29]. Since Ber does not form micelles on its own, 
it was not possible to determine the interactions parameter between Ber and surfactants 
molecules in the micelles formed in the ELP + Ber and RH40 + Ber systems (M1 and M2). 
Considering the complexes of ELP + Ber and RH40 + Ber as individual compounds, it 
was possible to determine the interactions parameter (BM) of Ber with ELP and RH40 in 
the mixed micelles of mixture M3 using the Hua and Rosen equation (Equation (27)). As 
follows from Table 3, the ВМ parameter, calculated based on CMC determined from the 
surface tension measurements, is negative for the RH40 + Ber complex at each temperature 
but for ELP + Ber complex only at T = 293 K. A similar situation can be observed for 
the values of the excess Gibbs energy of micelle formation per mole of the surfactant 
mixture (G™). In fact, the values of ВМ depend on the temperature for both complexes. 
The negative BM parameter suggests that there is synergism in the CMC. Unfortunately, 
this suggestion cannot be confirmed by the second condition for the existence of synergy 
for the abovementioned reason. It should also be noted that the values of the activity 
coefficient of ELP in the mixed micelles (Table 3) decrease, but for RH40, they increase with 
the T increase. 


3.4. Thermodynamic Parameters of the Adsorption and Micellization 


The standard Gibbs free energy, standard enthalpy, and entropy are useful to deter- 
mine the tendency of the surfactants to adsorb at different interfaces and to form micelles, 
as well as the reason for this tendency. In the literature, there are many different methods 
used to determine the standard Gibbs free energy of adsorption (AG®,.) and micellization 
(АСО, ), the standard enthalpy of adsorption A(HO аз) and micellization (АНО гс): and the 


mic 
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standard entropy of adsorption (AS0,.) and micellization (AS? ..) of the single compounds. 
Regarding cases in which the mixed monolayer at the interfaces and mixed micelles is 
formed, it is difficult to find in the literature methods strictly connected with thermody- 
namics rules. Therefore, the criteria that should be fulfilled to calculate the real values of 
the thermodynamic parameters of adsorption and micellization are presented above. Un- 
fortunately, it was impossible to determine the activity coefficients of the components in the 
mixed monolayer at the W-A interface based on the Hua and Rosen concept. Therefore, the 
AGO gs Values for Ber, ELP, and RH40 in the M1-M3 mixture was calculated from Equation 
(18), assuming the activity coefficients to be close to unity. The obtained results showed 
that, in the range of surfactants concentration at which they are present in the solution in 
the monomeric form, the AG? js Values are almost constant (Table 4). At the surfactants 
concentration higher than CMC the AG? js Values increase as a function of concentration. 
However, it should be mentioned that the surfactant molecules adsorb at the water-air 
interface only in the monomeric form. At surfactant concentrations higher than the CMC, 
their concentration in the monomeric form is constant. For this reason, the AGO di values 
are not real. 


Table 4. Thermodynamic parameters of the adsorption process of ELP, RH40 and Ber (kJ/mol). 


Т-293К Т-303К Т-313К 

AGqads АН ав TAS; AGgds AH ав TAS nds AG; АН ав TAS; 
ELP —46.91 —1.64 45.27 -47.97 -1.16 46.81 —50.00 —1.64 48.36 
ELP from M1 —44.62 —0.23 44.39 —46.16 —0.26 45.90 —47.65 —0.23 47.42 
ELP from M3 —43.64 2.65 46.29 -45.25 2.62 47.87 —46.80 2.65 49.45 
КН40 —48,12 4.62 52.74 —50.11 4.43 54.54 -51.72 4.62 56.34 
КН40 їгош М2 —45.69 2.54 48.23 —47.57 2.30 49.87 —48.98 2.54 51.52 
RH40 from M3 -44.45 —8.85 35.6 —45.72 —8.91 36.81 46.88 84.91 38.03 
Ber —32.00 —1.97 30.03 —33.10 —2.04 31.06 —34.05 —1.97 32.08 
Ber from M1 —31.63 —1.89 29.74 —32.62 —1.87 30.75 —43.66 —11.89 31.77 
Ber from M2 —31.42 -3.15 28.27 —32.32 —3.08 29.24 —33.45 -3.25 30.20 
Ber from M3 —31.02 -1.28 29.74 —32.02 =1.27 30.75 —32.95 —1.18 31.77 


To compare the AGO js Values obtained from Equation (18) to those determined using 
other methods, the calculations of AG?,. were performed using the Langmuir equation 
modified by de Boer [29,43]. This equation has the following form: 


А) A Ci AG isi 
A120 9 0 oe REN 2 


where A; and A9 are the areas occupied by one molecule of the i-th component of the 
mixture in the mixed monolayer and the limiting one. 

The AG? 4, Values calculated from Equation (39) based on the concentration of particu- 
lar components in the mixed monolayer at the W-A interface, determined from the Frumkin 
equation (Figure 59) in the range of surfactant concentrations smaller than CMC, are similar 
to those obtained from Equation (18) (Table 4, Figures 59-511). From the obtained AG? ds 
values, it can be stated that the Ber presence decreases the ELP and RH40 tendency to 
adsorb at the W-A interface. This outcome confirms our suggestion that the surfactant + Ber 
molecule complexes can be adsorbed at this interface. As mentioned above, there are some 
differences in the hydration degree of the molecule complexes in comparison to the individ- 
ual ones. In such cases, there is also a reduction in the difference between the tail-water and 
tail-air interface tension, which according to Equation (28) decreases AG... The standard 
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enthalpy of adsorption of particular components of the studied mixtures, determined based 
on Equations (20) and (21) in most cases, has a small absolute value except for the standard 
enthalpy value for RH40 adsorption at the W-A interface from the aqueous solution of 
mixture МЗ (Table 4). The small absolute values of ЛНО,, indicate that our suggestion 
about the adsorption of surfactant + berberine complexes is the most probable. 

In many cases, the absolute values of the standard enthalpy of micellization are greater 
than those of the standard enthalpy of adsorption (Tables 4 and 5). In most studied systems, 
the AH? гс is positive, indicating that some bonds are broken during the micellization 
process. In fact, it should refer to hydrogen bonds. There are some differences between 
the behavior of ELP and RH40. In solutions including RH40, the greatest changes of the 
standard enthalpy are observed in both adsorption and micellization processes compared 
to in the absence of RH40. 


Table 5. Thermodynamic parameters of the micellization process of ELP, RH40, and Ber (kJ/mol). 


T = 293 К T = 303 К Т = 313 К 

AG с АН? ic TAS с АС? с АН? ic TAS) с АС? i. АНО ic TAS? ic 
ELP —35.97 4.17 40.14 —37.33 4.18 41.51 —38.71 4.17 42.88 
ELPfromELP+ | 45g 6.68 42.49 87.33 6.61 43.94 -38.71 6.68 45.39 

RH40 

ELP from M3 —32.96 1.47 34.43 —33.96 1.64 35.60 —35.31 1.47 36.78 
RH40 23321 8.98 42.19 —36.80 6.83 43.63 —38.44 6.63 45.07 
RF RE AG —36.19 -8:23 32.96 —36.80 —2.71 34.09 —38.44 —3.23 35.21 
RH40 from M3 —31.28 16.92 48.20 39.88 89.72 49.84 —43.71 7.78 51.49 
Ber from M3 -29.15 4.55 33.70 —31.57 3.28 34.85 —32.72 3.28 36.00 


4. Materials and Methods 


Kolliphor® ELP (ELP) (Cremophor®ELP, Polyoxyl 35 hydrogenated castor oil, polyoxyl- 
35 castor oil), Kolliphor® RH 40 (RH40) (Cremophor® RH 40, macrogolglycerol hydroxys- 
tearate, PEG-40 castor oil, polyoxyl 40 hydrogenated castor oil) (Sigma-Aldrich (St. Louis, 
MO, USA) and berberine chloride (Ber) (Alfa Aesar, Kandel, Germany) were used without 
further purification. The doubly distilled and deionized water used for the preparation 
of the aqueous solutions of ELP + Ber (M1), RH40 + Ber (M2), and ELP + RH40 + Ber 
mixture (M3, the mole fraction of ELP (œ) in the bulk phase equal to 0.8) was obtained from 
a Destamat Bi18E distiller ([nkom Instruments, Warsaw, Poland). The surfactant solution 
concentration was from 1 х 1076 to 1 х 1072 mol/dm%, and the Ber concentration in the 
surfactant solutions was equal to 1 x 107^ mol/dm?. Because Ber is poorly soluble in 
water, it was introduced into the aqueous solutions in the form of an ethanolic solution, 
similar to other research [44-46]. The concentration of ethanol in the aqueous solution was 
constant and equal to 1 x 107? mol/dm?. 

The surface tension (угу) measurements of the aqueous solutions of the M1, M2, and 
M3 mixtures were performed at temperatures of 293, 303, and 313 K using a Krüss K100 
tensiometer (Krüss, Hamgurg, Germany), which was calibrated before the measurements, 
according to the platinum ring tensiometer method (du Nouy's method). The calibration 
was performed at 293 K using water and methanol, the surface tension values of which at 
this temperature were equal to 72.8 and 22.5 mN/m, respectively. The surface tension mea- 
surements for each concentration and composition of the studied solutions were repeated 
at least ten times. The standard deviation of the results obtained from the measurements 
was +0.1 mN/m, and the uncertainty was in the range of 0.3% to 0.9%. 

Measurements of the advancing contact angle (0) were performed using the ses- 
sile drop method with a DSA30 measuring system (Krüss, Germany) in a temperature- 
controlled chamber. For 0 measurements on pressed Ber water (Destamat Вї18Е), for- 
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mamide (>99.5%, Sigma-Aldrich (St. Louis, MO, USA) and diiodomethane (>99%, Sigma- 
Aldrich, St. Louis, MO, USA) were used. As previous studies showed, the value of the 
contact angle depends on the difference between the interface pressure and the hydrostatic 
pressure of the drop [47]. Therefore, in the contact angle measurements, the droplet sizes 
were different for diiodomethane, formamide, and water and equaled 4, 5, and 6 ст, 
respectively. Теп drops for each studied system were used, and the standard deviation was 
in the range of 1 to 1.5°. 

The conductivity measurements were performed using a Mettler Toledo™ Seven Multi 
with accuracy of +0.5%. 

The density of the studied solutions was measured with a U-tube densitometer (DMA 
5000 Anton Paar). The precision of the density and temperature measurements given by 
the manufacturer is +0.000005 g ст? and +0.001 K. The uncertainty was calculated to be 
0.01%. The densitometer was calibrated regularly with distilled and deionized water. 

Steady-state fluorescence measurements were performed using a Hitachi F-2700 Fluo- 
rescence spectrometer. Fluorescence excitation was determined at 450 nm, and the emission 
spectra were recorded in the range of 350-650 nm at a scan speed of 300 nm/min. The 
excitation and emission slit widths were 2.5 nm. 


5. Conclusions 


The measurements and discussion of the obtained results based on the thermodynamic 
rules allow for drawing interesting conclusions. Berberine (Ber) is a bipolar compound, 
and its surface tension results from the Lifshitz-van der Waals and acid-base components. 
However, the contribution of the Lifshitz-van der Waals component is considerably larger 
than that of the acid-base one, indicating that Ber has poor solubility in water. In turn, the 
electron-acceptor-parameter of the acid-base component is considerably smaller than that 
of the electron-donor one. 

Berberine reduces the water surface tension to a small extent, and its maximal Gibbs 
surface excess concentration is considerably smaller than those of ELP and RH40. 

The minimal area occupied by one Ber molecule is close to its contactable area at a 
parallel orientation toward the water-air interface. 

The components and parameters of the ELP and RH40 surface tension at the orientation 
of their molecules towards the air phase can be determined from their contact angle on the 
PTFE and PMMA surfaces. These components and parameters are similar to those of the 
TX165 surface tension at its orientation by the hydrophilic part towards the air phase. 

The molecules of ELP and RH40 adsorbed at the water-air interface reduced only the 
LW component of the water surface tension. 

The contribution of oxyethylene groups to the surface tension of ELP and RH40 is 
similar to that of Triton X-165. 

The surface tension isotherm of the aqueous solution of Ber mixture with ELP, RH40, 
and ELP + RH40 can be described by the exponential function of the second order. 

The LW components of the Ber, ELP, and RH40 tail surface tension, as well as the 
hydration degree of Ber molecules and the tail and head of ELP and RH40 molecules, 
are decisive regarding the concentration and composition of the mixed monolayer at the 
W-A interface. 

The composition of the mixed monolayer at the W-A interface, determined based on 
the surface tension isotherm of aqueous solutions of Ber, ELP, and КН40, is close to that 
obtained from the Rosen and Rubingh equations. This composition allows for determining 
the surface concentration of particular components in the mixed monolayer, as well as their 
concentrations in the mixed monolayer at the W-A interface using the Frumkin equation. 

At large concentrations of ELP, RH40, or the mixture of ELP + RHAO, the mole fraction 
of Ber in the mixed monolayer is larger than in the bulk phase. 

The concentration of the particular components in the mixed monolayer at the W-A 
interface for all studied mixtures is smaller than those of the individual components in the 
absence others. 
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The sum of the concentrations of Ber and the surfactant or of Ber and the ELP + RH40 
mixture in the mixed surface monolayer is smaller than that of the surfactant or surfactant 
mixture in the absence of Ber. 

The CMC values of the surfactant + Ber mixtures, determined from the surface tension 
isotherms, conductivity, density, and fluorescence emission spectroscopy, are different even 
for the same mixture. 

The concentration of Ber in the mixed micelles is higher than in the bulk phase. This 
outcome indicates that the tendency of Ber toward solubilization in the micelles of ELP, 
RH40, and ELP + RH40 is greater than its tendency to adsorb at the water-air interface. 

Using our thermodynamic considerations and the Hua and Rosen concept, it was pos- 
sible to determine the standard thermodynamic parameters of adsorption and aggregation. 


Supplementary Materials: The following supporting information can be downloaded at: https: 
/ / www.mdpi.com/article/10.3390/molecules28073115/51, Figure 51: A plot of the surface tension 
(ту) of M1 (a), M2 (b), and МЗ (c) aqueous solutions vs. the logarithms of ELP (log Cg; p) and RH40 
(log Сандо) and that of the sum of their concentrations (log(Cerp + Скндо)) at T = 293 K, 303 K, and 
313 K. Figure S2: A plot of the constant y (a), A1 (b), Aa (с), t1 (d), and tą (e) in Equation (35) vs. the 
temperature (T). Figure 53: A plot of the mole fraction of ELP and Ber in M1 in the surface layer 
and in the bulk phase vs. the logarithm of ELP concentration (log Cg; p) at T = 293 K, 303 K, and 
313 K. Figure 54: A plot of the mole fraction of RH40 and Ber in M2 in the surface layer and in 
the bulk phase vs. the logarithm of RH40 concentration (log Срнло) at T = 293 K, 303 K, and 313 К. 
Figure 55: A plot of the mole fraction of ELP, RH40, and Ber in M3 in the surface layer and in the bulk 
phase vs. the logarithm of the sum of ELP and RH40 concentrations (log(Cgr.p + CRH4o) at T = 293 К, 
303 K, and 313 K. Figure 56: A plot of the surface concentration (T) calculated from Equation (5) vs. 
the logarithm of ELP concentration (log Cg; p) at T = 293 K, 303 K, and 313 K. Figure 57: A plot of 
the surface concentration (T) calculated from Equation (5) vs. the logarithm of RH40 concentration 
(log Cgu40) at T = 293 K, 308 K, and 313 K. Figure 58: A plot of the surface concentration (Г) of the 
M3 mixture components calculated from Equation (5) vs. the logarithm of the sum of ELP and RH40 
concentrations (log(Cgr p + Сандо) at T = 293 К, 303 К, and 313 K. Figure 59: A plot of the Gibbs 
standard free energy of adsorption (AG? 4,) calculated from Equation (39) vs. the logarithm of ELP 
concentration (log Cg; p) at T = 293 K, 303 K, and 313 K. Figure 510: A plot of the Gibbs standard free 
energy of adsorption (AG? ds) calculated from Equation (39) vs. the logarithm of RH40 concentration 
(log CpH4o) at T = 293 K, 303 K, and 313 K. Figure 511: A plot of the Gibbs standard free energy of 
adsorption (AG? ds) Of the M3 mixture components calculated from Equation (39) vs. the logarithm of 
the sum of ELP and RH40 concentrations (log(Cgr p + Свндо)) at T = 293 K, 303 K and 313 K. 


Author Contributions: Conceptualization, M.S., K.S., A.Z. and B.J.; methodology, M.S. and K.S.; 
software, K.S.; validation, K.S., A.Z. and B.J.; formal analysis, K.S., A.Z. and В.Ј.; investigation, M.S. 
and К.5.; resources, K.S.; data curation, K.S., A.Z. and В.Ј.; writing—original draft preparation, K.S., 
A.Z. and B.J.; writing—review and editing, K.S., A.Z. and B.J.; visualization, K.S., A.Z. and BJJ.; 
supervision, В.Ј.; project administration, K.S., A.Z. and B.J.; funding acquisition, B.J. All authors have 
read and agreed to the published version of the manuscript. 


Funding: This research received no external funding. 
Institutional Review Board Statement: Not applicable. 
Informed Consent Statement: Not applicable. 


Data Availability Statement: The data presented in this study are available in Supplementary 
Materials. 


Conflicts of Interest: The authors declare no conflict of interest. 


Sample Availability: Samples of the compounds are not available from the authors. 


References 


1. Zhang, X; Dai, K.; Liu, C.; Hu, H.; Luo, E; Qi, Q.; Wang, L.; Ye, F.; Jin, J.; Tang, J.; et al. Berberine-coated biomimetic composite 
microspheres for simultaneously hemostatic and antibacterial performance. Polymers 2021, 13, 360. 
2. Patel, P. A bird’s eye view on a therapeutically ‘wonder molecule’: Berberine. Phytomed. Plus 2021, 1, 100070. [CrossRef] 


Molecules 2023, 28,3115 23 of 24 


10. 


11. 


12. 


13. 


14. 


15. 


16. 


17. 


18. 
19. 


20. 
21. 


22. 
23. 
24. 
25. 
26. 
27. 
28. 
29. 


30. 
31. 


Pang, B.; Zhao, L.H.; Zhou, Q.; Zhao, T.Y.; Wang, H.; Gu, C.J.; Tong, X.-L. Application of berberine on treating type 2 diabetes 
mellitus. Int. J. Endocrinol. 2015, 2015, 905749. [CrossRef] 

Xue, M.; Yang, M.-X.; Zhang, W.; Li, X.-M.; Gao, D.-H.; Ou, Z.-M.; Li, Z.-P.; Liu, S.-H.; Li, X.-J.; Yang, S.-Y. Characterization, 
pharmacokinetics, and hypoglycemic effect of berberine loaded solid lipid nanoparticles. Int. J. Nanomed. 2013, 8, 4677-4687. 
[CrossRef] [PubMed] 

Younis, F.A.; Saleh, S.R.; Abd El-Rahman, S.S.; Newairy, A.-S.A.; El-Demellawy, M.A.; Ghareeb, D.A. Preparation, physicochemical 
characterization, and bioactivity evaluation of berberine-entrapped albumin nanoparticles. Sci. Rep. 2022, 12, 17431. [CrossRef] 
Wang, L.; Ye, X.; Hua, Y.; Song, Y. Berberine alleviates adipose tissue fibrosis by inducing AMP-activated kinase signaling in 
high-fat diet-induced obese mice. Biomed. Pharmacother. 2018, 105, 121-129. [CrossRef] 

Lee, Y.S.; Kim, W.S.; Kim, K.H.; Yoon, M.J.; Cho, H.J.; Shen, Y.; Ye, J.-M.; Lee, C.H.; Oh, W.K.; Kim, C.H.; et al. Berberine, a natural 
plant product, activates AMP-activated protein kinase with beneficial metabolic effects in diabetic and insulin resistant states. 
Diabetes 2006, 55, 2256-2264. [CrossRef] 

Zhang, Y.; Li, X.; Zou, D.; Liu, W.; Yang, J.; Zhu, N.; Huo, L.; Wang, M.; Hong, J.; Wu, P.; et al. Treatment of type 2 diabetes and 
dyslipidemia with the natural plant alkaloid berberine. J. Clin. Endocrinol. Metab. 2008, 93, 2559-2565. [CrossRef] 

Pirillo, A.; Catapano, A.L. Berberine, a plant alkaloid with lipid- and glucoselowering properties: From in vitro evidence to 
clinical studies. Atherosclerosis 2015, 243, 449-461. [CrossRef] 

Wang, Y.; Huang, Y.; Lam, К.5.; Li, Y; Wong, W.T.; Үе, H.; Lau, C.W.; Vanhoutte, P.M.; Xu, A. Berberine prevents hyperglycemia- 
induced endothelial injury and enhances vasodilatation via adenosine monophosphate-activated protein kinase and endothelial 
nitric oxide synthase. Cardiovasc. Res. 2009, 82, 484—492. [CrossRef] 

Sahibzada, M.U.K.; Sadiq, A.; Faidah, H.S.; Khurram, M.; Amin, M.U.; Haseeb, A.; Kakar, M. Berberine nanoparticles with 
enhanced in vitro bioavailability: Characterization and antimicrobial activity. Drug Des. Dev. Ther. 2018, 12, 303-312. [CrossRef] 
Behl, T.; Singh, 5.; Sharma, N.; Zahoor, I.; Albarrati, A.; Albratty, M.; Meraya, A.M.; Najmi, A.; Bungau, S. Expatiating the 
pharmacological and nanotechnological aspects of the alkaloidal drugberberine: Current and future trends. Molecules 2022, 27, 
3705. [CrossRef] 

Guo, J.; Xia, Y.; Liu, Y.; Liu, S.; Zhang, L.; Li, B. Microscopic adsorption behaviors of ionic surfactants on lignite surface and its 
effect on the wettability of lignite: A simulation and experimental study. J. Mol. Lig. 2022, 345, 117851. [CrossRef] 

Corona, R.R.B.; Sad, C.M.S.; da Silva, M.; Lopes, D.L.; Leite, J.S.D.; Viegas, G.M.D.F.; Goncalves, G.R.; Filgueiras, P.R.; de Castro, 
E.V.R. Adsorption of anionic surfactant in graphite oxide: A study for treatment of laundry wastewater. J. Environ. Chem. Eng. 
2021, 9, 106858. [CrossRef] 

Abu-Huwaij, R.; Al-Assaf, S.F.; Hamed, R. Recent exploration of nanoemulsions for drugs and cosmeceuticals delivery. J. Cosmet. 
Dermatol. 2022, 21, 3729-3740. [CrossRef] 

Yu, E; Miao, Y.; Wang, M.; Liu, G.; Yuan, L.; Geng, R.; Qiu, Q.; Ni, C.; Kay, M. Predicting nanoemulsion formulation and studying 
the synergism mechanism between surfactant and cosurfactant: A combined computational and experimental approach. Int. J. 
Pharm. 2022, 615, 121473. [CrossRef] 

van Oss, C J.; Constanzo, P.M. Adhesion of anionic surfactants to polymer surfaces and low-energy materials. J. Adhes. Sci. 
Technol. 1992, 4, 477-487. [CrossRef] 

van Oss, C.J. Interfacial Forces in Aqueous Media; Marcel Dekker: New York, NY, USA, 1994. 

van Oss, С.Ј.; Good, R.J. Surface tension and the solubility of polymers and biopolymers: The role of polar and apolar interfacial 
free energies. J. Macromol. Sci. Chem. 1989, 26, 1183-1203. [CrossRef] 

van Oss, C.J.; Chaudhury, M.K.; Good, R.J. Monopolar surfaces. Adv. Coll. Interface Sci. 1987, 28, 35-64. [CrossRef] 

Szymczyk, K.; Lewandowski, A.; Zdziennicka, A.; Szaniawska, M.; Jańczuk, B. Behavior of Auramine O in the aqueous solution 
of two Kolliphors and their mixture. Molecules 2022, 27, 8493. [CrossRef] [PubMed] 

Szaniawska, M.; Szymczyk, K.; Zdziennicka, A.; Jańczuk, B. Adsorption properties and composition of binary Kolliphor mixtures 
at the water-air interface at different temperatures. Molecules 2022, 27, 877. [CrossRef] 

Markarian, S.A.; Sargsyan, H.R.; Grigoryan, G.S.; Ghazoyan, H.H.; Chaban, V.V. Dimethyl sulfoxide fosters larger micelles of 
docusate sodium in the mixed solutions with water. J. Mol. Liq. 2023, 369, 120960. [CrossRef] 

Lalthlengliani, J.; Gurung, J.; Pulikkal, A.K. Solubilization of aqueous-insoluble phenothiazine drug in TX-100 micellar solution 
and interactions of cationic/anionic surfactants with phenothiazine-TX-100 system. J. Mol. Lig. 2022, 354, 118823. [CrossRef] 
Adamson, W.; Gast, A.P. Physical Chemistry of Surfaces, 6th ed.; Wiley Interscience: New York, NY, USA, 1997. 

Atkins, P.; de Paula, J. Atkins’ Physical Chemistry; Eighth Edition; Oxford University Press: New York, NY, USA, 2006. 

Rosen, J.M.; Hua, X.Y. Surface concentrations and molecular interactions in binary mixtures of surfactants. J. Colloid Interface Sci. 
1982, 86, 164-172. [CrossRef] 

28. Rubingh, D.N. Solution Chemistry of Surfactants; Mittal, K.L., Ed.; Plenum Press: New York, NY, USA, 1979; Volume 3, pp. 
337-354. 

Rosen, M.J. Surfactants and Interfacial Phenomena, 3rd ed.; Wiley-Interscience: New York, NY, USA, 2004; pp. 34-178. 

Butler, J.A.V. The Thermodynamics of the Surfaces of the Solutions. Proc. Roy. Soc. Ser. A 1932, 138, 348-375. 

Jańczuk, B.; Méndez-Sierra, J.A.; González-Martin, M.L.; Bruque, J.M.; Wojcik, W. Properties of decylammonium chloride and 
cesium perfluorooctanoate at interfaces and standard free energy of their adsorption. J. Colloid Interface Sci. 1997, 192, 408-414. 
[CrossRef] [PubMed] 


Molecules 2023, 28,3115 24 of 24 


32. 


33. 
34. 


35. 


36. 


37. 


38. 


39. 


40. 


41. 


42. 


43. 


44. 


45. 


46. 


47. 


Zdziennicka, A.; Krawczyk, J.; Szymczyk, K.; Jańczuk, B. Macroscopic and microscopic properties of some surfactants and 
biosurfactants. Int. J. Mol. Sci. 2018, 19, 1934. [CrossRef] 

Groszek, A.J. Selective adsorption at graphite /hydrocarbon interfaces. Proc. R. Soc. Lond. A Math. Phys. Sci. 1970, 314, 473-478. 
Beyer, K. Phase structures, water binding, and molecular dynamics in liquid crystalline and frozen states of the system Triton 
X-100-D>O: A deuteron and carbon NMR study. J. Colloid Interface Sci. 1982, 86, 73-89. [CrossRef] 

Desai, T.R.; Dixit, S.G. Interaction and viscous properties of aqueous solutions of mixed cationic and nonionic surfactants. 
J. Colloid Interface Sci. 1996, 177, 471-477. [CrossRef] 

Zdziennicka, A.; Szymczyk, K.; Krawczyk, J.; Jańczuk, B. Some remarks on the solid surface tension determination from contact 
angle measurements. Appl. Surf. Sci. 2017, 405, 88-101. [CrossRef] 

Szymczyk, K.; Zdziennicka, A.; Jańczuk, B. Properties of some nonionic fluorocarbon surfactants and their mixtures with 
hydrocarbon ones. Adv. Colloid Interface Sci. 2021, 292, 102421. [CrossRef] 

Bhadra, P; Siu, S.W.I. Effect of concentration, chain length, hydrophobicity, and an external electric field on the growth of mixed 
alkanethiol self-assembled monolayers: A Molecular dynamics study. Langmuir 2021, 37, 1913-1924. [CrossRef] 

Hubbe, M.A.; McLean, D.S.; Stack, K.R.; Lu, X.; Strand, A.; Sundberg, A. Self-assembly of alkyl chains of fatty acids in papermaking 
systems: A review of related pitch issues, hydrophobic sizing, and pH effects. BioResources 2020, 15, 4591-4635. [CrossRef] 
Kloubek, J. Orientation of aliphatic hydrocarbons at the liquid-air interface. Colloids Surfaces 1990, 48, 323-335. [CrossRef] 
Fowkes, F.M. Attractive forces at interfaces. Ind. Eng. Chem. 1964, 56, 40-52. [CrossRef] 

Zdziennicka, A.; Szymczyk, K.; Krawczyk, J.; Jańczuk, B. Critical micelle concentration of some surfactants and thermodynamic 
parameters of their micellization. Fluid Phase Equilibria 2012, 322—323, 126-134. [CrossRef] 

De Boer, J.H. The Dynamic Character of Adsorption; Oxford University: Oxford, UK, 1953. 

Rathod, K.; Ahmed, H.; Gomte, S.S.; Chougule, S.; Prabakaran, A.; Dethe, M.R.; Patel, R.J.; Bharadwaj, D.; Alexander, A. Exploring 
the potential of anti-inflammatory activity of berberine chloride-loaded mesoporous silica nanoparticles in carrageenan-induced 
rat paw edema model. J. Solid State Chem. 2023, 317, 123639. [CrossRef] 

Sailor, G.U.; Ramani, V.D.; Shah, N.; Parmar, G.R.; Gohil, D.; Balaram, R.; Seth, A. Design Of Experiment Approach Based 
Formulation Optimization Of Berberine Loaded Solid Lipid Nanoparticle For Antihyperlipidemic Activity. Indian J. Pharm. Sci. 
2021, 83, 204-218. [CrossRef] 

Xiong, W.; Sang, W.; Linghu, K.G.; Cheang, W.S.; Li, J.; Hu, ҮЈ.; Yu, H.; Wang, Y.T. Dual-functional Brij-S20-modified nanocrystal 
formulation enhances the intestinal transport and oral bioavailability of berberine. Int. J. Nanomed. 2018, 13, 3781-3793. [CrossRef] 
Zdziennicka, A.; Szymczyk, K.; Jańczuk, B.; Longwic, R.; Sander, P. Adhesion of canola and diesel oils to some parts of diesel 
engine in the light of surface tension components and parameters of these substrates. Int. J. Adhes. Adhes. 2015, 60, 23-30. 
[CrossRef] 


Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual 
author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to 
people or property resulting from any ideas, methods, instructions or products referred to in the content. 


ZALACZNIK 4a 


MATERIAL UZUPEENIAJACY [SM4] 


M. Szaniawska, K. Szymczyk*, A. Zdziennicka, B. Jańczuk, Thermodynamic 
parameters of berberine with Kolliphor mixtures adsorption and micellization, 


Molecules, 2023, 28(7), 3115 


Supplementary Material 
Thermodynamic parameters of berberine with 


Kolliphor mixtures adsorption and micellization 


Magdalena Szaniawska, Katarzyna Szymczyk*, Anna Zdziennicka and Bronisław Jańczuk 


Department of Interfacial Phenomena, Institute of Chemical Sciences, Faculty of Chemistry, Maria Curie-Skłodowska University 
in Lublin, Maria Curie-Skłodowska Sq. 3, 20-031 Lublin, Poland 
*Correspondence: Katarzyna.szymczyk©mail.umcs.pl; Tel.: +48-81-537-55-38 


ty (mN/m) 


түү (mN/m) 


тү (mN/m) 


6 -5 -4 


logC 


(ELP + RH40) 


Figure 51. A plot of the surface tension (у;у) of M1 (a), M2 (b) and МЗ (c) aqueous solutions vs. the logarithm of 
ELP (logCg;p) RH40 (logCRn4o) and that of their concentration sum (log(CeLp + CRH4o) ) at the constant 
temperature equal 293 K (points 1 and curve 1’), 303 K (points 2 and curve 2’) and 313 K (point 3 and curve 3'). 
Points 1 — 3 correspond to the measured values, curves 1’, 2’ and 3’ to the values calculated from Eq. (35). 
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Figure S2. A plot of the constant y? (a), A, (b), A; (c), t4 (d) and t; (e) in Eq. (35) vs. the temperature (T). 
Curves 1 - 6 correspond to the aqueous solutions of ELP, RH40, M1, M2, binary mixture of ELP and 
КН40 at the mole fraction of ELP in the bulk phase equal to 0.8 [22] and M3 mixtures, respectively. 
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Figure S3. A plot of the mole fraction of ELP and Ber in M1 in the surface layer (curves 1 and 2) and in 
the bulk phase (curves 1’ and 2’) vs. the logarithm of ELP concentration (logCg;p) at the constant 
temperature equal 293 K (a), 303 K (b) and 313 K (c). Curves 1 and 1’ correspond to ELP, curves 2 and 
2’ correspond to Ber. 
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Figure 54. A plot of the mole fraction of RH40 and Ber in M2 in the surface layer (curves 1 and 2) and 
in the bulk phase (curves 1’ and 2’) vs. the logarithm of RH40 concentration (logCpy49) at the constant 
temperature equal 293 K (a), 303 K (b) and 313 K (c). Curves 1 and 1’ correspond to RH40, curves 2 


and 2’ correspond to Ber. 
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Figure 55. А plot of the mole fraction of ELP, RH40 and Ber in M3 in the surface layer (curves 1 -3) 
and in the bulk phase (curves 1 — 3’) vs. the logarithm of the sum of ELP and RH40 concentration 
log(Cg_p + Свндо) at the constant temperature equal 293 К (a), 303 К (b) and 313 К (c). Curves 1 and 1” 
correspond to ELP, curves 2 and 2’ correspond to RH40 and curves 3 and 3’ correspond to Ber. 
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Figure S6. A plot of the surface concentration (T) calculated from Eq. (5) vs. the logarithm of ELP 
concentration (logCgrp) at the constant temperature equal 293 К (a), 303 К (b) and 313 K (с). Curves 1 
and 2 correspond to the T values of ELP and Ber in M1 mixtures, curve 3 to their sum. Curves 1’ and 2’ 
correspond to the Г values for ELP and Ber in their individual aqueous solutions. 
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Figure S7. A plot of the surface concentration (T) calculated from Eq. (5) vs. the logarithm of RH40 
concentration (108Сднао) at the constant temperature equal 293 К (a), 303 К (b) and 313 K (c). Curves 1 
and 2 correspond to the T values of RH40 and Ber in M2 mixtures, curve 3 to their sum. Curves 1’ and 
2' correspond to the I values for RH40 and Ber in their individual aqueous solutions. 
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Figure S8. A plot of the surface concentration (Г) calculated from Eq. (5) vs. the logarithm of 
concentration (logC) at the constant temperature equal 293 K (a), 303 K (b) and 313 K (c). Curves 1-4 
correspond to the Г values for ELP, RH40 and Ber in M3 mixtures and their sum. Curves 1’, 2’ and 5 
correspond to the Г values for ELP and RH40 in the 0.8ELP + RH40 mixtures and their sum. Curves 1’, 
2’ and 3’ to theT values for ELP, RH40 and Ber in their individual aqueous solutions. 
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Figure S9. A plot of the Gibbs standard free energy of adsorption (AGQ,;) calculated from Eq. (39) vs. 
the logarithm of ELP concentration (logCg; р) at the constant temperature equal 293 K (curves 1 and 2), 
303 K (curves 1’ and 2’) and 313 K (curves 1” and 2"). Curves 1, 1' and 1" correspond to ELP, curves 2, 


2' and 2" to Ber in M1 mixtures. 
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Figure S10. A plot of the Gibbs standard free energy of adsorption (AG2,,) calculated from Eq. (39) vs. 
the logarithm of RH40 concentration (ІовСрнло) at the constant temperature equal 293 K (curves 1 and 
2), 303 K (curves 1’ and 2’) and 313 K (curves 1” and 2"). Curves 1, 1' and 1" correspond to RH40, 


curves 2, 2’ and 2” to Ber in M2 mixtures. 
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